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General Comments 

Chapters are divided into sections, for example §1.3, which are usu- 
ally partitioned into subsections, §1.3.2, which are, in turn, occasionally 
divided into further subsections, §1.3.2b. Tables and Figures are num- 
bered sequentially through each chapter as Fig. 1-1, Fig. 1-2,...; however, 
within a chapter the chapter number is suppressed, e.g., Fig. 1, Fig. 2,.... 
Equations are numbered sequentially through each sections, e.g., (1), 
(2),...; however, when referring to an equation in another section, the sec- 
tion number is included, e.g., (1.3-1), (1.3-2),.... Pages are numbered 
sequentially for each chapter. Most chapters and some sections start 
with a summary of the major points. 
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Summary: 

It is assumed that all students reading this material have had some 
course (e.g., the traditional semester of a junior-level physical chemistry 
course) presenting the basic elements of quantum mechanics with some 
treatment of the hydrogen atom, the harmonic oscillator, and angular 
momentum. This course will concentrate on the explanation of the struc- 
ture and reactivity of molecules using quantum mechanical ideas. The 
explanations will stress qualitative and semi-quantitative considerations 
with the emphasis on developing principles (based on quantum mechan- 
ics) that can be used to make reliable predictions on new systems 
(rather than merely rationalize known results). 

Chapter 1 is a review of materals that all students should have had 
previously, but with an emphasis on those points that will be important 
later in the course. 

The basic principles of quantum mechanics are summarized in § 1.1. 
A key idea here is that in the classical description of an atom, the elec- 
tron would collapse into the nucleus. The critical difference with the 
quantum description is that the kinetic energy is proportional to the 
average value of the square of the gradient of the wavefunction, T = 
Yi<\V<p\ 2 >. Consequently, for an electron sitting on the nucleus, the 
kinetic energy is infinite (since V<p is infinite). This forces the electron to 
remain distributed over a finite region surrounding the nucleus and 
prevents the collapse of the electron into the nucleus. Thus the quantum 
description is essential for stability of atoms. We will find in later 
chapters that modifications in the kinetic energy (due to superposition of 
orbitals) also plays the key role in the formation of chemical bonds. 
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Throughout this course we will be searching for qualitative ideas con- 
cerning the sizes and shapes of wavefunctions and for simple ways of 
predicting the energy ordering of the states of a system. A useful con- 
cept here is the nodal theorem described in §1.3. Basically, this theorem 
tells us that the ground state of a system is everywhere positive [no nodal 
planes (zeros) interior to the boundaries of the system]. 
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§1.1 BASIC PRINCIPLES OF QUANTUM MECHANICS 

In the following section we highlight the basic concepts of quantum 
mechanics relevant for this course. All of these ideas should be familiar 
to you; good references for reviewing these topics and for outside reading 
during the first part of Ch 120a are: 

I. N. Levine, Quantum Chemistry (Allyn and Bacon, Boston, 1974), 

Second Edition. 

H. Eyring, J. Walter, and G. W. Kimball (EWK), Quantum Chemistry 

(Wiley and Sons, New York, 1944). 
§1.1.1 The Need for Quantum Mechanics 

In order to see why quantum mechanics is so important to chemis- 
try, let's examine the classical mechanical description of the hydrogen 
atom 
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T + V 


where the kinetic 


energy is* 







(1) 



* Actually, the total kinetic energy of the hydrogen atom has two terms, 

T = L^ v * + L^ v * = -L-w + ^-(p.). . {3 ) 

However, considering the case where there is no net motion (i.e., no net inertia or 
momentum) leads to 

P P + 7> 8 = 
and hence (p p ) 2 = (p a ) 2 > so that (3) becomes 
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T = T™* - 5JT"' • < 2 > 

where m, v, and p are the mass, velocity, and momentum of the electron, 
and the potential energy is 

v = MfL. = _ *L . ( 6 ) 

r r ■' 

where q e = -e and qp = +e and the charge of the electron and proton and r 
is the distance between them. 

The ground state is when the system has its lowest possible total 
energy. Any other state (higher energy) is referred to as an excited state. 
Generally, systems in excited states will eventually decay to lower energy 
states,' and we will be interested in the stable (ground) states. The 
lowest kinetic energy occurs for p = 0, leading to T = 0, while the lowest 
potential energy occurs for r = 0, leading to V = -*>. Thus, in the classical 
description, the ground state of the hydrogen atom has the electron 
standing (or sitting) on the nucleus, leading to 

E = -» . 

Since the charges cancel and the atom has a radius of zero, these atoms 
would not combine to form molecules. Thus, in classical mechanics the 
atom is not stable! If classical mechanics provided the proper 

(4) 



(5) 
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Since rrip - 1836 rrig, then fi — 0.9995 m B , and for our purposes we can consider 
just the kinetic energy of the electron as in (2). 

' For systems containing charges, this is accompanied by emission of light. 
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description of atoms, and we would not be here pondering the universe. 

The solution to this problem is provided by quantum mechanics, as 
will be discussed below. Essentially the conclusion is that electrons must 
be described in terms of wavef unctions y>(r), where the shape of the 
wavefunction simultaneously determines both the kinetic energy and the 
potential energy. [In classical mechanics we can independently adjust r 
and p.] The result is that the state of the system with lowest potential 
energy (r = 0) has an infinite kinetic energy preventing the atoms from 
collapse. 

§1.1.2 INTERFERENCE AND DIFFRACTION OF LIGHT 

Before proceeding to a discussion of quantum mechanics, we will 
review some relevant features concerning the properties of light. 

The early controversy upon the nature of light between Newton, who 
considered light as corpuscles, and Huygens, who considered light as 
waves, was settled partly on the basis of the fact that (coherent) light 
waves interfere, a property difficult to explain except on the basis of 
waves. Basically, the idea is that 

(i) light is described by a -wavefunction 

iKx.t) 
that depends upon x and t, for example, 



ip{x,t) = cos 



-r—x — 2-nvt 

A 



(where \ is the wavelength and v is the frequency), 

(ii) detection of the light is proportional to the square of the wavefunc- 
tion (called the intensity ) averaged over a time long compared with 
the frequency 
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I(x) = <[i,(x,t)f> t 

where the brackets indicate an average and the subscript t indicates 
that the average is over t, 

(iii) superimposition of two wavefunctions leads to a new wavefunction 

where the amplitudes add, and 
(iv) the intensity for two superimposed waves is [from (ii)] 
I™ w (x) = <[i/f d + ^| w ]> f 

= <[ff d ] 2 > t + <|>! w ] 2 > f + 2 <i f w if t d > t 
= if* + if* + /« , 

where the /f and /§f* are the intensities of the component waves and 
lint is a new interference term present only when the component 
waves are present simultaneously. 

The interference term in (iv) may be nonzero and can lead to complete 
cancellation of the other terms. Particularly impressive interference 
phenomena are the diffraction effects found for such uniformly spaced 
scatterers as diffraction gratings as illustrated in Fig. 1. 



(a) 
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wavelength X 



(b) 
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T 
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-> 



screen to 
observe intensities 




Figure 1-1. 

With a set of uniformly spaced scatterers, the observed intensities are 
sharp spikes at particular angles tf n , where* sintf n = n — From measure- 
ments of sintfobs one can calculate — - Therefore, knowing the X of light, 

we can determine the spacing a or vice versa, knowledge of the spacing a 
can be used with sin tf obs to determine A. A comparison of the observed 
intensity with that expected if light did not interfere is given in Fig. 2. 



* To have maxima for $ n ^ 0, we must have nX<a; the wavelength must be small- 
er that the spacing of the scattering. 
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Figure 2-2. 



A particularly exciting application of these ideas occurred around 

1912. By that time a number of scientists believed that x-rays were elec- 
ts 
tromagnetic waves like light but with very short wavelengths, \ ~ 1 A. If 

so, they should exhibit diffraction, if a grating could be found having 
equally spaced scatterers with uniform separations of ~1 A . In addition, 
by 1912 a number of scientists were convinced that atoms do exist 
(rather than being just theoretical constructs) and that crystals might 
consist of uniformly spaced atoms having separations of a few angstroms. 
F. von Lave, an expert on diffraction theory, suggested the experiment of 
exposing a crystal to a beam of x-rays and looking for diffraction spikes. 
After a couple of years of work the experiments were successful, proving 
both the wave nature of x-rays and the existence of ordered atoms in cry- 
stals. Since then, such x-ray diffraction studies have led to enormous 
advances in our atomic-level understanding of matter. 
§1.1.3 Electrons 

The critical experiment* establishing the wave nature of the elec- 
trons is that a crystal diffracts a beam of electrons in exactly the same 
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was as it diffracts a beam of x-rays, as illustrated in Fig. 3. Thus, elec- 
trons must be described as waves. 



x-rays 



(a) 



electrons 



xtal 



(b) 



xtal 



Figure 1-3. 

Indeed, from these observations one can determine the relation between 
the wave properties and particle properties of the electrons, namely, 



P = F' 



(!') 



where p = ^/2m.E is the momentum of the electron and E is the kinetic 
energy of the electron, h is a constant (Planck's), and X is the wavelength 
of the electrons (obtained from the spacing of the diffraction peaks). 

Based on this and other experiments, we now know that electrons 
should be described as wavef unctions, 

where the probability, P, of observing the electrons at some point x is 
proportional to the square of the wavefunction, 

P(x)=<[^(x,t)f> t . 



* This experiment was carried out in 1928 as a test of the ideas arising from the 
theorists developing quantum mechanics. Earlier experiments had, in fact, ob- 
served what is now recognized as diffraction; however, the experiments were not 
properly interpreted. 
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The consequences of this will be pondered in the next section. 

§1.1.4 The Sehrodinger Equation 

In the 1920's, a number of experiments, such as electron diffraction, 
showed that matter exhibits interference phenomena just as does light. 
This led to the idea that matter, like light, should be described by an 
amplitude function, 

called a wavefunction, such that superposition of two systems fa and fa 
leads to superposition of the amplitudes 

it = fa + fa ■ (1) 

but such that the probability of finding the total system with particular 
coordinates r and t is given by the (absolute) square of this amplitude. 

P(r,t) = |#r,OI 2 = 1<*(r,t)far,t) . (2) 

Combining (1) and (2) leads to 

P(t ,t) -fa *fa +■ fa *fa + fa *fa + fa *fa , 

and hence interference effects— as observed. 

Putting these ideas together leads to the basic postulate of quantum 
mechanics: Every physical system is described in terms of a wave/unc- 
tion f coritairiirig all observable information about the system. This 
wavefunction is probability amplitude, meaning that a superposition of 
states of the system leads to a superposition of the amplitudes, 

* = Eft • (3) 

i 

As part of this basic postulate, we assume that if fa and fa are two 
acceptable wavefunctions of a system, then 
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f = c 1 i/ l + c z f z . (4) 

(where c l and c 2 are constants) is also an acceptable wavef unction (this is 
called the principle of superposition). 

As part of this basic postulate, the probability of the system having 
particular coordinates at a particular time is taken as the absolute 
square of the wavefunction |^| 2 as in (2). Since the total probability of 
the system being somewhere is 1, we have the normalization condition: 

1 = f P(r,t)dT - ff*i/dr = <f|^> (5) 

(dr is the incremental volume element and this integration goes over all 
regions of space). The normalization of the wavefunction ip is indepen- 
dent of time (being always unity). 

The above postulate implies that anything we can know about the sys- 
tem must be extracted from the wavefunction. Thus the wavefunction at 
some future time, 1 + 6t , is completely determined by the form of the 
wavefunction at time t . In other words, there must be some rule or for- 
mula relating i/(r,t + 6t) = p>i(r) to i/{r,t ) = <pq{t). Such an association 
of functions is generally called a transformation and is denoted as 

Pi(r) = Aio^oW , 

where 4 10 is referred to as the operator effecting the particular transfor- 
mation from state tp to state ^. Similarly, the time derivative of the 

wavefunction at time t Q , — ^ must be determined by the form of the 

at 

wavefunction at time 1 , and hence we can write 



d± 



= S^(r,t Q ) , (6) 



dt 
where B is called the time evolution operator. For convenience we 
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replace B with H = ihB, where H is referred to as the Hamiltonian. 
Thus (6) becomes 



««--/» 



h^ . 



(?) 



which is known as the time-dependent Schrbdinger equation. In (7), i = 
V-i, ft is Planck's constant* (1.054589 10 -27 erg sec), and H has the 
dimensions of energy. Since (7) must also apply to any superposition of 
wavefunctions (3), H must be a Linear operator J 

H(f 1 +f 2 ) = Hift + Hfr . (8) 

We find that the operator H depends upon the nature of the system and 
that it is, in general, a function of both position r and time t. 

If the Hamiltonian H is independent of time, then the solutions of (7) 



have the form 



where 



y(r,t) = tp{r)T(t) 



dT 
dt 



= -iuT{t) 



(9) 



(10) 



* Actually, the original Planck's constant 6. is 

h - Z-nh. . 

However, we will use only hand refer to it as Planck's constant. 

* From (4) 



ih w i, = Cl 






ih di/ * 



dt 



= CtHiti + c z His z . 



where (7) was applied to t{/ 1 and i/ 2 , respectively. Applying (7) directly to ^ leads 
to 

a 



and hence (8). 
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and 

H{t)<p(t) = &>p<r) . (11) 

Equation (10) has the solution 

T(t) = e~ iut , . (12) 

so that (9) becomes 

#r.O = <p{t)2^ , (13) 

where p(r) is yet to be determined from (11). 

At this point we recall the quantum mechanical interpretation of two 
experiments. In the photoelectric experiment, light behaves as a stream 
of particles (called photons), each having a quantum of energy 

E = hj , (14) 

where 

u = 2nv = 2nc/X 

is the angular frequency of the light. This suggests that the hu in (11) be 
considered as the energy in (14). In electron diffraction, the diffraction 
pattern for electrons with momentum p and energy E is equivalent to the 
diffraction pattern for light with wave vector 

k = 2n/\ (15) 

and angular frequency w, where u is given by (12) and k is given by 

p = ** . (16) 

Thus we postulate that the energy and frequency are always related by 
(14), leading to 

f(r,t) = <p(r) b-**'* (17) 
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H(rMr) = E<p(r) 



(18) 



The latter equation is known as the tome-independent Schrbdinger equa- 
tion and is the fundamental equation determining chemical bonding. 

§1.1.5 The Form of the Hamilton! an 

In (18) we see that there is a relationship between the operator H 
and the total energy of the system, E. In classical mechanics, the opera- 
tor associated with the total energy of the system is the Hamiltonian, H cl , 
which is given by 

(for nondissipative systems), where T cl and V 01 are the kinetic and poten- 
tial energies. We will postulate that there are quantum mechanical 
operators f and V, corresponding to the classical quantities T and V, such 
that the quantum mechanical operator H is given by 

H = f + V , 

and we will refer to H as the Hamiltonian operator. For a system in 
which the classical potential V* 1 is velocity-independent [that is, a func- 
tion of the coordinates of the particle only], we will postulate that the 
quantum mechanical operator corresponding to 7(r) is just the classical 
function 

V(r) = V{r) cl . (19) 

Thus for the hydrogen atom, 

V(r) = - 



e 2 



For a particle moving in a potential V(r) ci , the kinetic energy (classically) 

is 
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2m ^ 



(20) 



where p = my is the momentum of the particle. We will postulate that the 
quantum mechanical Hamiltonian operator corresponding to (20) is 



2m * 



(21) 



where p is the quantum mechanical operator corresponding to the 
momentum. 

Now we need the form of the quantum mechanical momentum opera- 
tor, p. A plane wave of wave vector k and angular frequency a has the 
form 

<p{x,t) = b***-** , (22) 



and hence the wave vector is given by 



k = — 

9 



1 dtp 



i dx 

From the diffraction experiments it was found that p = fik, leading to 



1 
P = — 



ft dp 



i dx 



Thus we postulate that the momentum operator p z corresponding to 
momentum in the x direction is given by 



Px - 



i dx 



(23a) 



and similarly for the other directions 



Py = 



i By 



(23b) 



- = 5_ JL 
Pz i dz 



(23c) 
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Just as the classical momentum is a vector quantity, the three quantities 
in (23) are considered as the three components of a vector operator 



p = -v 



(24) 



where V is the gradient operator. 

Now we construct the kinetic energy operator. Since 



y2 = V ■ V = 



a 2 a 2 a J 



dx' 



dy 



dz' 



(25) 



s obtain 












P 2 = P P = 










= - 


-&\ 


7* 


- - 


ff 


a 2 a 2 a 2 

dx 2 dy z dz 2 



and hence 



T = 



—P 2 = 
2m F 



2m 



■V 2 



(26) 



From (18), (19), and (26) we obtain 



H = 



& nZ 



2m 



V 2 + V(r) 



(27) 



as the explicit form of the Hamiltonian for a particle of mass m moving in 
a potential 7(r). 

Basically, the Schrbdinger equation (8) arises from considering the 
time evolution of a system, and the Hamiltonian H describes how the sys- 
tem changes with time. If we change the system, say, by applying an 
electric or magnetic field, this change is manifested by a change in the 
Hamiltonian H. Such changes in H lead to changes in if. With suitably 
ingenious experiments, it is often possible to determine something about 
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how iff changes in response to the field and thereby something about the 
form of i> before changing H . In this way we can determine various pro- 
perties of i/. Ultimately each physical property can be related somehow 
to some type of change in the Hamiltonian of the system and hence to 
some (Hermitian) operator, 

Op = A# . 

§1.1.6 More on the Schrodinger Equation 

§ 1.1.6a The Hilbert Space 

Given any two functions fa and fa, "we can generate from (3) an 
infinite number of wavefunctions 

^ = difi + c 2 fa 

by using various Ci and C 2 . In addition, there is an infinite number of 
choices for the functions fa and fa. Even so, the postulates of quantum 
mechanics lead to constraints on the functions, and hence we need not 
consider every wavef unction. For example, from (5) we need consider 
only wavefunctions for which the integral of the square of the wavefunc- 
tion is unity <if/\i/> = 1. Of course, given some wavef unction i/ with 

<^ |^> = J d t \fa z - a 

with finite (nonzero) a, we can always define a new function* 

^ = ^/Va 

that is normalized, i.e., 

<f\ip> = 1 . 

On the other hand, we need not consider any wavefunctions iff for which 



* Note that a can never be negative. 
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the integral f dr \i(/\ z does not converge. That is, we need deal only with 
square-integrable functions. The set of all possible such functions (satis- 
fying whatever boundary conditions are being imposed) is referred to as 
the Hilbert space (for systems having this particular set of boundary con- 
ditions). Thus the Hilbert space is merely the collection of all possible 
wave/unctions for our system. 

§ 1.1.6b Hermitian Operators 

In Appendix A we consider the implications of requiring that the 
norm of the wavefunction be unity, 

<f\i/> = 1 

and hence independent of time for any superposition of wavef unctions, 

it = ifti+ltj ■ 

The conclusion is that for all possible functions ^ and pj the Hamiltonian 
operator if must satisfy the condition 

fdrfHiti)*^ = fdrit^CHitj) , 
which we denote as 

<\Hit i )\it j > = <Viffr|ty>„ 
Such an operator is called Hermitian. 
The expectation value 

<i/ 1 f> 

of a Hermitian operator is always real (see Appendix A). Hence the 
energy 

= <j,\H\j,> 

<i/\i/> 
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in the Schrbdinger equation must be real. 

In Appendix A we show that the momentum operator, p = (fi/i)V, and 
the kinetic energy operator, f = (l/2m)p 2 , are hermitian. Similarly, any 
function of coordinates, V(r), is hermitian, so that the Hamiltonian in (27) 
is also hermitian. 

§1.1.7 Analysis of Kinetic Energy and Potential Energy 

In the above sections we have established the Schrbdinger equation 

H<p = Ey , (18) 

where 

H = f + V 



f = - -^-V 2 
2m, 



V = V(r) . 
Multiplying both sides of (18) by <p* and integrating leads to 

<<pfH\<p> = E<cp\<p> = E , 
where 

<<p(H\tp> = fdr(p*H(p 
and 

<<p\?> = 1 . 
Defining the numbers T and V as 



T = <<p\f \tp> - f d t <p*{r) 



2m 



72 



9(r) (28) 



V = <<p\ V\<p> = f d r<p*(r)V(r)<p{r) , (29) 
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we see that the total quantum mechanical energy E can be written as a 
sum of quantities 

E = T + V (30) 

interpreted as a kinetic energy (T) and potential energy (V). 
The quantity (29) can be rewritten as 

V = fdrp(r)V(r) , 

where 

p{r) = <p*{r)<p{r) 

is the probability of finding the system in the volume element dr near 
configuration r. Thus V corresponds to the average of the classical poten- 
tial energy, weighted by the probability of the electron being at any par- 
ticular position. 

As written in (6), T does not seem to bear much relation to the classi- 
cal kinetic energy. However, in Appendix B we show that 

<(p\-^\<p> = <|Vp| 2 > , (31) 

so that (28) becomes 

" 2 



f= ^<w> 


2m < 


-r-Vc0 


7, we see that 

f = 1 

1 Onrt 


„ <\P9\ 2 > . 



(32) 



which can be compared with the classical kinetic energy, T l = 



1 ~ —P 2 , 
2m F ' 



suggesting that <\p<p\ z > corresponds to the square of the classical 
momentum. Throughout this course we will find (32) to be a useful way to 
think about kinetic energy. This expression says that big -gradients or 
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slopes lead to large kinetic energy, and hence the best kinetic energy 
occurs for the smoothest functions. Thus, comparing the wavefunctions 
in Fig. 4 (all normalized), we see immediately that <p c has the highest f, 
while (p a has the lowest.* 





+< 



Figure 1-4. 



*, 



The essential difference between classical mechanics and quantum 
mechanics is that in classical mechanics the kinetic energy and the 
potential energy are independent (one is determined by momentum, the 
other by position), whereas in quantum mechanics f and V are simul- 
taneously determined by the ivavefunction, with the kinetic energy pro- 
portional to the average square of the gradient of the amplitude function! 
It is the balance of trying to find a wavefunction leading to both the 
lowest T and the lowest V that is responsible for the stability of quantum 
mechanical atoms. 



* Of course, (32) assumes that <(p | <p> = 1. 
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§1.2 The Ground State of Hydrogen Atom 

In this section we consider the ground state of the hydrogen atom, 
that is, an electron with mass m and charge -e interacting with a nucleus 
of infinite mass and charge +Ze. Classically, the energy is given by 

2 Ze 2 



E = E-- 



2m r 

where r is the distance of the electron from the nucleus. Thus the 
ground state (lowest energy) is for r = and p = 0, leading to E = -«. 
That is, the classical H atom collapses to a point. 

Quantum mechanically, the Hamiltonian is 

and the energy is obtained by solving the Schrbdinger equation, 

H<p{r) = Etp(r) . (2) 

We will find that the quantum mechanical form of the kinetic energy 
keeps the electron from collapsing into the nucleus. 

In these sections we will obtain the wavefunction tp(r) for the ground 
state of H atom. The result is that 

<p{r,^,<p) = N e~< T , 

where 

N Q = Vp/^r 

a 



E = -^Z 2 £- 
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a = ft 2 / me 2 . 

In §1.2.3 we define atomic units where 

. ft' = 1 , ra = i , |e | = 1 . 

In these units, the unit of length is 

1 bohr = la = ^/me 2 = 0.529177 A , 
and the unit of energy is 

lhartree = l/i = e 2 /ao = me 4 / ft 2 - 27.2116 eV = 627.510 kcal/mol. 
In these units, the Hamiltonian for H atom becomes 

H = -&V 2 - ^-, 

T 

the energy becomes 

E = -%Z 2 , 
and the scale parameter becomes 

t = Z . 

Before going into the details of the wavef unctions of the hydrogen 
atom, we will consider why such an atom can exist. 

§1.2.1 Atoms Exist! 

In §1.1.1, we found that the classical description of the atom leads to 
collapse, 

T ~ 7T-P 2 -* as p -> 

2771 



e 2 
V = > — oo as r -» 



Therefore, the lowest energy state is for the electron sitting on the 
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nucleus. Since the charges cancel, this is like not having an atom. Now 
we will look at this problem with quantum mechanics. A major difference 
in quantum mechanics is that both f and V are determined by the same 
quantity, the wavefunction, whereas in classical mechanics, T and V 
involved independent quantities p and r. Thus, 



< 


P 2 


V - <* 


B 


< 


r 



> 
> 



(3) 



tf 



t - £^<m*> . 



Consider now the sequence of similar wavef unctions in Fig. 5.* 



* To be specific, consider the normalized function 



$ = N e 



-r/H 



(4) 



where** N = 1/ v TrR 3 . Thus, if R is very large, # decreases very slowly with r, 
leading to the delocalized function a in Fig. 5a, while with very small R , <p de- 
creases rapidly to zero for small r, leading to the localized function in Fig. 5c. 



fT oo 



<$ 2 > = y gin tfdtf y cttf J $2 r 2 dr = 4ttN 2 J r z dr e ~ 2 ( r/ *) 



- 4ttN 2 



R_ 

2 



f p z dp2'P = nN 2 R 3 = 1 



where we used J e p p m dp = m\ 



z§- 



(a) 



■(b) 



R 

de localized 
_ large R 
V slightly negative 




(c) 




localized 
_ small R 
V very negative 

T very positive 



T slightly positive 

Figure 1-5. 

Clearly, V becomes more and more negative (lower energy) as the elec- 
tron is localized closer and closer to the nucleus (just as in classical 
mechanics), and in the limit the wavefunction leading to the best V is 
localized at the nucleus (R = 0). However, this localization of the elec- 
tron near the nucleus now leads to a very large and positive T. Since V 
and f have opposite effects as the electron is concentrated near the 
nucleus, we need to be a little more quantitative in the analysis. 

First we define an average radius R as 



< 
< 



> 1 



leading to 



V = - 



R 



(5) 



Consider now some wavef unction, say b, in Figure 5 as the reference 
wavefunction (with R = 1 in some units) and let 

?! and 7*! 

be the energies for this wavefunction. Using this reference point, we will 
examine how V and T change as the wavefunction is squeezed or 
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expanded.* 

From (5) we see that 



Vt> = 



< 



> 



R 



(6) 



In order to see how kinetic energy changes, note that each term has the 
form 



< 

<l 



dx 



> 

J > 



so 



that* 



7V> = 



1_ 
R 



>2 



(7) 



Thus, f becomes small for delocalized smooth functions (large R), and T 
becomes large (and positive) for localized functions (small R). From (6) 
and (7) we see that 



Ik 



T l 



R V x 

(note that f is always positive and V is always negative) . 
Consider first the case as R -* », then from (6) and (7) 

T -> 
V -» 

and 

E - T + V -» , 



(a) 



* The technical term is scaled. 



1 



t Basically, the gradient is proportional to ^-and hence the gradient squared is 
f^2 R 



proportional to 



R 
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as expected. For sufficiently large R (that is, R » | 7V VM), we see from 
(8) that 



and hence the total energy 



V\ » \T 



E - T + V 



must be negative. However, for very small R (that is, R « | TV Vi\), we 
see from (8) that 

|f | » |7| , 

and hence the total energy must be positive. Thus, the energy of the 
wavefunctions in Fig. 5 must behave as in Fig. 6 as a function of R (i.e., as 
a function of the size of the wave function). That is, the lowest energy 
(corresponding to the ground state of the atom) occurs at a finite size, 
R = Ropt- l n quantum mechanics the hydrogen atom is stable! 



E = 




>R 



Figure 1-6. 

In the above example, we considered just the stretching and 
compression of the one function considered in Fig. 5. However, the same 
result is obtained independent of the shape (namely, the optimum energy 
occurs for finite R), and hence trying all possible shapes we will eventu- 
ally find the optimum wavef unction and its optimum R. This optimum 
wavefunction is discussed in the next section. 
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Summarizing the above discussion, we find that: The potential 
energy wants the wavefunction to be localized at the nucleus. Thus, 
starting with a delocalized wavefunction (Fig. 5a), the total energy drops 
as the wavefunction is localized closer to the nucleus. This localization 
that aids the potential energy leads concomitantly to a more repulsive 
kinetic energy; however, for sufficiently diffuse wavefunctions, potential 
energy always wins.* On the other hand, the kinetic energy increases 
quadratically as the wavefunction is compressed, while the potential 
energy only drops linearly, so that eventually the increase in kinetic 
energy will prevent any further contraction of the wavefunction. For the 
optimum wavefunction there is a balance in these potential energy and 
kinetic energy terms. One might say that kinetic energy provides a pres- 
sure that keeps the atom from collapsing. 

§1.2.2 The Ground State Wavefunction 

Now we wish to obtain the wavefunction p(r) of the ground state of H 
atom, 

H<p(r) = E<p{r) , (2) 

where 

*- -£-*-¥-■■ (1) 

The Hamiltonian in (1) is independent of orientation of the atom in 
space, and hence the eigenf unctions will have the form 

fir)z(*#) , (9) 

where f is a function of r only and Z(fl,p) is a function of angular coordi- 



* We are assuming here Coulombic attractions. 
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nates only. Since kinetic energy favors having smooth wavef unctions, the 
ground state wavef unction should be as devoid of wiggles as possible. 
Thus we will take Z{-d-,<p) as a constant leading to 

"£: v2/(r) ~T i/(r) = Ef{r) (10) 

for the Schrbdinger equation. 

There are straightforward mathematical techniques for solving (9); 
for example, see Appendix D. Here we will use a physically oriented 
approach to examine some features of the solutions. At r = », the poten- 
tial in (1) is zero; thus the bound states of (2) have negative energy, 

E < . 

Now consider a very large r so that the Coulomb term is negligible, 

f L « 1*1 • (II) 

In this case the Schrbdinger equation reduces to 

ff -V 2 /(r) = Ef(r) 



2m 

or 

V 2 /(r) = +ff(r) , (12) 

where 

(note that E is negative and hence f is real). Consider a point along the 
positive x axis. Since r is very large, (df/dy) w and (df / dz) a 0. 
Thus (12) becomes 



dx z 



- yl 



+ t5t= cf • 
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Consequently, 

(/ = e+* z is also a solution, but this function is not nor maliz able). Since 
f is spherically symmetric, the wavef unction at very large r is of the form 

f(r) = e~< r . (14) 



Energy 



Region 

satisfying 

(11) 



energy= 



X 



energy = -E 



classical turning point 




Region 
-^satisfying 
(16) 



Figure 1-7. 

In Appendix D we show that the wavefunction (14) is an eigenfunction 
of (10) for all r if f is chosen so that 

(15) 




where 



From (13) and (15) we have 



F - - *$ t-2 - 



a n = & 2 /me 2 



(16) 



2tti clq 



2„2 



1_ ZH 
2 a. 



and hence 
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E = -%Z* 


do 





Normalizing the wavefunction ( 14) leads to (see Appendix D) 



^(r.tf.p) = N e Zr/a ° , 



where 



(17a) 



(17b) 



AT = VzVJra£ 



(18) 



This wavefunction is plotted in Fig. 8. 



*(0,0,s) 







>z 



z=-a o z=0 z = a„ 
(a) Line plot of the H atom orbital. 



♦ = 0. 05 



* = 0. 10 




z = 

(b) Contour plot of the H atom orbital. 
Adjacent contours differ by 0. 05 a. u. 



Figure 1-8. 



The average radius of the wavefunction (16) is 



■=- = <i(/ 
R < 



> 7 

y> = -£- 

> oo 



so that 



R = 



Z ' 



(19) 



where a = ft/me 2 is referred to as the Bohr radius (or, more simply, the 
Bohr) in honor of Niels Bohr. Substituting into (15) leads to 



*■-»* 



(20) 
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which can be compared with 

< 

V = <f 



Thus 



and 



Ze z \> 
r > 



E = ViV 



Ze< 



R 



T = E -V = -%V = % 



Zb> 



R 



Equation (20) provides an easy way to remember the proper energy 
expression (it is just half the total potential energy). 

For the hydrogen atom (Z = 1), the above equations become 

,2 



E = -Vi 



a, 



(21) 



V = - 



a 



(22) 



T = + & — = + 



ff 1 



0L f 



2m 



d§ 



where 



<l/r> = 



a ( 



< [vp? > = rr • 



„2 
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Energy 



V = 




Average Distance of Electron from Nucleus 

(|z|) =a ) 
Classical Turning Point , i , „ . 

(|z| = 2a ) 



distance along z axis 



energy of ground state 

(E = 4e 2 /a ) 



Potential energy (V = e /r) 



Figure 1-9. 

In Fig. 9 we show the potential energy as a function of distance. Imagin- 
ing a classical particle with this same energy E moving in the potential 
V(r) = - (e 2 /r), we would find the particle bouncing back and forth from z 

2 2 Z 

= -2a to z = +2a . The kinetic energy isT = E-V = E + — = - ■§• — + — . 

7* cCLq T 

For r > 2a , we would have T < 0, but this is not possible classically since 
T = Yittlv 2 must be positive. Thus the classical limit for the motion of the 
electron is |z| = 2a , at which point the velocity has reduced to zero. In 
the quantum description there is a finite probability (but not large) of the 
electron being farther than 2a from the nucleus. 

Note that the wavefunction (17a) is positive for all finite x, y, z. Since 
it is never zero for finite distances, we say that the wavefunction is node- 
less. 

§1.2.3 Atomic Units 

— Ctn 

As mentioned above, the "size" of the atom is R - -^-; thus a natural 



unit of length for atomic problems is the bohr radius 



a = 



me' 



(23) 
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The quantity (e 2 /ct = meVtf) in the energy expression (17b) has 
units of energy. It is referred to as the hartree (in honor of D. R. Hartree 
who first suggested the atomic system of units in 1927)* and is denoted 
as h , 

e 2 me 4 



a tf ' 

so that the energy of the hydrogen atom (17b) becomes 

E = -&Z 2 /i . (25) 

Throughout this course will encounter quantities such as a = f?/me z 
corresponding to a length, h = e 2 /a = me 4 / ft 2 corresponding to an 
energy, and dimensionless quantities such as 

e 2 1 

a = 5T = 137.03604 ^ 

(the fine structure constant). 

We will find it convenient to use a particular set of units, called 
atomic units,** where 

h = i 

m. = 1 (27) 

|«| = l . 

We use these units because they simplify many of the equations of quan- 
tum mechanics and give a reasonable order of magnitude (near unity) for 
the properties of molecular systems. Some useful conversion constants 



* D. R Hartree, Proc. Cambridge Phil. Soc, 24, 89 (1927). 

** Sometimes called hartree atomic units tO distinguish from other occasionally 
used atomic units; see Appendix E. 
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are included in Appendix F. 

With atomic units, the Hamiltonian for the hydrogen atom becomes 

H = -^V 2 -f (r) 

T 

and the ground state wavef unction become 

i/{r,V,<p) = N e- Zr , (17) 

where 

N = VFT^ (28) 

and the ground state energy becomes 

E = -YzZ* . (17') 

§ 1 . 2 . 3 a Conversion Factors 

In order to compare the results of quantum mechanics (expressed in 
atomic units) with those of experiment (expressed in cgs units), it is 
necessary to become facile at converting between those units. To do this 
most simply it is sufficient to remember a few basic conversions (e.g., 
energy and distance) and to rewrite the expressions involving other quan- 
tities in terms of these basic units before converting. Thus, from (17) the 
ionization potential of the ground state of hydrogen atom is 



IPh = 3* 



a 



In atomic units this quantity is 

IP H = 0.5 h 
and experimentally it is known that* 



* The numbers in parentheses indicate the estimated limits of error in the last di- 
git quoted, the standard deviation. 
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IP H = 13.605805(37) eV ** 
Thus, 1 " 

o 2 

lh Q = — = 27.21161 eV . (29) 

Similarly, using the known values of fi, m, and e 

b. = 1.0545887(56) x 10 -27 gm cm 2 sec -1 

m e = 9.109534(47) x 10" 2a gm (30) 

e = 4.803242(14) x 10~ 10 grn^cm 3/2 sec -1 , 
we find that the atomic unit of length is 

U = -^— = 0.52917704(44) A . (31) 

me* 

Another useful relation is the fine structure constant 

ff_ = 1 

a he 137.03604(11) ' 

a dimensionless constant. In atomic units, e = 1 and &= 1, and hence the 
speed of light is 

c = 137.03604 a.u. (32) 

in atomic units. From the cgs value of c, 

c = 2.99792458(1) x 10 10 cm/ sec , 

we find the atomic unit of time (denoted r ) 

t = 2.41888 x 10 -17 sec . (33) 



** This is corrected to correspond to a nucleus of infinite mass. 

t An underline indicates that there is uncertainty in the underlined number. 
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The biggest disadvantage in using atomic units is that the various 
quantities such as ft, e, and m e will be missing from the equations, making 
it difficult to convert to cgs units. The best way to convert is to rewrite 
the quantity of interest in terms of energy, length, and velocity quantities 
and then to use the above conversions. For example, what if we want to 



evaluate the Coulomb interaction between Na + and Cl~ at 10 A? We con- 
vert R to bohr 

R = — = 18.897 a a . 

0.529177 ° 

Since |e| = 1, the energy of interaction in atomic units is 

E = - 4" = "0.052918 h Q . 

IX 

We can now convert back to electron volts, 

E - -0.052917x27.2116 = -1.43998 eV . 
In general, then 

*(eV) = ^fp-. (34) 

It would be instructive to calculate other quantities in terms of 
atomic units, e.g., what is the average momentum of the electron in the 
ground state of H? In Bohr's model, how long does it take the electron to 
orbit the atom? 
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§1.3 THE NODAL THEOREM 

Even without carrying out detailed calculations of the eigenfunctions 
for a system, it is often possible to make some general conclusions con- 
cerning the ordering of the states by considering their nodal structures 
(i.e., the loci of points for which f = 0). Here we will develop some of the 
general considerations for such analyses. 

First we consider a one- dimensional, one-particle system with Ham- 
iltonian 

H = -w£r+v(*). CD 

where M = m/6? and V(x) is some function of x depending only upon the 
spatial coordinates.* If this potential contains bound states, then we can 
prove the nodal theorem: 

i) The ground state wave function does not change sign (i.e., has no 
nodes); 

ii) The bound state with n sign changes (n nodes) has a lower energy 
than the state with n + 1 sign changes (n + 1 nodes). 

That is, 

Eq < E n 

(2) 
E n < E n +i , 

where n is the number of nodal points (internal to the boundaries). In 
the case of a sufficiently singular potential, some inequalities in (2) may 
be equalities. 



* [V(x) is independent of momentum and spin and is not an integral operator]. 
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The nodal theorem is proved in Appendix C. Here we will provide 
some intuitive reasoning concerning this theorem. 

§1.3.1 The Ground State is Nodeless 

Consider first f lt the lowest eigenstate of the Hamiltonian 

Hf, = e^ 2 
having one sign change as in Fig. 10a. 





Figure 1-10. Illustration for the Nodal Theorem. 

From Tp x we can form a wavef unction <p = 1^ | (not necessarily an eigen- 
function of H) that has no sign changes, as in Fig. 10b. Since it l is nor- 
malized, then cp is also normalized, 

<<Po\<Po> = fdc\M x )\ 2 = f&c\1<i(x) 
and the energy of <p Q becomes 



12 _ 



1 , 



£o = <<Pq\H\<Pq> - <( Po 



£r£r\*> + <*>\y\*>> 



The potential energy of <p is the same as that of ^ lf 

<tp \V\^ > = fdx V(x)\<p (x)\ 2 = fdx V{x)\^{x)\ z = <^|7|Vi> ■ 0) 
From Appendix B the kinetic energy of ^ can be expressed as 
1 d 2 



< 
< 



2771 dx z 



*l* 5S-<I **>!■> 
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1 + ~ 
= / dx 



2M 



dx 



But 



toft** 



dx 



= %fdx 



dy 
dx 



since the integrands are equal except at one point (r ). Thus the kinetic 
energies of <p and fi are equal, 



< 
< 



d< 



ZM dx 2 



> < 

fi> = <<Po 

> < 



ZM dx 2 



> 

<Po> 



(4) 



and consequently from (3) and (4) the total energies of ^ and <p are 
equal, 

That is, given any eigenfunction fx of H that changes sign, we can con- 
struct a function <p which does not change sign and yet has the same 
energy. 

Now consider a new function ^ identical to <p Q except that it is 
smoothed in the region very close to the position of the node c. If the 
potential is not singular at this point, the function ^ can be chosen to 
have the same potential energy (and normalization) as <p , 



<Po\H\<Po> = <<Po\V\<p6> 
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(b) 



Figure 1-11. 



However, since 



dx 



« 



dx 



in the region near c, <p Q ' will have a smaller kinetic energy than <p , 



2M < 


d<p 
dx 



2 > 
> < 



1 K < 

2m < 


dcp 
dx 



(5) 



Consequently, the energy of 9? is lower than that of <p Q 



£n < £n = .EN 



The best (i.e., lowest energy) nodeless (i.e., non-negative) wavefunction 
has an energy, E , no higher than ~s and hence 



En < E, 



(6) 



Similar arguments can be used to derive the other relations. Thus, for a 
general potential we expect the bound solutions to increase in the 
number of nodes as E increases, as in Fig. 12. 
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Figure 1-12. Illustration of the nodal patterns of successive states of a 
general (one-dimensional) potential. 

§1.3.2 Multi dimensions 

In two dimensions, a wavefunction that changes sign will have a line 
of points with f = (a nodal line), and for three dimensions there will be 
a surface of points with if/ = (a nodal surface). Just as in one dimension, 
the ground state will always be nodeless. However, for multidimensions 
one can no longer use the nodal theorem to order all states. Thus, in two 
dimensions we can construct three orthogonal wavefunctions (all orthogo- 
nal to the ground nodeless state), each with one nodal surface, as illus- 
trated in Fig. 13.* 



* We will, in this section, use the same notation Til as for states of the three- 
dimensional H atom. 
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(b) 2p function 



(c) 2p function 



(a) 2s function 

Figure 1-13. Two-Dimensional States with One Nodal line. 



If the potential energy is independent of angle, the wavef unctions in Fig. 
13b and 13c will have the same energy; however, the wavefunction in Fig. 
13a may be higher or lower than the other two, depending on the exact 
form of the potential. Even worse, we cannot use the nodal theorem to 
determine whether the 3s function in Fig. 14a is above or below the 2p 
functions of Fig. 13b and 13c. 




y 

+ 

+ 



(a) 3s function 



(b) 3p v 



(c) 3p v 



(d) 3d 



xy 



Figure 1-14. Two-Dimensional States with Two Nodal Lines. 



The clue to which comparisons can be made and which cannot is 
apparent from the way that the one-dimensional theorem was proved in 
the previous section. Start with the optimum wavefunction of some nodal 
structure, say Fig. 14b, and change the sign on opposite sides of a single 
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nodal surface to obtain either Fig. 15a or 15b, each of which has exactly 
the same energy as Fig. 14c. 



(a) 
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Figure 1-15. 

The wavefunction in Fig. 15a is an upper bound on the 2p y wavefunction of 
Fig. 13b, and the wavefunction in Fig. 15b is an upper bound on the 2s 
wavefunction of Fig. 13a. 
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Figure 1-16. 

Similarly, starting with the 3s wavefunction of Fig. 14a, we see that 
the wavefunctions in Fig. 16 have the same energy and are upper bounds 
to the 2s wavefunctions of Fig. 13a. However, there is no wavefunction to 
compare the energy of the 3s wavefunctions with those of the 2p 
wavefunctions. Continuing in this way, we can derive the following rela- 
tions: 
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ls < 2s < 3s < 4s 

Is < 2p < 3p < 4p 

2s < 3p ... 

3s < 4p ... 

2p < 3d < 4d ... 

3p < 4d .... 



etc. 
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§1.4 VIBRATION AND ROTATION 

Throughout this course we will focus upon the electronic wavefunc- 
tions for molecules. Thus for an N electronic wavefunction, we determine 



**(1,2 N) 



with energy 



E = 



<* el \H\* el > 



where H is the Hamiltonian for the system. The electronic wavefunction 
and its energy will depend upon the geometry of the molecule. For each 
geometry, we solve for the optimum wavefunction and energy at that 
geometry. For a diatomic molecule the result is a total energy that is a 
function of R (internuclear distance), as indicated in Fig. 17. 



e = e -> 




E = -D. 



Figure 1-17. 

As the nuclei move together or move apart, we imagine the electrons 
readjusting at each instant to reoptimize for that particular R. For a clas- 
sical system, if we started at some particular R, say point b, the nuclei 
would move apart until they reached point c and would then come 
together till point b and would continue oscillating between these points 
(assuming no friction). Starting at point d, the R would continue increas- 
ing until R = °°. On the other hand, if we started at point e, the system 
would stay still. Thus point e is called the equilibrium bond distance 
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(R e ). Starting with the molecule at equilibrium, R e , the energy to pull it 
apart (to break the bond) is called the bond energy, D e . 

For energies below the limit at infinity, we can think of the system in 
terms of two masses (each corresponding to a proton) connected by a 
spring of length R e . However, in quantum mechanics, this spring can 
never be completely at rest. The nuclear motions are described in terms 
of wavef unctions, just as are the electrons, and the kinetic energy of the 
nuclear motions depends on how localized the wavef unctions are. To 
localize the nuclei at exactly R = R« would imply an infinite kinetic 
energy. The result is that for the ground state the nuclear wavef unction 
has the form 



,nuc 



(R) 



JL 



>R 



R. 



Figure 1-18. 

That is, the most likely R is R^ , but the nuclei have a finite probability of 
being found at other R near R e . The result is that the energy of the 
molecule is higher than the absolute minimum (E = -D„) in the energy 
curve by an amount referred to as the zero-point energy. This lowest 
state, Fig. 18, is referred to as the ground vibrational state ($ vib , with v = 
0), and one thinks of the molecule as vibrating back and forth with a fre- 
quency Vq. 

At the bottom of a potential curve, the slope of the energy curve is 
zero and the curvature is positive, so we can write 

E(R) = E(R e ) + ^k(B-R e ) z , (1) 
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where k, the curvature* at the bottom of the well, is called the force con- 
stant. In this approximation (called the harmonic oscillator approxima- 
tion), the vibrational frequency is given by 

"° = w VF- (2) 

where fj, is the reduced mass 

and Mi and M 2 are the masses of the two nuclei. In this case, the zero- 
point energy is given by 

fchi/o (4) 

(where h is Planck's constant). Thus the energy of the ground vibrational 
state is 

E = -A> . 

where 

Dq = D e -%hv . (5) 

The quantity D is the actual energy to break the bond starting with the 
molecule is in the ground vibrational state, and it is the quantity that 
would be measured experimentally. 

In QM, the excited vibrational wavefunctions must be orthogonal to 
the ground wavef unction, leading to the form in Fig. 19, 



, d z E 

k = =- 

dR 2 i? e 
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A vib 



^ 



R e 
Figure 1-19. 

for the first excited vibrational state (v = 1). The excitation energy is 

E 1 -E Q = h i/ • (6) 

(in the harmonic approximation). The separations between vibrational 
states, as in (6), can be determined experimentally, thereby providing 
experimental values for the zero-point energy (4) and for the force con- 
stant k. 

So far we have considered the molecule to lie along the z axis. In 
fact, the axis of the molecule can be oriented along any direction in 
space. Generally, the molecule will be rotating, but the ground rotational 
state is the one for which all orientations are equally likely. This is analo- 
gous to the L = or s state for electrons and is denoted as the J = rota- 
tional state. Excited rotational states have energies of 

Erot = %~ J (J + 1) . 

where / = jaR? is called the moment of inertia. This is analogous to the 
classical rotational energy 

jpclass — _± rZ 

■t-'rot ~ OJ ' 

where L is the rotational angular momentum. Experimentally, the bond 
distance of molecules is often obtained by measuring the rotational ener- 
gies and thereby deriving I and hence /?„. 
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For H 2 and H 2 + , the vibrational energies are 

hv = 4401cm" 1 = 0.546 eV = 12.5B kcal for H 2 

hv = 2322 cm" 1 = 0.288 eV = 6.64 kcal for H 2 + 
and the rotational energies are 

Ej = BJ{J + 1) , 
where 

B - 60.85 cm" 1 for H 2 

B = 30.21 cm -1 for H 2 + . 
In these systems, the total bond energies are 

A> = 36117 cm -1 = 4.478 eV = 103.3 kcal for H 2 

D = 21382 cm" 1 = 2.651 eV = 61.1 kcal for H 2 + . 
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Appendix 1-A. Hermitian Operators 

In §1.1, we found that the basic postulate of quantum mechanics 
implies that the wavefunction i> is normalized, 

<i?\i/> = i , (i) 

and that the time derivative of the wavefunction is determined by the 
relation 

i&ff=fff. (2) 

Here we will show that these conditions imply that H is a Hermitian 
operator 

for all (allowed) functions ^ and %, This roperty results from the 
requirement that the total integrated probability (1) not change with 
time for any superposition of (allowed) functions. 

1-A.l Notation 

First we must establish some notation. For any operator B and any 
functions ^ and fk, we define the jk matrix element of B as 

B jk = <^|B|^ k > = /dr^*(B^ k ) . (3) 

The Hermitian conjugate of B is defined as the operator B ' 

(Bt) jk = <^|B r |^k> = /dT(B^)*tic = <(B^)|(^> (4) 

for all ^j and ^ k (of the Hilbert space). From (3) we see that 

<(B^ 3 )|^ k > = <*c|Stfj>* = B% , (5) 

and hence (4) can be written 

(Bt) jk = B% . (6) 
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If B is equal to its Hermtian conjugate, 

<f j |B t j^> = <^|B|^> , (?) 

we say that B is Hennitian and write 

B r = B . (8) 

1-A.2 Hermitivity of H 

From (1) the total probability of finding the particle somewhere, 

<i\i> = 1 , 

is independent of time. Thus, taking the derivative with respect to time, 
we have 







= < t3H*> + <y £-> = s*\*th + ** If • 



< at r > ' < r at > J ~' at r ' T at 
Substituting the Schrbdinger equation (2) here leads to 

= (ifi)- 1 f dr[-(Hi/)* ^ + if*(Hf)] 



or 



= \-<1/\H\it>* + <i/\H\i»\ 



(9a) 



which implies that the quantity 

E = <i/\H\i/> 
(referred to as the energy) is real. 
Consider now the superposition 

i - Cj^j + c k ^ 

(where C] and C k are numbers, possibly complex) of two states ^ and ^ k 
that are orthogonal 

<^jlfk> = 
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at some time t . Then, since <f\i/>, <^j|^j>, and <^kl^k> are all unity and 
independent of time, it must be that 

C*kCj<Hfc|^> +.C* j C k <f j |^ k > (10a) 

= <^|^>-C* j C J <^ j |^ j >-C* k C k <Vicl*k> 
is also independent of time. Similarly, considering 

V = i Cj^ + C k ^k 

(where i = v^-l), we find that for <f\f> to be independent of time 
requires that 

i C\C j <^ k |^> - i C* 3 C k <^ \i*> (10b) 

also be independent of time. Combining (10a) and (10b) leads to the con- 
dition that 

is independent of time. This leads to 

< a^ ^ < a^ k > 

< at ^ k > < Yj at > 

and hence 

<i>- } \H\^> = <^|£ r |fk> (9b) 

(which also applies to j = k). This relation (9b) must apply to all possible 
pairs of functions ^ and ^ k > and hence the Hamiltonian operator, H, must 
be a Hermitian operator. From this derivation we see that the Hermitian 
property of H results from the assumption that the total integrated pro- 
bability of any superposition of functions is independent of time (con- 
servation of normalization) . 
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1-A.3 The Momentum Operator 

An example is appropriate here. Consider a one-dimensional system 
with coordinates in the range 

< x < a . 

Is the operator 

^ _ h. d 
p ' F dx" 

Hermitian? To find out we consider 

<^j|p r l^k> = /dx(p^)*f k = -|-/dx^-^ k . 
Integrating by parts, this becomes 

<^|p r l^k = - 7-W(aty k (a) -Vj*(0)V*(0)J + ^-/dx^* ^ 

1 1 q Q.X 

or 

<^|p r !^k> = ~t^*i(a)i^(a)-V *j(0)^(0)] + <^|p|^> . (11) 

Thus, the operator p is Hermitian if and only if the boundary conditions 
are such that 

^(atyk(a) = ^(O)rtc(O)- (12) 

for all allowed functions ^ and ^ k . Thus, it is nonsense to say that 
p = — - — is an Hermitian operator; rather, one must say that the opera- 
tor is Hermitian given such and such boundary conditions. Some accept- 
able boundary conditions for the above case are 

(i) ^(a) = 0, ^(0) = (this is the case of a particle in a box), 
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(ii) ^(a) = -^(0) (this corresponds to periodic boundary conditions 
where the point x = a is physically equivalent to the point x = 0; a 
common example is for angular coordinates <p where a = 2tt is 
identically the same point as x = 0), and 

(iii) for systems of infinite range - » < x < +«=, then for <p to be normal- 
ized, <<p \xp> = 1, it must be that <p -» as x -» ±», leading hence to 

9?j*(a)^k(a) -> °° as a -» ±°° . 

If the boundary conditions are such that p is Hermitian, then the 
kinetic energy operator 

T = ar* 1 

is necessarily Hermitian. This follows by applying the Hermitian proper- 
ties of p sequentially, 

<(PxPxfj)l^k> = <(P X ^j)IPx!^k> = <^jlPxP^k> • 
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Appendix 1-B. The Kinetic Energy 

Summary: We will show in this section that the matrix element 

tab = <?al -&V 2 \tp h > (1) 

can be written as 

tab = )£<Vp a -Vp> b > , (2) 

where the dot product is indicated. Thus the kinetic energy of an orbital 
becomes 

taa = £<|Vp a J 2 > (3a) 



in atomic units or 



taa = -if<l^a| 2 > (3b) 



in other units. 



Consider first the integral 

a 



<P 



a I 



s dx z 
Integrating by parts, this becomes 



> 



n * = f f f dxd y dz< p*z ax 2 



dy* & 9^ b 



Q(x= +«) -Q( X = —) -fffdxdydz -^-g^- . (4) 

where 



Q = f/6ydztp\^ 



3x 

(this is referred to as Green's theorem). However, for cp & and p b to satisfy 
the basic postulates of QM, we must have 



and hence 



Thus 



and we obtain 
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<?alP>a> = 1 
<PblPb> = 1 . 

9? a -* as x -» ±< 

Q(x= +«=) = 



Q( x = _«) = o 



< 

<<?a 
< 



ax 2 



> 



< 3^a 

< 3x 



dx > 



Proceeding similarly for the other terms of (1), we obtain 
tab = ViJ dr [V? * a • Vp b ] = % <V? a ■ V? b > , 
and letting ^ a = <p b leads to (3). 



(5) 



(6) 
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Appendix 1-C. The Nodal Theorem 

In this section we consider the eigenfunctions p n for a general one- 
dimensional system, 



where 



H<Pn - EjrfP, 



H = - ■^r-T^r + V(x) 



2M dx 2 

is real (so that the eigenfunctions are real) and 

M = m/& . 

Letting n denote the number of nodes (internal to the boundaries), 
we will show that 



and 



Eq < En 



E n < En+i . 



(1) 



(2) 



That is, the ground state of the system is nodeless and the state with n 
nodes has a lower energy than the state with n+1 nodes. For sufficiently 
singular potentials, the inequalities in (1) and (2) become equalities. 

1-C.l. The Inequalities 

Consider first the functions p and p lt which are the eigenstates of H, 

H? = (?+V)p = E <p 



(3) 



H 9l = (f+V)?! = E^! 



with zero and one node, respectively, as in Figure 20. Here a and b are 
the boundaries of the system (they may be at ± ») and zeroes at the 
boundaries are not counted. 
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,b — 



Figure 1-20. <p and <p l for the proof of the Nodal Theorem. 

First we will show that 

Ej > E . 
Letting c be the location of the node in <p h we consider the region 

a < x < c (4) 

so that both <p x and <p Q are positive. Then from (3) we have 

E = V+ -f-f p 

(5) 

El = v+^t* 

for all points in region (4). Thus the energy difference is given by 

E 1 - E = +-V 9l - X-f ? = -V[?o(t 91) " Pi(?Po)] • (6) 

The integral over all space of the term in brackets is zero since T is Her- 
mitian,* 

b 

/ dx[? (T jpO - ^(f ? )] = <^olT| 9 3 1 > — <t ^ol^i> = ; 

a 

however, the integrand is generally not zero. 



* For a one-dimensional system, the orbitals can always be taken as real. 
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To estimate the sign of (6), we multiply by p^o, then integrate from a 
to c, and then divide appropriately to obtain 



(E,-E ) = 2-. 



where 



(7) 



B = fdxtp f f>\ -fix (pit 



fo 



(8) 



A = / dx jo^o • 



(9) 



Integrating by parts, the first term of (8) becomes 

2M~ 



jfdx^T^ = -^r/dx 



Po 



&9\ 



[ax 2 J 



= ■ + 



2M 



•/dx 



d<p Q 



dx 



d<p ] 



dx 



> _ ' 



2M 



■<Po 



d<P\ 



dx 



2M I ax dx 2 



9\ 



2M 



d?c 



dx 



<P\ ~ <Po 



d? : 



dx 



Combining (7)-(10) together, we obtain 



(Ei - Eo) = 



Since 



1 

J. 

2A~ 


9i 


d^ ' 
i 3x j 


- 


d(p l 
dx 


<Po 



Pi(c) = 



tp (c) > 



(10) 



dx 



< 



and 



we obtain 
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A = f dx ipiVo > , 



Ei - Eq - - 



2A 



■Po(c) 



dx 



> 



(ID 



that is, 



Ei - E > 



(12) 



Thus the nodeless wavefunction has a lower energy than the wavefunction 
with one node. The same proof shows that <p has a lower energy than any 
wavefunction with more than one node. Hence, the ground state of the 
system has no nodal points (inside the boundaries). 

Similarly, the above proof can be applied to the comparison of E n and 
En +1 , that is, the energies for wavefunctions having n and n+1 nodes, 
respectively. The result is that 

E n+1 > E n , (13) 

and hence the eigenstates of a system have energies increasing in the 
same sequence as the number of nodes. 

1-C.2. Singular Potentials 

To obtain (12) we assumed in (11) that 

Po(c) * (14a) 

and that 



dy. 



dx 



* . 



(14b) 



Usually these conditions (14) are satisfied; however, there can be cases 
where the potential is such that one of the quantities in (14) is zero. In 
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this case we have 

Ex - E = , 
and hence the general condition (12) should be 

Ei > E , (15a) 

and (13) should be 

E n+ i ^ E n . (15b) 

For example, consider the case wherein the potential V(x) is so 
strongly repulsive at some point c that all solutions of finite energy must 
have a node at c. In this case, the functions p> and <pi in Figure 20 will 
have the shapes in Figure 21. 





(a) 

Figure 1-21. The first two solutions for a potential sufficiently singular 
at point c. 

From (11) this leads to 

Ei - E = , 

so that E and Ei are degenerate. In this case, the functions <p and <p l will 
have the same shape in each region. At x = c, one of them changes sign 
but since both must be zero there, they have the same energy. 

If the potential is not singular at c, the function <p can generally lead 
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to a lower energy by being positive near c (thereby obtaining a smaller 
kinetic energy). 

The presence of cases such as in Figure 9a complicated our notation 
for the states (in referring to the number of nodes). To keep things sim- 
ple, n refers to the number of points at which the wavef unction changes 

sign. 

Singularities in the potential would also lead to equalities in the ener- 
gies for some excited states, 

1-C.3. A Singular Example 

As discussed above, the equal sign in (1) would occur when the poten- 
tial V(x) is sufficiently singular that the best non-negative wavefunction ? 
has a node at some point. For example, if 

V(x) = , 1 . , 
the potential energy is 

<*<x) | V(x) | rfx)> = / dx i^Wi' + f dx J*& 

Expanding 

<p(x) = <p(xo) + (x - x ) <p'(x ) + 
we see that the dominant term in the integral 



2 



X-Xo| ^ |x-x 



/dx 



■p(x) 



X -Xq 



IS 



?(xo) 2 / |y dx v = ^(x ) 2 ln0 

|X Xq I 
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and hence the energy diverges if ^>(x ) * 0. 

Thus, for <p to yield an E < « (and hence to describe the ground state) 
it must be that 

tp(x) -» as x -» x . 

In this case, 

Eo = Ei 

[assuming no other singularities in V(x)J. 
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Appendix 1-D. The Ground State of Hydrogen Atom 

Summary 

The ground state of the hydrogen-like atom (with nuclear charge Z) 
has the wavef unction 



^(r.tf.p) = N e ^"o , 



(19) 



where 



N f 



■V 



nag 



(23) 



and 



ft 

a ~ 2~ 

me tf 



(16) 



. a °, 



u Zei 
34 R 



(18) 



is denned as the Bohr. The energy of this wavefunction is 

E = -&Z 2 
where 

R " Z 

is the average size of the atom, where a is denoted as the Bohr radius. 

1-D.l. Solution of the Schrodinger Equation 

We will solve for the ground state of the hydrogen atom, that is, the 
lowest solution of 



** v 2 -^- 



2m 



V'dvfl.p) = Ynp{T^,(p) 



(1) 



Since the potential term is independent of angle, the wavefunction has 
the form (see §1.1.1) 



*,*,?) = f(r)Z(fl,p) 



(2) 
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where the angular function Z(tf,p) is a constant for the ground state,* 

Z(tf,p) = 1/V4tT. (3) 



Thus, our chore is to solve 

- #"V 2 f(r) - — f(r) = Ef(r) . 
2m r 

First we must express ^(r) in terms of spherical coordinates, 

z = r cos tf 



(4) 



Since 



x = r sintf cos ?? 



y = r sin tf sin <p 



r 2 = x 2 + y 2 + z 2 . 



(5) 



dr 



dz 



ix.y 



= A-V^+y^z 2 = - 
3z J r 



(where the subscript x,y indicates that x and y are fixed), we see that 



F 



dz 



■f(r) 



t \ f \ 



*.y 



3r 



dz 



k >x,y k y 



3f 



3r 



= -f'(r) • 
r 



Thus, letting 



f"(r) = 


d 2 f 
dr 2 


f(r) = 


df 
dr 



♦The normalization condition is 

21T 1T 

f dtp f sintf dtf [Z(tf,<?)] 



2 - 



= 1 



(6) 



leading to (3). 
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we obtain 



dz< 



•f(r) = 



dz 



-f'(r) 
r 



1 ,,/ \ , dz d 
— f(r) + z ^— -=— 



— f(r) 
r 



= r-f'(r) + 



±-f"(r) - -L-f(r) 



~ fH '« + ^r-r(r) 



q2 a2 

Combinding with -r-^-f(r) and — g-f(r) leads then to 

V 2 f(r) = f'(r) + |-f' (r) . 



(7) 



Substituting (7) into (4) leads to 



2m 


d 2 f 2 df 
dr 2 r df 


- — f(r) = 

r 


= Ef(r) . 



(8) 

as the Schrbdinger equation. Since the potential in (1) goes to zero as r 
-» «, only the states with E < are bound, and hence we take 



E < 



Consider now a sufficiently large r that 



(9) 



Ze" 



« IE 



(10) 



and 



df 



dr 



« 



d 2 f 



dr' 



In this case, (8) reduces to 



d!L- ,z 



dr 2 



= <f 2 f(r) 



(11) 
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where 

^~*0. (12) 

The solution of (1 1) is 

f(r) = e~< T . (13) 

Thus all bound state solutions of (8) must necessarily go to zero exponen- 
tially.* 

Consider now the substitution of the exponential function (13) into 
the Schrbdinger equation (8), 

From (11) the first term of each side (the long-range terms) cancel at all 
r, leaving 

2_df_ = _2L fW = _ 2mZe 2 1_ ( , 
r dr r l{r) tf r { } ' 

Thus the exponential function (13) is an eigenfunction of the Schrbdinger 
equation if 

f = ^fr 1 ■ (15) 

Since £ has the units of inverse length [e.g., see (12)], it is convenient to 
consider the length quantity 

ao = fi?/ me 2 (16) 

(referred to as the Bohr radius or simply the Bohr) as the fundamental 



* Note that in the limit of very large r, the function 

X s1 e-t* 

would also satisfy (8). 
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atomic length leading to 



t = Z/a . 



From (12) the value of £ is related also to the energy, 



E = - 



2m 



■(* 



Thus, 



E = - 



tf 



2m 



a 



= -^z 2 



a 



Summarizing, we find that the wavef unction 



^oM.p) = f(r) = N e _Zr/a ° , 



is an eigenf unction of the Schrbdinger equation with an energy of 



Eo = -^Z 2 


[e»] 
[aoj 


■ 



where 



ao = #Vme 2 . 



(17) 



(18) 



(19) 



(18) 



(16) 



Since the wavefunction (19) is nodeless, we know from the nodal theorem 
that this eigenfunction of H is the ground state of the hydrogen atom. 

For the wavefunction (19), the average value of — is 



R < r 



Thus, the average potential energy is 

Ze 2 



i/> = — 

> ao 



(20) 



< 



> Ze 2 

$> = 

> r 



Zfej 



(21) 
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and the total energy can be written as 



E = fcV = -fc^-= -fc^L. (22) 



In order to normalize the wavefunction (19), note that 

<^|^/> = f r 2 dr / sintf dtf / d<p[f(r)f = 4n f r 2 drf(r) £ 

z 

where the angular integral is 

f dQ = / sin tf dtf / dp = 4rr . 
o o 

Since 



we see that 



No = V? = V 1 ^ (23) 

1-D.2. Analysis of the Wavefunction 

A plot of the orbital along the z axis is given in Figure 22. Note that 
the slope in the wavefunction is discontinuous at z = 0. This singular 
behavior is referred to as a cusp and results from the singular behavior 
in the potential energy at this point. 
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Figure 1-22. The ground state wavefunction, if/ 0t for the H atom (plotted 
along the z axis). 

The Schrodinger equation (3) says that 

Ef(r) 



is equal to 



J^M-isL^) 



for every point r. But, as r -* 0, the term 

- — f(r) 

r 

goes to - oo. Thus, since Ef(r) is finite, the Schrodinger equation requires 
that 



2m 



■V 2 f (r) 



goes to +M as r -» 0. The cusp in the wavefunction leads to a V 2 f(r) that 
goes to +oo as r -> and exactly cancels the negative singularity in the 
potential term. 
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Appendix 1-E. Rydberg Atomic Units 

Another set of atomic units used occasionally employs as the unit of 
energy the ionization potential of the hydrogen atom, 

p _ me 4 /^ 

This quantity is called a Rydberg and is related to the hartree by 

1 Rydberg = )<£ hartree. 

With this choice for the unit of energy (me 4 / 26 s = 1), we cannot use the 
convenient sets of units in §1.2-23. If the unit of length is still taken as 
the Bohr (h 2 / me 2 = 1), then Rydberg units lead to 

me 4 = 2& 2 = 2me 2 , 

and hence we must choose 

|e| = vr 

and 

fi 2 = 2m . 
If we take m = 1, then 

h = V2". 

These units are sometimes used by scattering theorists since the kinetic 
energy of a plane wave, (ft 2 / 2m)k 2 , becomes simply k 2 . Also, some work- 
ers reporting band calculations on solids use Rydberg units. The series of 
books by Slater also uses these units. However, the regular atomic units 
(or hartree atomic units), as described in §1.2.3, are more convenient and 
more common, and we will always use them. 
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Appendix 1-F. Units and Conversion Factors 
1-F.l. SI Units 

In an effort to bring some order to the proliferation of units that con- 
tinues to occur in the sciences, and international group adopted (in 1960) 
what is referred to (in English) as the International System of Units, or 
more commonly as SI units. 

In this system there are seven fundamental units: 



Unit 



Abbreviation 



Physical Quantity 



meter m 

kilogram kg 

second s (or sec) 

ampere* A 

kelvin* K 

mole mol 

candela cd 



length 

mass 

time 

electric current 

thermodynamic temperature 

amount of substance 

luminous intensity 



* Note that these units are not capitalized even though they are derived 
from the names of people. 

From these fundamental units can be derived a number of combined 
units that prove quite useful. Thus, from Newton's Law, F = ma, we know 

that force has units of mass x e *¥? — . an £ n j s convenient to define the 

(time) 

unit of force (newton*) as 

1 newton = 1 kg m sec -2 . 

Similarly, a constant force F exerted over a distance I does an amount of 
work W, so that the unit of energy (joule) is 
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1 joule = 1 newton meter = 1 kg m 2 sec 



-2 



Some common derived units in SI are 



Unit 



Abbreviation 



Definition in Terms of 
Fundamental Units 



Physical 
Quantity 



liter 


I 


newton 


N 


joule 


J 


watt 


W 


pascal 


Pa 


coulomb 


C 


volt 


V 


ohm 


Q 


hertz 


Hz 



10" 3 m 3 




volume 


m kg sec -2 




force 


Nm = m 2 kg sec 


-2 


energy 


J sec -1 = m 2 kg 


sec -2 


power 


Nm -2 = m _1 kg sec -2 


pressure 


A sec 




electric 
charge 


WA" 1 = 

m 2 kg sec -3 A -1 




electric 
potential 


VA" 1 = 

m 2 kg sec" 3 A" 2 




electric 
resistance 


sec -1 




frequency 



The acceptable multiples or fractions to be used for the basic SI units are 
designated by the following prefixes: 
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Fraction 


Prefix 


Symbol 


1CT 18 


atto 


a 


10 -15 


femto 


f 


10" 12 


pico 


P 


10" 9 


nano 


n 


10" e 


micro 


M 


10" 3 


milli 


m 


10" 2 


centi 


c 


1CT 1 


deci 


d 


10 


deka 


da 


10 2 


hecto 


h 


10 3 


kilo 


k 


10 6 


mega 


M 


10 9 


giga 


G 


10 12 


tera 


T 



Although the above fundamental units are convenient for a number of 
quantities, they are quite inconvenient for others. Examples include 



Physical Quantity 



Abbreviation 



Definition in Terms 
of SI Units 



Charge on an Electron 

Atmospheric Pressure 
(at sea level) 



atm 



1.602189 x 10" 19 C 



101.325 Pa 
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Other units are not necessarily superior to the SI quantities, but 
their use is so widespread that scientists must be facile with their use. 
Examples include: 



Unit 



Abbreviation 



Definition 



Relation to 
Fundamental Units 



Kilocalorie kcal 

Angstrom A 

Electron Volt eV 



original definition was: 1 kcal = 4.184 kJ 
energy to heat 1 kg of (defined to be exact) 
H 2 by a temperature 
of 1°K (at 15°C) 



1 A = 1CT 10 m = 10" £ 



cm 



energy change upon 
moving a charge of 1 
electron through an 
electric potential 
of 1 volt 



1 eV = 96.483 kJ mol 



-i 



1-F.2. Units for Coulombs Law 

Conversion between units can sometimes get confusing for coulomb 
interactions. Coulombs law states that the force between two charges Q x 
and Q 2 , separated by a distance R is 

QiQz 



Total Force = 



and hence the enegy of interaction is 



i? 2 



E = /F- dr = ^p- 



In cgs units, we define the electrostatic unit of charge (esu) as that 
charge which leads to a force of 1 dyne when the charges are separated 
by 1 cm. Thus 

Force (dynes) = \^[\ z 
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and 

Energy (ergs) = l£fpil]L . (1) 

i?(cra) 

In SI units, the unit of charge (the Coulomb) is defined in terms of a 
current (1 ampere = 1 Coulomb /sec), and the units of current are 
related to force and distance through a different force law (magnetic 
induction). The relationship between Coulombs and esu turns out to be 

1 Coulomb = 



(2.998 10 9 ) ' 

where the 2.998 comes from the speed of light (2.998 10 8 m/sec). Since 

Force (newtons) = 10~ 5 force (dynes) 

R (meters) = 10 -2 R (cm) , 
the Coulomb law becomes 

Force (newtons) = [8.988 iQ^]l£j^^^- . (2) 

and 

Energy (joules) = [8.988 10 9 ] ^(coulomb)] 2 

K\cm) 

In SI units, the constant in this expression is generally written as 

8.988 x 10 9 = — = — , 

47TS 

where 

s = 8.854 10" 12 C z N- l m- z (3) 

is called the permittivity of a vacuum. This leads to 
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for the coulomb energy. 

In this course the Q x and Q 2 are always some multiple of the funda- 
mental charge on a proton (or electron), 

le = 1.602 10" 19 C 

= 4.803 10" 10 esu . 

Thus we will write 

Q\ = q x e and Q z = g 2 e , (5) 

where q l and g 2 have no units. In addition, we will often write R in terms 
of Bohr radii, e.g., 

R = r a , (6) 

where r has no units. In this case the coulomb energy becomes 

?1<?2 



Energy = 



ttr 



(7) 



where no Arnso factor is included. In (7) the unit of energy is 



e 2 
1 h = 1 hartree = — 

a 



= 27.2116 eV - 2625.5 kJ/moi = 627.51 keal/mol 



(8) 



The fast way to calculate atomic level coulomb energies in various units is 
to first express all distances and charges in atomic units, as in (5) and 
(6), calculate the energy using (7), and then to convert from atomic units 
to SI units using (8). 

Example: Calculate the interaction of two protons at a distance of R = 5 

o 

A. 

Answer: , 
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E(h ° ) = R^r= 9k= 01058 ^ 



= 2.88 eV = 278 kJ/ mole = 66.4 kcal . 

Useful conversion factors here are 

e 2 = 14.3998 eV A = 332.059 (kcal/mol) 1 

1-F.3. Units for Mass 

In atomic units, the mass of the electron is unity; however, this is not 
to be confused with the atomic mass unit (amu) which is the standard for 
relating the masses of atoms. The modern convention is to define the 
dominant isotope of C (i.e., 12 C) as having a mass of 12.0000. In these 
units, the mass of the hydrogen atom is 

1.00783 amu, 
the mass of a proton is 

1.00728 amu, 
and the mass of an electron is 



amu. 



1822.89 
Thus, in (hartree) atomic units, the mass of the proton is 

1.00728 x 1822.89 = 1836.16 
The conversion to SI units is 

1 amu = 1.660566 x io~ 27 kg. 

1-F.4. Energy Quantities for Photons 

The wavelength (A) and frequency {v) of light are related by the speed 
of light (c) 
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Xu = c . 

Thus, since c = 2.99792458 x 10 10 cm/sec, the frequency for yellow light (A 
= 600 nm) is 

v = g x ^ = 5 x 10 16 Hz 
6 x 10 -7 

(1 Hz = 1 cycle per sec), and the wavelength for KZLA (v = 94 MHz) is 

X = JX 1U = 319 cm = 3.19 m . 
94 x 10 6 

In the quantum description of light, the energy of a photon is given 
by 

E = hi/ 

(where h = 2-nh. is the original Planck's constant); thus the energy of a 
photon of light can be expressed as 

E = hi/ = ^- = hcT> , 

A 

where v - —is called the inavenumber (and denoted as cm -1 ). Substitu- 

A 

tion for the known values of h and c leads to 

E ( eV > = Wn^ 5 = 8065 - 48 V (Cm_1) 

v(cm x ) = 



X(nm) 

Thus, when an electron decreases its energy by 1 eV (dropping into a 
lower energy state), it may emit this energy as a single photon with 
wavelength X = 1240 nm = 1.24 microns or wavenumberT> = 8065 cm" 1 . 

1-F.5. Other Energy Relations 

Chemists often use the energy quantity kilocalories per mole which 
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we will abbreviate as kcal, 

1 eV = 23.06036(14) kcal. 

Recently, emphasis has been placed on SI units in which the kilojoule per 
mole (denoted as kJpm) is the energy unit, 

1 kcal = 4. 18400 kJpm. 

The atomic unit of energy, the hartree, is kind of large for convenient 
use and we will often use the millihartree, denoted as mh. Relations 
between these units are 

1 eV = 36.7490 mn = 23.0604 kcal = 96.4847 kJpm = 8065.48 cm" 1 
1 mh = 0.27212 eV = 0.627511 kcal = 2.62550 kJpm = 219.475 cm" 1 . 

The average thermal energy of an oscillator at room temperature is 

kT = ^-eV = 1 mh = 0.6 kcal = 200 cm" 1 = 2.4 kJpm. 

The strength of the H 2 + bond is 

2.5 eV = 92 mn = 58 kcal = 241 kJpm. 
The vibrational energy 6w e of H 2 + is 

3000 cm H = 14 mh = 8.6 kcal = 0.37 eV = 36 kJpm. 

1-F.6 Examples 

The fundamental constants are experimentally determined and 
hence the best values for them change with time. For aid in calculating 
these constants in the future, we summarize the procedure. 

1 (eV/atom) x 1.6021892 x 10" 22 (kJ/eV) x ^^ x 6.022045 



x 10 23 (atoms /mol) = 23.06036 kcal mol" 
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Therefore, 1 eV/atom = 23.06036 kcal mol" 1 . For force constants we use 

k = Mu 2 = _mass_ = Energy 
(Time) 2 (Length) 2 

, h 27.21161 » 1.6031B93 * 10~ 12 erg 

i ail — — — — p 

ao 2 (0.52917706 * 10" 8 cm) 2 

= 1.556919 * 10 8 dyne/ cm 

= 15.56919 mdyne/ A . 

1-F.7. Conversion Factors 

Included herein are the fundamental constants as of 1973. The refer- 
ence is E. R. Cohen and B. N. Taylor, J. Phys. Chem. Ref. Data, 2, 717 
(1973), and some tables from this reference are included on the following 
pages. 

Commonly used constants are:* 

1. lbohr = 0.5291771 (4) A 

2. 1 hartree = 27.21161 (7) eV 

1 hartree = 627.5096 (4) kcal/mol = 219474.7 cm" 1 

3. 1 eV = 23.06036 (14) kcal/mol = 96.483 kJ/mol 
1 eV = 8065.48 (2) cm" 1 

leV= 1239.652(3) 
\(nm) 



1 kcal = 4. 18400 kJ = 349.755 cm 



-l 



* A number in parentheses indicates the standard deviation in the error for the 
last quoted digit. 
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5. Atomic Masses: 1 amu = 1822.887 (1) m e 

This is used to convert atomic masses (mass 12 C = 12.000) to har- 
tree atomic units 

6. Force Constant: 1 au = 15.56919 mdyne/A 

7. Gas Constant: R = 1.98719 (7) cal mol" 1 K" 1 
Gas Constant: R = 82.0568 cc atm mol" 1 K" 1 

8. Avogadro's Constant = 6.022045 (31)* 10 23 molecules /mol 

9. Fine Structure Constant: 1/a = hc/e 2 = 137.03604 (11) 

10. Dipole Moment: 1 au = 2.541765 (8) Debye 

Dipole Moment: 1 au = 2.541765 (8) * (lCT ia esu*cm) 

11. Quadrupole Moment: 1 au = 1.345044 * 10" 28 esu cm 2 
Quadrupole Moment: 1 au = 1.345044 Buckingham 

12. Electric Field Gradient: 1 au = 3.241391 * 10 15 esu cm" 3 

n n™,^™*. TPr^^c e2 f a \r\ - 14.3998 e 2 / Vrtol >> _ 332.059 

13. Coulomb Energies: ^-(eV) = R ^ , , -g-(kcal) = R ^\ 
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Table 1-F.I. Current Values for Selected Fundamental Constants. 



Quantity 




Symbol 


Value 


Uncer- 


Units 




tainty 
(ppm) 


SI 


cgs ° 


Elementary charge 




e 


1. 6021892(46) 
4. 803242(14) 


2.9 
2.9 


io-^c 


10" 20 emu 
10" 10 esu 


Planck constant 




h 

fi = h/2* 


6. 626176(36) 
1.0545887(57) 


5.4 
5.4 


10" 34 J-s 
10" 34 J-s 


10"" erg- s 
10" 2 erg-s 


Speed of light in vacuum 


c 


299792458(1. 2) 


0.004 


10 8 m»s" L 


it\W -i 
10 cm-s 


Fine- structure constant, 
[,i c 2 /4*](eytic) 


S- 


7. 2973506(60) 
137. 03604(11) 


0.82 
0.82 


10" 3 


to" 3 


Electron rest mass 




m e 


9. 109534(47) 
5. 4858026(21) 


5.1 
0.38 


1<C kg 
10^ u 


io:r g 

10 4 u 


Proton rest mass 




m p 


1.6726485(86) 
1.007276470(11) 


5.1 
0.011 


10- 2T kg 
u 


lO" 24 g 
u 


Ratio of proton to 
electron mass 




V m e 


1836. 15152(70) 


0.38 






Atomic mass unit, 
10" 3 kg. mol" 1 N A l 




u 


1. 6605655(86) 


5.1 


10" 27 kg 


io- M g 


Ratio amu to in 
e 




u/m e 


= 1822. 887(l)m e 








Bohr radius, 
[^ c 2 /4»]- l (BVm e e 2 


) = a/AwE^ 


a o 


5. 2917706(44) 


0.82 


10-" m 


10" B cm 


Rydberg constant, 
[n cV4*] 2 (m e eV4«ll 


3 c) 


R oo 


1.097373177(83) 


0.075 


lO'm" 1 


10 5 cm" 1 


Avogadro constant 




N A 


6.022045(31) 


5.1 


10 33 mol* 1 


10 23 mol" 1 


Faraday constant, N*e 


F 


9. 648456(27) 
2. 8925342(82) 


2.8 
2.8 


tO^'mol" 1 


10 3 emu* mol" 
10 14 esu -mol" 


Bohr magneton, [c] (efi/2m c) 


^B 


9. 274078(36) 


3.9 


lO-^J-T" 1 


lO^erg-G" 1 


Nuclear magneton, 
[c] (efi/2m p c) 




**N 


5. 050824(20) 


3.9 


lO'^J-T" 1 


10" 24 erg-G" L 


Molar volume of 
ideal gas at s. t. p. 




m 


22.41383(70) 


31 


10" 3 m 3 -mol" 1 


10 3 cm 3 -mor 1 


Molar gas constant, 


PoV m /T 


R 


8.31441(26) 


31 


J'mol'^K" 1 


10 T erg -mol" 1 -K" 1 


(T = 273.15K;p = 
= 1 atm) 


101325 Pa 




8.20568(26) 
= 1. 98719(7) kcal/mol 


31 10" 
K 


' m 3 *atm.'mol" -K' 


10 cm -atm-mor^K" 1 


Boltzmann constant, 


R/N A 


k 


1.380662(44) 


32 


10"* 3 J-K" 1 


lO'^erg-K" 1 



Note that the numbers in parentheses are the one standard- deviation uncertainties in the last digits of the quoted value 
computed on the basis of internal consistency, that the unified atomic mass scale U C = 12 has been used throughout, that 
u = atomic mass unit, C = coulomb, F = farad, G = gauss, H = henry, Hz = hertz = cycles/s, J = joule, K = kelvin 
(degree Kelvin), Pa = pascal = N- m" a , T = tesla (10* G), V = volt, Wb = weber = T • m 2 , and W = watt. In cases where 
formulas for constants are given (e.g. , R ), the relations are written as the product of two factors. The second factor, 
in parentheses, is the expression to be 



charge in electrostatic units. 



used when all quantities are expressed in cgs units, with the electron 



H 

In order to avoid separate columns for "electromagnetic " and *telectrostatic" units, both are given under the single 
heading cgs Units. " When using these units, the elementary charge e in the second column should be understood to 
be replaced by e m or e , respectively. 
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Chapter 2. The Chemical Bond: Hf and H2 
Summary: 

In this chapter we consider the two states of H 2 + 

<Pg = Xl+Xr 

9u - Xi ~Xr 

(the LCAO wavef unctions) arising from bringing a proton up to the ground 
state of hydrogen, and we consider the two states of H 2 

$g = XlXr + XrXl 
$u = XlXr ~XrXl 

(the VB wavefunctions) arising from bringing together two hydrogen 
atoms each in the ground state. As expected from the nodal theorem, 
the g state (symmetric) is the ground state for both systems. Indeed, in 
each case we find that the g state leads to bonding, while the u state 
leads to a repulsive potential curve. The g state of H 2 + leads to an 
increase of the electron density in the bond region; however (contrary to 
popular belief), this leads to an increase in the electrostatic interactions, 
thus opposing bond formation. A bond is formed because of a very large 
decrease in the kinetic energy due to the molecular orbital having a 
significantly decreased gradient in the bond region. The bonding of the g 
state of H 2 arises from the same term (modified by an additional overlap 
factor due to the second electron). 

The potential curves for both states of both molecules are dominated 
by exchange terms of the form 



.X — 



' fi 1 + S 



E X = - 



1 -s 

for H 2 + and 



E* = 
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2St 



E* = - 



1 + S 2 

2Sr 

1-S 2 



for H 2 , where S is the overlap of the atomic orbitals. The quantity r is the 
quantitative manifestation of the decreased kinetic energy (and 
increased potential energy) arising from interference of the xi and \t 
orbitals. It has the form 

for large R. Thus, at large R the bonding of H^ is proportional to S, while 
the bonding of H 2 is proportional to S 2 . Consequently, for large R the 
bond energy of H 2 exceeds that of H 2 . For small R , where S * 1, the bond 
energy of H 2 is approximately twice that of H^. The u states are far more 
repulsive than the g states are attractive (due to the 1 ± S and 1 ± S 2 
terms in the denominators of £* and E x ) . 

We also examine the molecular orbital (MO) wavefunction for H 2 

*g m (1.2) = <P g (D?g(2) , 

which provides a simple description of the ground and excited states for 
small R. For large distances, the ionic terms implicit in the MO wavefunc- 
tion lead to an improper description. 
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2.0 The Chemical Bond in Hg + and Hg 

Many atoms will combine with other atoms to form a strongly bound 
molecule. The point of this chapter will be to establish the origin of the 
chemical bond for the simplest one- and two-electron systems. 

We will observe the following conventions on notation in this and fol- 
lowing chapters: Lower-case Letters will be used for one- particle 
wavef unctions (93) and energies (s), while upper- case letters will be used 
for many- particle wavef unctions ($) and energies (E). 

2.1 Origin of the Bond in H 2 + 

We first consider the smallest possible molecule, H2, consisting of one 
electron plus two protons separated by a distance R. This system is 

sketched in Figure 1, where the two protons are denoted as a and b. 

e- 
A 
' I 

'a.' ; 

-' i 

a * ff----*b 

Figure 2-1. Coordinates for H£. 

2.1.1 LCAO Description 

Consider first the case with R = *>. With the two protons infinitely far 
apart, the ground state is obtained by placing the electron in the Is orbi- 
tal of one of the other of the two protons. This leads to the two states, 

H H + (la) 

and 

H* H (lb) 

which are described by the wavef unctions,* 
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■ <p = xi = Ne" ra (2a) 

and 

9? = Xr = Ne" rb , (2b) 

respectively, where xi and x r denote hydrogen Is orbitals centered on the 
left and right protons. 

For finite R, the exact wavef unctions no longer have the atomic form, 
but useful approximate wavefunctions can be obtained by allowing the 
wavefunction to be a (linear) combination of the atomic orbitals in (2), 

9 = CiXt + C rXr ■ (3) 

This simple type of wavefunction is often referred to as LCAO for a linear 
combination of atomic orbitals. We will find that the optimum LCAO 
wavefunction is the symmetric combination, 




9g = tta+Xr)/D 9 y x S V (4) 



(where D g is a normalization factor). The other combination of the orbi- 
tals (5.3) is the antisymmetric combination, 




<Pu = (Xl ~Xr)/A. X S~ (5) 



(A* is the normalization factor). 

The energies for the wavefunctions cp g and <p u in (4) and (5) are shown 
as a function of R in Figure 2. Here we see that the g state is strongly 



* N is the normalization factor. 
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zd 4.0 ao 8.0 

INTERNUCLEAR COORDINATE (BOHR) 



10.0 



Figure 2-2. The energies of the LCAO wavefunctions for Hg*. 




-4.0 



-2.0 2.0 

2 COORDINATE (BOHR) 



Figure 2-3. The densities p g and p u for the LCAO wavefunctions of Hi com- 
pared with the superposition of atomic densities, p*. 
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bonding (that is, the energy drops as the nuclei are brought together), 
while the u state is strongly antibonding (the energy increases as the 
nuclei are brought together). The objective of this section will be to 
understand the origin of the bonding and antibonding charactered exhi- 
bited by the tp g and <p u states. 

2.1.1a Ele ctrostatic Energy 

First we consider the electron density, 

P9 = ■*# = -hriXi +Xr) 2 = ixf + X? + 2XiXr)/D£ . (6) 

Integrating p| over all space must give one electron 

<*!> = i (?) 

and similarly 

<Xi> = 1 
<X?> = 1 

(recall that these are just the Is orbitals of H atom). Thus (6) leads to 

1 + 1 + 25 
i = 

where 



1= AT 



S = <%Xr> (8) 

is called the overlap of the two atomic orbitals. Consequently, the nor- 
malization condition in (4) is 

D g = V2(T+-S") • (9) 

If there were no interference terms in (6), the density would be* 

P cl = %{X? + X?) (10) 



* The factor of }£ leads to the required condition <p ci > = 1. 
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but because of the interference terms the density near the bond mid- 
point is increased, as shown in Figure 3. This result has given rise to the 
prevalent idea that the chemical bond arises from the increase in the 
electron charge density in the bond region. The idea is that an electron 
in between the nuclei attracts both nuclei, holding them together to form 
the chemical bond, 



p + e p* 



(u) 



77iis reasoning is false as will now be demonstrated. The total potential 
energy is given by 



g 2 g 2 g 2 



7(r) = - — + 



r b R ' 



(12) 



as sketched in Figure 4. 



V(r) 



V(r) = Z- 
R 



V 




Figure 2-4. The nuclear attraction potential V(r) for H2". 

as sketched in Figure 4. Here we see that the best place for the electron 
(i.e., lowest energy) is at a nucleus (r tt = or r b = 0), not at the bond mid- 
point. From Figure 3 we observe that the increase in charge at the bond 
midpoint is at the expense of charge near the nucleus. Thus, in forming a 
bond, the charge is transferred from a low energy region (near the 
nucleus) to a high energy region (the bond midpoint), an effect that 
should operate against bond formation. Indeed, this is the case, as 
shown in Figure 5, where 
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UJ h- 

h- cr 

Si 

h-o: 

< UJ 

-J z 

cr 



0.20 



0.10- 



-0.10 - 




2.0 4.0 6.0 8.0 

INTERNUCLEAR DISTANCE (BOHR) 



10.0 



Figure 2-5. The relative potential energies, V g and V,, for the LCAO 



wavefunctions of H£. The absolute values are obtained by noting that 
, B = -1.0atR = «. 



V e = V, 




4.0 



2 COORDINATE (BOHR) 



Figure 2-6. Comparison of the <p g and <p u LCAO's of Hg + with the hydrogen 
atom orbital, <p la . All wavefunctions have been normalized. 
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V g = <p|7(r)> 
is the total potential energy for the <p g wavefunction. 

Our conclusion then is that transfer of electron charge into the bond 
region leads to repulsive electrostatic interactions. The fact that the 
bonding state leads to such a transfer indicates that the origin of the 
bond lies in the other contribution to the energy, the kinetic energy, as 
will be discussed next. 
2.1.1b Kinetic Energy 

A qualitative prediction of changes in kinetic energy upon bond for- 
mation is easy. The kinetic energy is the (average) square of the gradient 
of the wavefunction, T = (n 2 /2m) <|V^| 2 >. Superimposing two atomic orbi- 
tals symmetrically as in <p g leads to a large decrease in the slope in the 
bond region (see Fig. 6), and hence a large decrease in the kinetic energy 
(see T g of Fig. 7), 

T g = (n 2 /2m)<(Vp ff ) 2 > 

resulting in a strong bond. On the other hand, the antisymmetric combi- 
nation in <p u leads to a large increase in the slope in the bond region (see 
Fig. 6) and hence the kinetic energy opposes bond formation (see T u of 
Fig. 7). 

The resulting total energies are given in Fig. 2, where we see that <p g 
is strongly bonding while <p u is strongly antibonding. 
2. 1 .2 Bonding to p Orbitals 

Above we found that it is the change in the kinetic energy that dom- 
inates the energy changes in the LCAO description. Basically, if two 
atomic orbitals are superimposed so that no new nodal planes are 
created, as in Fig. 8a, then the kinetic energy drops significantly due to 
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Figure 2-7. The changes in the total kinetic and potential energies for the 
g and u LCAO wavefunctions of Hi. The actual values at R = « are 
Ta(-) = T g («) = + % an d v,(-) = V^-) = -i. 
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Figure 2-8. (a,b) Symmetric and antisymmetric superposition of Is 
atomic orbitals. (c,d) Symmetric and antisymmetric superposition of 2^ 
orbitals (oriented along the axis). 
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the decrease in the gradient of the orbital in the internuclear region. 
This is a general phenomenon and depends only on the fact that in the 
bond region the gradients of the atomic orbitals are in opposite direc- 
tions (contragradient) so that (symmetric) superposition of the orbitals 
leads to a decrease in the gradients. 

Conversely, superposition of orbitals so as to lead to a new nodal 
plane, as in Fig. 8b, gives rise to much steeper gradients and a much 
larger kinetic energy, leading to antibonding potential curves. 

These results are not limited to superimposing Is orbitals. Consider, 
for example, a bond between p z orbitals on two atoms (assuming z to be 
the internuclear axis), as in Figure 9. The plus combination leads to a new 
nodal plane (higher gradients) and antibonding (see Fig. 8d), while the 
minus combination leads to lower grandients and bonding (see Fig. 8c). 
Similarly, bonding of the p x orbitals leads to Figure 10. Now the minus 
combination leads to a new nodal plane and antibonding, while the plus 
combination leads to bonding. 

2. 1.3 The Optimum Distance for Bonds 

There is a natural optimum range for the effects that dominate bond- 
ing: 

i) If R is very large (near «), there is a large region in which the gra- 
dient is decreased; however, at each point, one or the other of the 
two orbitals has a very small gradient, so that the decrease in the 
gradient is very small (and goes to zero as R -» «=). The result is a 
small bonding contribution for large R. 

ii) If R is very small (near 0), there is a large decrease in the gradient; 
however, the region of this large decrease is only the small region 
between the nuclei (which goes to zero as R -> 0). [The latter effect is 
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Figure 2-9. Bonding between p* orbitals. (a) and (b) are schematic 
diagrams of the shape of the orbitals in the xz plane, (c) and (c) are plots 
of the orbitals along the z axis. 
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Figure 2-10. Bonding between p, orbitals. 
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illustrated in Figure 11, where the left side is for R near optimum and 
the right side is for small R.] 

Thus the optimum bond is formed at an intermediate distance where the 
gradients are large and opposite (contragradient) for a large region. For 
the hydrogen Is orbital, the optimum distance is about 2a , which is just 
the sum of the atomic radii. For a p orbital, the optimum decrease in the 
gradient occurs when the outer lobes are separated significantly, as illus- 
trated in Fig. lie. 
2.1.4 Symmetry Considerations 

The H2" molecule has a great deal of symmetry. In quantum mechan- 
ics., symmetry in the molecule generally leads to symmetry in the 
wavefunction, and knowledge of these symmetries can aid us both in solv- 
ing for the wavef unctions and in reasoning qualitatively about them. For 
the time being we will concern ourselves with only one of the symmetries 
in H2", namely, the inversion symmetry. 

2.1.4a The HamiLtonian 

First we need to consider the form of the Hamiltonian for H 2 *. Using 
the coordinate system of Figure 1, the full Hamiltonian for H 2 + is 

**> ■ - fr '• - fr* -£*-¥--¥-* ^ ■ "* 

We will simplify (13) by assuming the nuclear masses to be infinitely 
heavy (M a = M b = °°), by taking the nuclear charges as unity (as appropri- 
ate for H2") and by using atomic units (h = m = e = 1). This reduces (13) 
to 

*W) = -W-J--i- + F- (14) 

We will group together all the terms depending upon the coordinates of 
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Figure 2-11. Illustration of the effect of R on the contragradience of orbi- 
tals. In each case the R for the left case is near optimum, while the R for 
the right case is too small. 
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one (and only one) electron as 

Ml) = -J*Vf + v(l) (15) 

(referred to as the one-electron Hamiltonian), where 

v(l) = - -^- i- (16) 

is the nuclear attraction term (arising from the attractive electron- 
nuclear interactions). This leads to 



£(H 2 + ) = h(l) + £ 



(1^) 



The exact (electronic) wavefunction of Hg" is obtained by solving 

H<p{l) = E*(l) , (18) 

where H is given by (17). Substituting (17) into (18) and rearranging, we 
obtain 

hsp = b<p , (19) 

where 

e = E - 1- (20) 

is referred to as the electronic energy. Although (19) may appear to 
involve only the electronic coordinates r, the internuclear coordinate R is 
involved implicitly since it determines the spacing of the attractive terms 
in v, (16). In solving for the wavefunction of Hg", we choose an R and solve 
(19) to obtain the electronic wavefunction <p{r) and the electronic energy 
s. We then choose a new R and again solve (17), obtaining a new <p(r) and a 
new s. The result is an electronic wavefunction <p (r) and an electronic 
energy s, each of which is parametric ally dependent upon R. (This pro- 
cedure is referred to as the Born-Oppenheimer approximation.) 



-14- 

2. 1 . 4b Inversion Symmetry 

The operation of inversion through the origin of a coordinate system 
leads to the changes 

x -» — x 

y -» -y (21) 

z -» — z 

in the coordinates, and will be denoted as L 

Taking the origin of the coordinate systems as the bond midpoint in 
Fig. 1, the inversion of the coordinates of the electron leads to Figure 12. 
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Figure 2-12. The effect of inverting the coordinates of the electron. 

The electron is now r tt from the right nucleus (b) and r b from the left 
nucleus (a). However, since the nuclear charges are the same, the poten- 
tial terms in the Hamiltonian are the same. 

Upon inversion, the kinetic energy terms in H are also unchanged 

v* = JL + .sL + .2L -> * + _JL_ + _2i_ = v2 
ax 2 ey 2 az 2 a(-x) 2 a(-y) 2 d(-z) 2 

and hence the Hamiltonian is invariant upon inversion of the electronic 
coordinate (through the bond midpoint). 

Now consider that we have solved (18) to obtain eigenstates of H 2 \ 

H<p = E<p . (22) 
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and apply T to both sides of (22). The result is 

!(£$) = E(?9p) , (23) 

which we could write as 

H(-r)<p(-r) = E?(-r) , (24) 

indicating the result of inversion. But H is invariant under? 

H{-r) = H{t) (25) 

so that (24) becomes 

H(t)?{-t) = E?(-r) 
or 

H(r)(ltp) = E(T?>) . (26) 

Equations (22) and (26) state that <p and I p are each eigenfunctions of 
exactly the same Hamiltonian with exactly the same energy. There are 
two possibilities here: 

(i) The state is nondegenerate, in which case I<p and <p must be propor- 
tional to each other, or 

(ii) the state is degenerate, in which case l<p and <p may be linearly 
independent* functions. 

First we consider that the state is nondegenerate. In this case 

T9? = \p , (27) 

where \ is some constant. But, applying T twice leads to 

x -► x 

y -* y 

2 -> Z , 



* i.e., not proportional. 
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and thus must return the original function 

!Mr) = T<p(-t) = <p(r) , (28) 

whereas applying T to (27) leads to 

I 2 ? = XI? 
and using (27) on the right side leads to 

I 2 ? = \ z y> . (29) 

Combining (28) and (29) leads to 

*(r) = X 2 ?(r) 

or 

X 2 = 1 

leading to 

X = ±1 . 

That is, nondegenerate states of H 2 + must be either symmetric under 
inversion (X = +1) or antisymmetric (X = -1). Wavef unctions with these 
symmetries are denoted with g (for gerade or even in German) or u (for 
ungerade or uneven), as in ? g or p u . 

Consider now the case (ii) of a degenerate state with % not propor- 
tional to (p. We can form two new functions, 

<p u = <p- \<p 

such that each function is still an eigenfunction of H (with the same 
energy) H<p g = E(p g , H<p^ = E?> u but such that one function is gerade 

l?>g = <p g , 

while the other is ungerade 
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Thus, in this case also the eigenf unctions of H are g or u J 

The same procedure can be used for higher degeneracies, and hence 
the conclusion is that for any H invariant under inversion, each eigen- 
state can be taken as either g or u. Examples are given in Figures 8, 9, 
and 10. 

2.1.4c The Nodal The orem 

An ungerade wavefunction for H 2 + necessarily must change sign at the 
plane passing through the bond midpoint. Consequently, from the nodal 
theorem we know that the ground state of Hi vjill be ag state. * 
2.1.5 The Exchange Energy 

There is a direct relationship between the bonding observed in p g and 
the antibonding observed in <p u , both being dominated by changes in 
kinetic energy as the bond is formed. We will now obtain an explicit form 
for this relationship. 
2.1.5a The Classical Energy 

Consider first the wavefunction for H 2 + with no super imposition of 
atomic orbitals, 

9* = Xi ■ (30) 

We refer to this as the classical wavefunction because it does not have 



' If a certain state is doubly-degenerate with wavefunctions $0 a and y^, then start- 
ing with just one function, say <p&, we generate both a g function and a u function, 
^ga = <P& + l^a a nd tp u& = tp & — I^ a . If these functions are both nonzero, then (p^ 
will be a linear combination of 9? ga and $ff ua , and nothing need be done with it. 
However, if <p & were already g or u, then (p\, is needed to generate the second func- 
tion. 

* Since there is no singularity at the nodal point, the inequality in the nodal 
theory applies, resulting in E g < E u . However, for R = °°, even the u wavefunction 
is zero at the midpoint, and hence the lowest g and u states are degenerate. 
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interference effects arising from superposition of atomic orbitals. The 
energy of this wavefunction, 



is nonbonding, as shown in Figure 13. Using (17) we obtain 



£ ci = <xi\H\ X i> = <Xi 

< 






> 

Xi> 
> 



(31) 



£atom + < Xl 
< 



_1_ 



> 

Xi> 

> 



where 



E atom ~ < Xi 

< 



_^V2-A- 



Xi> = -^(e 2 /a ) 
> 



Thus 



where 



£ ci = 8«tom+*^ eI • 



6 v ' = 



<xf 



rb 



> + 



R 



is repulsive. 

2.1.5b The Exchange Energy 

Now we consider the wavefunction <p g with energy 
<l + r I H I i + r> 



£a = 



<J + rU + r> 



<l\H\l + r> 

<l\l + r> 



Since 



<J|i7|i+r> = <l\H\l> +<i\H\r> 



= £ ci + <l\H\r> 



(32) 



(33) 



(34) 



(35) 



<l\l +t> = <l\l> + <i I r> = 1 + S , 
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we obtain 



= £ ct + <l\H\r> _ 4 + <l\H\r> -e c1 3 



1 + S 



1 + S 



g 1 + s 



where 



(36) 



r = <i!tflr> -Ss cl 



hir + S 



R 



-S 



^ + F 



or 



T = hj r - S hu ■ 



Similarly, 



e« = 



- <* -rlff|Z -r> _ <l \H\l -r> 



<Z -rU -r> 



<l I -r> 



(3?) 



- £ ci -<*i#[r> 

1 — 55 



or 



= *<* - 



1 -S ' 



(38) 



where r is again given by (37). Thus the interference resulting from 
superimposing the xi and Xr wavefunctions can be viewed as corrections 
upon the classical energy, 



s„ = s ct + s- 



S-u ~~ ** "*" £ti i 



(39) 



where the correction terms 






1 + 5 



(40) 
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are referred to as interference, exchange, or resonance terms. These 
energies are shown in Figure 13, where we see that e| favors bond forma- 
tion, while e£ opposes bond formation. 

The classical energy, as defined above, is the total energy of the sys- 
tem if the wavefunction is forced to remain an atomic orbital as R is 
decreased. The exchange part of the energy is the change in the energy 
due to the interference of xt a nd Xr. that is, due to exchange of the elec- 
tron between the left and right centers. As shown in Figure 12, e cl is 
weakly antibonding, and hence bonding in the g state of H^ results from 
the exchange energy ej. On the other hand, the exchange term e$ for the 
u state is strongly repulsive, resulting in a strongly antibonding potential 
curve. 

These quantities ef and e£ constitute a quantitative representation of 
the effects discussed qualitatively in §2.1.1. Thus, the decrease in kinetic 
energy for the g states resulting from the decrease in the average gra- 
dient in the wavefunction yields a large negative contribution to r. The 
increase in the potential energy for the g state arising from the shift of 
charge from the nuclear to the bond region yields a positive contribution 
to t. The net result is a negative r, leading to a negative value for 

£ | = t/(1 +s) 

and a positive value for 

*5 = -t/ (l -S) . 

2.1.5c Comparison of g and u States 

For Large R, where the overlap S is nearly zero, we see that (40) leads 
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to 



e# = t 

so that the bonding in the g state and the antibonding in the u state are 
equal. 

However, for small R the (1 + S) and (1 - S) terms lead to asym- 
metry, where the antibonding state is several times more antibonding 
than the bonding state is bonding. Thus, at R = 2.5ao = 1.32 A , we have 

S = 0.45B3 
t = -0.1063 (e 2 /ac) , 
leading to 

£* = -0.0743 (e 2 /a ) 
£* = +0.20939 (e 2 /a ) , 



whereas 



s cl = 0.00943 (e E /a ) 



2.1.5d Analytic Results 

Explicit evaluation of the various quantities involved in the energy of 
Hg" is carried out in Appendix A, leading to 



S = [1 + R + ±-R z ] e-* 



(41) 



*« = -JS+d + ^e-* 8 



T = 



3 R R 



e -R „ (! + ±) (i + R + A-jjS) e -3i? 

7" o 



2 p -Li 

3~* " F 



,-# 



where terms of order e" 3J? are neglected. Thus for large R 
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r * - |~S . (42) 

That is, £Jie quantity r dominating the bond in H2" 'is proportional to the 
overlap between the orbitals. At large R, this leads to a bond strength of 
the form 

r « - |-J? 2~ R ; (43) 

thus the bond energy decreases exponentially with internuclear distance. 

This simple relation between bonding does not hold for small R. We 
saw above that r is a niinimum (most negative) at R = 2ao, and the total 
energy is also a niinimum (bonding a maximum) around R = 2ao. On the 
other hand, the overlap continues to increase as R is decreased until S = 
1 at R = 0. 
2.1.5e Contragradience 

The above discussions indicate that the interference or exchange 
part of the kinetic energy dominates the bonding in H2 . This term is 
dominated by 

t* = fc[<(V X i) + (Vxr)> - S<(V Xi ) 2 >] . (44) 

which is large and negative in between the atoms. The region of space 
leading to negative Vxi ■ Vx r (and hence dominating the bond) is indicated 
for H 2 in Figure 14. 
2.1.5/ Historical Development 

H. Hellmann [after escaping from Hitler Germany into Russia in 
~1934 and before suddenly vanishing into Stalin Russia around 1937] was 
the first [ Z. PhysHc, 85, 180 (1933)] to suggest that bonding arises essen- 
tially from a decrease in kinetic energy. He suggested that the bond in 
Hg" results basically because the electron is allowed to delocalize over the 
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Figure 2-14. Contour plots of two hydrogen atomic orbitals for R = 2ao 
(contour increment 0.05 a.u.). The shaded region leads to negative values 
of VXi ' V;&. and hence to large contragradience. As a result, this region 
dominates the bonding. 
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region spanning two protons rather than just one. Using the uncertainty 
principle, he reasoned that a bigger "box" for the electron leads to a 
lower kinetic energy. Essentially the idea is as illustrated in Figure 15, 
where we see that for the <p g state the electron is distributed over a larger 
volume in Hi than in H atom. From the study of a p article-in- a-box, we 
know that the kinetic energy decreases as the box is made larger. Hence, 
because of a decrease in kinetic energy, the <p g state is expected to be 
stabilized with respect to H atoms. 

On the other hand, since the 9? u state has a node in the middle, the 
energy is just the same as if we had put the electron in either of two 
boxes, each of which is smaller than for H atom. This leads to an increase 
in the kinetic energy. 

Hellmann presented only very simple qualitative ideas and his view of 
bonding was largely ignored until Klaus Ruedenberg [ Rev. Mod. Phys., 34, 
326 (1962)] provided a more quantitative framework, showing (for specific 
cases) that interference terms resulting from the superposition of ampli- 
tudes leads to a significant decrease in the kinetic energy. Indeed, most 
workers before Ruedenberg argued that the bonding results from elec- 
trostatic interactions arising from increasing the density in the bond 
region. The development in this chapter is derived from a series of 
papers by C. W. Wilson, Jr., and "W. A. Goddard III [ Chem. Phys. Lett., 5, 45 
(1970); Theor. Chim. Acta, 26, 195, 211 (1972)]. Other somewhat related 
viewpoints have also been proposed: M. J. Feinberg and K. Ruedenberg, J. 
Chem. Phys., 54, 1495 (1971); M. J. Feinberg, K. Ruedenberg, and E. L. 
Mehler, Advan. Quant. Chem., 5, 28 (1970); R. F. Bader and A. D. Bau- 
draut, J. Chem. Phys., 49, 1653 (1968). 
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Figure 2-15. Illustration of the differences in the effective size of the box 
for the electron in the hydrogen atom and in the g and u states of E£. 
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2.2 The Molecular Orbital Description of H2 

We will now add a second electron to Hi to obtain the H 2 molecule. 
The simplest wavefunction for H 2 is to start with an electron in the best 
molecular orbital of H2" and to place a second electron in this p g -orbital. 
This leads to the molecular orbital (MO) wavefunction for H 2 , 

*£°(ri.ife) = ?g(ri)? g (r 2 ) (l) 

where 

n = Oft +Xr)/D g . (2a) 

and 

D g = V2(l + S) . (2b) 

With two electrons, the total wavefunction $(r lf r 2 ) must specify the proba- 
bility amplitude for electron 1 to have each possible value of its three 
coordinates (x^ y lt and z x , symbolized collectively as rj, and for electron 
2 to have each possible value of its three coordinates (x 2l y 2 , and z 2 , sym- 
bolized collectively as r 2 ). Thus, the wavefunction must be specified for 
all six simultaneous components of rj. and r 2 as in (1). 

First we will examine the meaning of the wavefunction (1). The total 
probability for electron 1 to be at some position r lt while electron 2 is 
simultaneously at some position r 2 is 

P( ril r 2 ) = l*" ^)! 8 = l^rOH^te)! 2 = P g ( ri )P g (r 2 ) . (3) 

This is just the product of the independent probabilities for electron 1 to 
be at position r x and electron 2 to be at position r 2 . Thus, the probability 
distribution for electron 1 is independent of electron 2.* Summarizing, a 
product wavefunction as in (1) implies that the electrons move indepen- 
dently of each other (no correlations in their motions) and vice versa. 



* Consider the analogous case of a red die (electron 1) and a green die (electron 
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In addition to using the <p g MO as in (2), we can construct wavefunc- 
tions of H 2 using the <p u MO, 



<Pu = (Xi -Xr)/^2(1-S) 
This leads to wavefunctions of the form 

*vg(1.2) = ?u(D?g(2) 

*gu(1.2) = <p g (l)<p u (2) 



(4) 



(5) 



(6) 



*uu(1.2) = «p u (l)^u(2) . (7) 

Since the <p u orbital is antibonding, the above wavefunctions of H 2 lead to 
much higher energies than (l) (except at large R), and we expect an 
energy level diagram as in Figure 16. 



£ 






9g9 u , ^ u ^g 



<t> $ 



Figure 2-16. Simple enei^gy diagram for MO wavefunctions of Hg. 

2.2. 1 Energies 

For H 2 we use the coordinate system of Figure 17. 



2). The probability of rolling a red 3 is 1/6 and the probability of rolling a green 5 
is 1/6 so that the total probability of getting both a red 3 and a green 5 is 

R~ X R~ = TfT ^ e < ^ ice are "^-dependent so that the probabilities multiply. 
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Figure 2-17. Coordinates for Hfe. 

Using the same conventions and assumptions as for H2" leads to the Ham- 
iltonian 

£(H 2 ) = h(l) + h(2) + ^-+ i- , (8) 

where is the Coulomb interaction between the two electrons, and 

r 12 

where 

h(i) = -HVf-j-- ±- (9) 

Lai Mri 

contains all terms depending only upon the coordinates of electron i. 
Consider now the energy of a product wavefunction 

*ab(1.2) = .*>.<l)*b(2) (10) 

and note that many two-electron integrals factor into products of one- 
electron integrals, e.g.,* 

<$ ab (l,2)|h(2)|§ ab (l,2)> = //d 3 r 1 d 3 r^ a *(l)^ b *(2)h(2)^ a ( 1)^(2) 

= [/d 3 r^a*(l)^(l)][/d 3 r 293b n2)h(2)^ b (2)] 

= <PalPa> <9> b ih|p b > (.11) 



* Note that the < > notation implies integration over however many electrons 
are in the wavefunction. 
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= <?blh|p b > 

since <<p & \<p a > = 1- In order to simplify the energy expressions, we will 
define 

The integral that does not factor is the one arising from the l/r 12 
terms in the Hamiltonian, which we will denote as 



Jij = < * a (i,2) 
< 



1 
ri8 



*n^)> - //«* A yi ' (1)i " (1)Pi>(3Vi(3) . (13) 

> r l2 



Note that 

Pi(l) = 9i*(l)^i(l) (14a) 

is the probablity density for finding the electron in orbital i at position r^ 
and 

Pj(2) = ^(2)*i(2) (14b) 

is the probability density for finding the electron in orbital j at position 
r 2 . Thus we can rewrite (13) as 

j, = //d^s^L. (is) 

But (15) is just the classical electrostatic interaction energy (the 
Coulomb energy) between the two charge distributions p s and p- y Thus we 
refer to J% as the Coulomb integral between orbitals & and <p y This term 
cannot be factored into a product of one-electron terms because of the 



ria = V(xj - xa) 2 + (yi - y 2 ) 2 + (zj - z 2 ) 2 
term. 

Since p{, p ]} and - L - are always positive, we see that the Coulomb 
integral is always positive 
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h * ° • (16) 

Using the above results we can write 

Egg = 2hgg + J gg + £- (17) 

E gu = Km + hgg + Jgu + g~ (IB) 

Eug = E gu (19) 

E uu = 2h uu + J uu + — . (20) 

For R = 2ao = 1.06 A (the Re for Hs? ), the MO wavefunctions lead to 

hgg = 

(21) 

hmi = 
Jgu = 

Juu = (22) 



Juu — 

1_ _ 

r 
resulting in 



Egg - 



Eim — 



Eug — E« u — (23 ) 



Thus the MO states are ordered as in Fig. 16. Note here that the Coulomb 
interactions are not negligible, but nonetheless, overall ordering of states 
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ean be predicted solely from considering the one-electron terms. 

2.2.2 Symmetries 

Before proceeding further in the discussion of H 2 , we will examine 
how symmetry can help us in sorting out the states. 

2.2.2a Inversion 

Starting with the configuration of particles in Figure 17 and invert- 
ing* the coordinates of electron 1 leads to Figure 18a, 

e2 

y* 

// r 12 \ r b2 

," \ 

a y/_„ R \ a ^ r ^ b 



v 



' ^ <n ^*~ / 



r bl '/ ~T~ r / \ al ^-^ ' 

'/ ^-'^al r bl / s^~- /r. 2 

el w el ^ - ^. \ / 

(b) n e2 

Figure 2-18. 

with electron 1 at a distance of rbi from nucleus a and r al from b). Just as 
for Ha", this does not change the nuclear attraction terms if the nuclei 
have identical charges. Even so, the total potential energy is changed 
since the distance between electrons 1 and 2 is changed. Thus, in order 
to preserve the same potential energy, we must simultaneously invert 
the coordinates of both electrons, leading to Figure 18. Thus we define 
the inversion operator, T, for H 2 as 



* Recall that the center of inversion is the bond midpoint. 



Xi -* "Xi 




X2 - -x 2 


yi -* -yi 


and 


y 2 -> -y 2 


Zl -* "Zi 




Z2 "► -Z 2 
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(24a) 

Z x -» -Zi Z2 -► -Z2 

or 

ri -* -ri and r2 -* -r 2 , (24b) 

and the total Hamiltonian is invariant under this inversion, 

T£(r lt r 2 ) = H(-r lt -r 2 ) = H(r lt r 2 ) . (25) 

The exact wavefunction for H 2 is obtained by solving 

H(l,2)$(l,2) = E*(1.2) (26) 

(don't panic; we won't try this yet). Since #(1,2) is invariant under the 
inversion T (25), we find that (26) implies 

£(I$) = E(I$) (27) 

and hence (just as for H 2 + ), the exact eigenstates of H 2 are either g or u, 

T* g (1.2) = + * g (1.2) 

(28) 
T$ u (l,2) = -* u (1.2) . 

Applying the T operator to the MO wavef unctions (1), (5)-(7) leads to 
the conclusion that 

*gg and $ uu are g states 

(29) 
and 

$ gu and $ug are u states. 

Since a u wavefunction must always have a nodal plane, somewhere we 
expect the ground state of H 2 to be g, just as indicated by Figure 16. 

2.2.2b PermutatTJonal Symmetry 
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The Hamiltonian for H 2 is unchanged if we renumber the electrons so 
that electron 1 becomes electron 2 and vice versa, that is, 

£(2,1) = #(1,2) . (30) 

To discuss such symmetries we define the transposition operator r 12 as 

ri -> r 2 and r 2 -* rj . (31) 

Thus, starting with the exact wavef unction (26), applying t 12 to both sides, 
and using (30), we find that 

£(t 12 $) = E(t 18 *) . (32) 

Since (t 12 ; transposes the electrons twice, taking us back to the original 
starting point, 

(T i2 ) 12 = t 12 t 12 = e (33) 

[where e signifies doing nothing, i.e., the unity or einheit (German) opera- 
tor], we can show (just as for inversion) that the exact eigenstates have 



the behavior 



r 12 $ s (l,2) = $ a (2,l) = + $ s (l,2) 



(33) 



or 



r 12 $ a (l,2) = $ a (2,l) = -$ a (l,2) 

under transposition. Examining the wavefunctions (1), (5)-(7), we see 
that 



Tig[Wl)p,(2)] = [dlk(2)] 



symmetric 



(35a) 



Ti 8 [sP g (lW2)] = [«»u(l)sPg(2)] 

Ti 2 [?u(D*g(2)] = [?g(l)?u(2)] 

Ti 2 [^(l)^ u (2)] = [^(1)^(2)] 



no symmetry 



symmetric 



(35b) 
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Thus the g states, $ gg and $ uu , have the proper permutational symmetry, 
but the u states, $gu and $ug, do not. In the next section we will fix up this 
problem with the u states. 

2.2.2c The u States 

Combining the u wavef unctions as follows* 

3 $u = («Wx-puPg) (36) 

2 *u = (pgPu'+PuPg) (37) 



leads to 



T12 3 *u = (<Pu? g - PgPu) = ~ 3 $ u (35c) 



T12 l K = (PuPg + <P e 9xd = + *u (35d) 

so that the combinations (36) and (3?) have proper permutational sym- 
metry. 

Now we must examine the physics behind the combinations (36) and 
(37). Aside from mathematical analyses indicating that the wavefunc- 
tions should have permutational symmetry, we also want to determine 
why one combination is favored in terms of achieving a lower energy. 

In order to carry out such an analysis, we will plot the two-electron 
wavefunction for the case where both electrons 1 and 2 are along the 
bond axis (z). In order to show the relative locations of both electrons, 
we will let the z coordinate of electron 1 be the ordinate (z^ and the z 
coordinate of electron 2 be the abscissa (z 2 ). This is indicated in Figure 
19 where some special points are indicated. 



* The notation will become clear when we discuss spin in Chapter 4. 
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Figure 2-19. Coordinates showing simultaneous positions for Zj and Zj. 

In order to see how to plot the wavef unction on Figure 19, consider the 
orbitals ip g (zj and <p u (z 2 ) in Figure 20ab. 



(a) * g (Zi) * 




1 h»Z, 
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(c) * gu (Zi,Z 2 ) 
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Figure 2-20. 
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Multiplying these orbitals leads to the wavefunction $gu(zi, z 2 ) with ampli- 
tudes at various points, as given in Figure 20c. Rather than listing 
numbers as in Figure 20c, we will draw contours of equal amplitude as 
indicated in Figure 20d, where solid lines indicate positive amplitude, dot- 
ted lines indicate negative amplitude, and long dashes indicate zero 
amplitude- Here one sees that the maximum positive amplitudes occur 
for z 2 = a and zj = a or b, while the amplitude remains large or positive 
between the above two points. The maximum negative amplitude occurs 
for z 2 = b and z l = a or b and remains large between the above two points. 
For z 2 = 0, the total wavefunction is zero, independent of z^ The reader 



should practice constructing .20d directly from Figures 20a and 20b 
(without going through Fig. 20c). 

In Figure 21 we compare the two-electron wavefunctions of $ gu , $ug, 
3 $ u , and ^u. All cases involve a single nodal plane, and if there were no 
electron-electron interactions, all these wavefunctions would have the 
same total energy. * The difference lies in the electron-electron interac- 
tions. In Figure 21 e we see that the 3 $ u wavefunction is zero along the 
line with zi = z 2 , whereas in Figure 2 If we see that the largest values 
(positive and negative) of l $ u occur for z x = z 2 . Since the electron repul- 
sion term 

e 2 

ri2 

is large and repulsive when the electrons are close, we see that 3 <$ u is 
favored and ! $ u is disfavored. Indeed, considering all possible combina- 
tions of $gu and $ug, the one with the lowest electron repulsion is just 3 $ u . 



* Note from (18) that both u states have the same total energy and that any com- 
bination of these states leads to the same one-electron energy, h^ + h gg . , 
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(a) 4> t 




(b) *, 




( c > *gu (1,2) = ^g (1) *u (2 > 




(d) #^(1,2) = <P u (l)<t> g (2) 



b -- 



a -- 




(e) 3 * u (1 > 2 > = Vu"^g 



(f) l » u (l,2) = <*> g * u + * u * g 





Figure 2-21. 
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Thus the energy diagram for the MO wavefunctions of H 2 becomes as 
in Figure 22. 



■uu 

■single excitation 




double excitation p u p u = $ x 

PgPu + 9u<Pg - ^u 

Pg?u - <Pu9g - 3 $u 

ground state ? g p g = <$ gg 

Figure 2-22, 

2. 2. 2d Quantitative Aspects of the Energies for the u States 

Now that we see the physics behind why the 3 $ u wavefunction (36) is 
the best u state, we will examine the quantitative energy expression. 

Ef 3$ ) = < 3 $ul#l 3 $u> = < g u-ugi#lgu-ug> 
1 < 3 § U \H\H U > <gu - ug | gu - ug> 



(38) 



Since 



_ 2<gulff lgu-ug> _ <gu|ff |gu> - <gu|#!ug> 
2<gu|gu-ug> <gu|gu> - <gu|ug> 



<gu]ug> = <glu><u|g> = 
<gu|h(l)jug> = <gjh|u><u|g> = 
<gu|h(2)|ug> = <g|u><u|h|g> = 



we are left with 



where 



E( 3 $ u ) = E gu -<gu|//|ug> , (39) 



< 

<gu|#|ug> = <gu 
< 


1 

r 12 


> 

ug> 
> 




„ r ,* ^*(l)^u(D^u*(8)^(2) _ 


Kgu 


^12 







(40) 
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where this two-electron term is called the exchange integral between 
orbitals xp g and <p u . The net result is 

E( 3 $ u ) = Egu - K gu 

(41) 
E( J $ U ) = Egu + Kgu . 

Since the previous section showed that 

EO^) > E( 3 $ u ) , 

we see that the exchange integral must be positive, 

K^ =* . (42) 

What is the physical significance of K gu ? It tells us how much the 
two-electron energy changes when we go from the wavefunctions in Fig- 
ure 21 c or 2 Id to the wavefunctions in Figure 21e or Z If. Thus the K^ is 
the quantitative representation of our earlier argument that 3 <$ u has a 
better two-electron energy than l $ u . It is better by precisely 2Kgu. 

Another way to look at this is to substitute (17) into (40), leading to 

E( 3 $ u ) = [h gg + h uu + ^ + [J gu - K gu ] 

(43) 
E( 1 $ u ) = [hgg + h uu +■ ^ + [Jgu + K^] . 

Here J^ - K gu is the total two-electron energy of Figure 2 If, while J gu + K^ 
is the total two-electron energy of Figure 21e, and J^ is the two-electron 
energy of Figure 21c and 2 Id. The two-electron energy of 3 $ u can also be 
written as 



< 3 $u 



s u g u ~~ 



8 , > 

Xd.2) 



< 3 K\ 3 K> JJ dr * dr * ri2 • {44) 



where* 



The ^- comes from < 3 $ u 1 3 $u> = 2- 
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3 Pu(L2) = |- 3 *u*(1.2) 3 $ u (l,2) . 

Since 3 P U is the absolute square of 3 # u , it is positive for all possible values 
of r { and r 2 . Since the integral in (44) is always positive, we see that the 
total two-electron energy of the 3 $ u state must be positive 

Jgu ~ Kgu ^ , (45) 

and hence 

Jgu ^ Kgu (46) 

the exchange integral is always less than the Coulomb integral. Combin- 
ing (46) with (16) leads then to* 

J^ > Kgu > . (47) 

2.2.3 Potential Curves 

So far we have discussed the MO wavef unction assuming that the 
bonding orbital <p g is much better than the antibonding orbital <p u . This is 
true for shorter internuclear distances R but does not remain true as the 
bond is broken. Thus, in Figure 23 we compare the energy of the MO 
wavefunction $ gg with the exact energy for the ground state of H 2 . This 
MO wavefunction leads to a good value for the bond length but a very bad 
description of the processes of breaking the bond. 

The origin of this problem can be seen by substituting the atomic 
orbital describtion of the MO (2) into the MO wavefunction (1), leading to 

$jf(1.2) = N[XiXl + XrXr + XlXr + XrXl ] 

(48) 

= 'covalent + * ionic • 

where N = [2(1 + S)]" 1 and 



This relation is true for any pairs of orbitals, as shown in Appendix 2-C. 
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$covalent = HXlXr + XrXl) (49) 

^ionic = ^(XlXl + XrXr) ■ (50) 

At very large R, the exact wavefunction will have one electron near the 
left proton and one at the right, as in (49), which we will refer to as the 
covalent part of the wavefunction. The other terms of (48) have both 
electrons near one proton and none near the other, thus an ionic 
wavefunction. At R = », these ionic terms lead to the energy of H" and H* 
rather than the energy of two hydrogen atoms. Since the MO wavefunc- 
tion must have equal covalent and ionic contributions, it yields terrible 
energies for large R. 

The basic problem with the MO wavefunction is that both electrons 
are in the same p g orbital and hence each electron has an equal probabil- 
ity of being on either center, regardless of the instantaneous location of 
the other electron. In the exact wavefunction, the motions of the elec- 
trons tend to be correlated so that if one electron is on the left, the other 
tends to be on the right. This correlation is necessarily ignored in the MO 
wavefunction, and the resulting error is often referred to as the correla- 
tion error. For small R, the two centers are close to each other and this 
neglect of correlation is not so important. At R = «, however, the correla- 
tion of electrons is of paramount importance and neglect of correlation 
leads to ludicrously poor wavef unctions. 

In the next section we will discuss a simple wavefunction, the valence 
bond wavefunction, that eliminates this problem of describing large R. 
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2.3 The Valence Bond Description of H2 

2.3. 1 The Covalent States 

We will now reexamine the ground state of H 2 molecule. However, 
rather than the approach of the above section (plopping electrons one by 
one into the orbitals of H 2 + ). we will instead start with the exact wavef unc- 
tion at R = 00. This, of course, consists of two hydrogen atoms infinitely 
far apart, say electron 1 on the left and electron 2 on the right as in Fig- 
ure 24. 

el # e2 

* \ 

* 1 

•'t. ' 

/ la \ r 2b 



a 



\ 



R = 00^, 



■** 



Figure 2-24. 

The wavef unction for Figure 24 is 

*a(ri,r 2 ) = Xi(ri)Xr(r2) . (l) 

where 

XiM = N2^ 

(2) 

(and N is the normalization factor). 

This wavef unction $„ says that the probability of electron 1 being at a 
particular position is independent of where electron 2 is and vice versa 
(since the atoms are infinitely far apart, the electrons should not be 
influenced by each other). 

There is a second wavefunction that is just as good (or as bad) as $ 
in (1), namely, 
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(a) # a (l f 2) = Xf (l)X r (2) (b) * b (l,2) = X ,(l)Xi(2) 



b -- 



a -- 




(c) * g = * a + * b 

= X.«X r "*■ XyXj^ 





(d) * u = * a -, * b 



= X *x- • x-x. 
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Figure 2-25. 
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*b(ri.r 2 ) = Xr(ri)*i(ife) , (3) 

where the electrons have been interchanged. This wavefunction $ 6 is 
different from $ since electron 1 is on the opposite sides of the universe; 
however, the energies of $ 6 and $ must be the same (since electrons 1 
and 2 have the same properties). 

"We will find it useful to combine $ and $& into two new wavefunc- 
tions, 

*,(1,2) = *• (1.2) + ** (1.2) = *(l)Xr(2)+Xr(DXi(2) (4a) 

*u(1.2) = *.<1.2)-#*(1.2) = Xi(l)Xr(2)-Xr(l)Xi(2) (4b) 

(unnormalized) because at finite R these are the optimum wavef unctions. 
Before examining the energies, we need to understand how to think about 
the relative locations of the electrons in these wavef unctions. 

In Figure 25, we plot the four wavef unctions, $ a , $ b , $ g , and $ u . Here 
we see that $ u has a nodal plane (coresponding to z l = z 2 ) while $ ff does 
not. Indeed, along the line between the two peaks in Figure 25c, we see 
that the gradient of the $ y wavefunction is smaller than that of $ tt or $ b , 
while the gradient of the $ u wavefunction is larger. This decrease in the 
gradient of $ g (and increase for $ w ) depends upon R with a bigger effect 
for smaller R; thus, based on kinetic energy, we would expect that $ g is 
bonding and $ u is antibonding, and indeed this is the case, as shown in 
Figure 26. 

23-2 Symmetry 

Inversion of the coordinates of the electrons r al ^-rbi and r a2 4>rb 2 
leads to [see (2)] 

Xi to) «» Xrfa) • 
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Figure 2-26. The energies Eg and Ea for the valence bond wavefunctions of 

H2. 



LU 

tr 

h- 

< 
x 

o 
tr 

LU 
LU 

LU 
> 

< 

-J 
LU 
CC 




2.0 4.0 6.0 8.0 

NTERNUCLEAR DISTANCE (BOHR) 



10.0 



Figure 2-27. The energy of the MO wavefunction for the ground state of Kb 
with comparison to the VB and exact energies. 
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Consequently, 



(5) 



7*0(^x2) = Xr(ri)xi(r 2 ) = <Mri,r 2 ) 
or 

As a result, 

I$ g = /(* B + $„) = ($ 6 + $ a ) = <£ ? 

I# u = /($„ -$«,) = * b -#. = -$ u , 

and we see that $ ff and $ u are indeed gerade and ungerade, respectively. 

2.3.3 Comparison of YB and MO Wavef unctions 

2.3.3a Ground State 

The MO wavefunction is (ignoring normalization) 

*ff(l.2) = Pg(l)?g(2) = [XlXr + XrXl]+[XlXl +XtXt\ . 

whereas the VB wavefunction is 

^(1.2) = [XlXr + XrXl] ■ 

The energies for these wavefunctions are compared in Figure 27 where we 
see that the VB is always better but that the difference becomes negligi- 
ble for small R. 
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Figure 2-28. 

The wavefunctions are compared in Figure 28, showing graphically 
how the VB wavefunction has smaller probability of having Zx = z 2 , leading 
to lower electron repulsion energies. On the other hand, the MO 
wavefunction is smoother, leading to smaller kinetic energies. For nor- 
mal bond distances, the electron repulsion effects dominate so that the 
VB wavefunction is better. However, for very short R, the kinetic energy 
becomes dominant so that the MO and VB wavefunctions lead to nearly 
identical total energies. 

2.3.3b The u States 

Expanding the MO's in terms of AO's (and ignoring normalization) 
leads to 

*gu = (I + r)(£ - r) = 11 + rt - tr - rr 



*ug = {I - r)(l + r) = U - rl + Ir - rr 



Thus 



$u = *«-*«* = 2(r* -ir) = $™ 



(6) 
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**« = gu + ug = 2(11 -rr) = $J 0N . (7) 

That is, the first excited state in the MO description, 3 $ u , is identical to 
the first excited state in the VB description. Both describe a covalent 
repulsive state that separates to two free H atoms, as indicated in Figure 
29. 

The second state, ^ u , leads for R = « to the wavef unctions xiXi 
corresponding to 

H" H + 

and XtXt corresponding to 

H + H~ . 

Thus we refer to this state as the ionic state. The energy curve is shown 
in Figure 29. 

2.3.3c The Second g State 

In the MO description* 

$£? = (I - r)(l - r) = [(U + rr) - (It + rl)]/ 2(1 - S) 



(8) 



$H° = (* + r )(* + r) = [(U + rr) + (Zr + rl)]/ 2(1 + S) . 
In contrast, the VB description leads to* 



*P = (lr + rl)/^2(l +S 2 ) (9) 

and 



*f N = (U + rr) / V2(i + S 2 ) 
for the covalent and ionic g states. 

The connections between these states are 



g • V 2(i + s 5 ) l U + b;$gg U b)?m 



* Now we include normalization factors. 
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Figure 2-29. Energies for the states of H2, using atomic orbitals 

ft" = 1-0). 
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(10) 



(1 + 5) 



V2(i + S 2 ) 
where 



*gg A *uu 



^ = rff- <") 



Thus we can fix up the MO wavefunction so that it behaves like the VB 
wavef unction by mixing together the $^° and $2? wavef unctions. This is 
related to the configuration interaction (CI) wavefunction, as discussed 
below. For large R, S = so that X = 1, whereas for R = 1.6ao = 0.8 A , S = 
0.7, leading to X = 0.18. Thus the CI is most important at larger R. 

For a more general description of these states, we would consider the 
wavefunction to have the form 

*S = Ci §J|° + C 2 *JSF (12a) 

or* 

*g = Di^+DaSf* , (12b) 

and choose the coefficients that lead to the best energy. This is called 
the configuration interaction or CI wavefunction and leads to the results 
shown in Figure 29. The excited g state, $g, can also be taken to have the 
form (12); however, it must be orthogonal to *fj, leading to 

<$g|*g> = . 

The overlap between the covalent and ionic g states is 

<$ ion, $ co V> = <U+rr\lr + rl> (13) 

gig 2(1 + S 2 ) v ' 

2S 



1 + S< 



* Equations (12a) and (12b) lead to an equivalent total wavefunction, as can be 
seen by comparing (8) and (9). 
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Thus, for S = 0.7 (R = 1.6ao = 0.8 A ) this overlap is 0.95, demonstrating 
just how similar are the ionic and covalent wavefunctions for small R. This 
creates a problem in describing the excited g state. The VB wavefunction 
is a close approximation to the $g wavefunction; however, except for S * 
0, the §[ 0N wavefunction is not a good approximation to the excited state, 
$.§. Instead, we must orthogonalize $j 0N to <$f ov , leading to new nodal 
planes and a much higher energy. This explains why the $£J state is 
always above the $i 0N state. Based on the nodal theorem, we would 
expect that §[ 0N (which has no nodal planes) would have a lower energy 
than $i 0N , and it does. However, the only nodeless state is the ground 
state $ lg which mixes whatever combination of fcf ^ and $] 0N gives the best 
energy. The excited g state, $g, necessarily has nodal surfaces since it 
must be orthogonal to the ground state. The result is that the ionic g 
state, $g, is always above the ionic u state. 

2.3.4 Quantitative Analysis of Bonding in Hg 

We will analyze the energies of the VB wavefunctions for the g and u 
states of H 2 in a manner very similar to that used for the LCAO wavefunc- 
tion of the g and u states of H 2 + . 

First we consider the energy of the simple product wavefunction 

$ c <(1.2) = »(DXr(3) = XiXr . (14) 

which is just part of the wavefunction for the g and u states, (4a) and 
(4b). We will refer to this wavefunction as the classical wavefunction and 
the energy 

<$ ci \H j$ ci> 



E ci = 



(15) 



- h« + hrr ■*- J(r + p~ 
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as the classical energy. 

The total energy of H 2 differs from the classical energy due to the 
presence of a second term in the wavefunctions (4a) and (4b). The 
second term has the electrons interchanged (exchanged) and hence is 
called the exchange term. 

The effect of the exchange term in the wavef unction, say (4a), is to 
change the energy from E ci to E g . We will refer to this change in energy 
as the exchange energy EJ, so that 



E g s E ci + E£ 



E = E ci + E *• 



(16) 



In Figure 30 we show the behavior of these quantities with R. Just as for 
H 2 + , we see that it is the exchange term that dominates the bonding 
energy. Indeed, we will find that it is the one-electron part of the 
exchange terms that provides the dominant interaction, just as for H 2 + . 
First we consider the detailed form of the various energy quantities. 

2. 3. 4a Analysis of E* 

At large R, the one-electron term 



hii = <Xi 
< 






> < 

Xi > + <Xi 

> < 



1_ 



> 

Xi> ■ 
> 



J 



atomic energy 



penetration term 



has the form 



Neglecting terms of order e _2R , and the Coulomb term has the form 



(17) 



>lT 
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1.0 2.0 3.0 4.0 5.0 6.0 7.0 
INTERNUCLEAR DISTANCE (BOHR) 



8.0 



Figure 2-30. The classical (E*), exchange (E*), and total (E) energies for 
the VB wavefunctions of Hfe. Note that each energy is references to the 
value for R = «, that is, E^) = JL(«) = EuH = -1.0 and E*H = EJ(«) = 
0.0. 








.0 2.0 3.0 4.0 5.0 6.0 7.0 
INTERNUCLEAR DISTANCE (BOHR) 



8.0 



Figure 2-31. The total exchange term (r) and the one-electron (tj and 
two-electron (r 2 ) parts. The potential (t t ) and kinetic (T t ) parts of fi are 
also shown. All results are for !%>. 
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( neglecting terms of order e~ 2R ). Thus the classical term is just twice the 
energy of an H atom (£i 3 ), 

E ci w 2e l8 

(neglecting terms of order e' 2 ®), with no net Coulomb terms. Including 
the additional penetration terms leads to 



E c< = 2e ls -H 



F + i + i*-F R2 



e" 2R , (18) 



corresponding to the interpenetration of the two atomic electron clouds. 
Although negative for R > 1.4a , this quantity is small, as shown in Figure 
30. Thus the bonding of H 2 cannot be explained as due to penetration of 
the charge clouds of two hydrogen atoms. 

2.3.4b The Exchange Terms 

Now we consider the energy of $| OT , (4a), 

E g = <* g |i?|* g >/<* g iv . 

By symmetry 

<<S g |$ g > = <XlXrl*g> + < XrXil$g> 
= 2<XiXrl*g> 

and 

<* g \H\* g > = <XiXT\H\$ g > + <XrXi\H\$ g > 

= 2<xiX r \H\* g > . 
Hence 

E _ <XlXgl#l$g> 

g ' <XlXr\$g> 

Evaluating the individual terms, we find 
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<XlXr\$g> = 1 + <XilXr><XrlXt> 

= 1 + S 2 
and 

<XlXr\8\*g> = <XlXr\H\XlXr>+ <XlXr\8 \XrXl> = E ci + $ , 

where 

$ = <XiXx\H\XxXi> (19) 

is referred to as the VB exchange term. Thus 

Eg = f7^= E ci +E| , (20) 

where the exchange energy is 

g (i + s 2 ) { * l) 

The same analysis for the Sj 8 wavefunction (4b) leads to 



= £Ct -B = gci + 



Eu = "_ s r = E°< + E* , (22) 



where 



^ (1-S 2 ) ' K 6) 

The close relationship between E£ and E£ is emphasized by defining* 

r = (^-S 2 ^) (24) 

so that 

E £ = /, ? c.2x (25) 

g (1 + S 2 ) v ' 



E5 = jf^- • (26) 
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From (25) and (26) the energy separation between the g and u states is 



E ff - E u - 



Zt 



,1 -S 4 



(27) 



From the nodal theorem, E g < Eu, and hence 

t < 
since S < 1. 

These results for H 2 are quite analogous to the case of H 2 + where the 
<p g and (p u state have energies 

1 



S, = E cl + 



1 + S 



(28) 



d, = £ 



cJ 



1 -S 



(29) 



£ g £ u - 



2T 



(i - s 2 ) 



with 



(30) 



e« = ha + f- 



(31) 



T = ht, - Shu 

2.3.4c Analysis of g£ 

The components of ^f (19) are 

<XiXr|h(l)|XrXi> = <Xi|k|Xr><XrlXl> 



(32) 



= hirS 



<XiXr|h(2)|XrXi> = <XiiXr><Xr|h|Xi> 



(33) 



* We use T here in order to distinguish this quantity for H2 from the r of Hi- 
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= Sh, 



•rt 



<XiXr\^HXrXt> = /dT»*(l)^(l)/dT 8 ^-Xp*(2)xi(2) 

= Kir 
<XlXr\ ^-\XrXl> = -pT ■ 

where S is the overlap and Ki r is referred to as an exchange integral. 
(Note the distinction between £F, the VB exchange term, and IQ r) the 
exchange integral.) Thus, 



S 2 



f = 2Sfa, P + K, r + |- 



2. 3. 4d Analysis of r 

Using E ct from (15) in (24), we find that 

f = f^Tz , 



(34) 



(35) 



where 



f x = [2Sh, r - S 2 (h.i +h n .)] = 2S[h, r -Shu] = 2Sr 



(36) 



T 2 ~ Kr " S Jf r ] 



(37) 



are the one- and two-electron parts, respectively.* These quantities are 
plotted in Figure 31 where we see thatt 2 has a smaller magnitude that tv 
Thus we can write 



t w T, 



Comparing (27) and (31a) we see that 



f l = 2St 



(38) 



(39) 



* The —terms cancel. 
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that is, the one- electronpartofrfor H 2 is related directly to the r of Hi, 
leading to 



r = f l = 2St , 



(40) 



where t is the quantity for H 2 + . Thus, for H 8 the bonding energy is deter- 
mined by 



ei - flw - *s 



(41) 



EJ * - . 3S L = EJ , 



1 -S a 



whereas for Hi it is determined by 



(42) 



Go — 



g 1 + s 



(43) 



E X = — 



1 -S ' 
These quantities are compared in Figure 32. 



(44) 
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Figure 2-32. The exchange energies (EJand EJ) for Hg" and Efe. In addition, 
the one-electron approximation E* to the E* for Ha is shown for each state. 
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2.3.42. Analysis off 1 = 2Sr 

Since the quantities f 1 and r dominating the bonding in H 2 and H 2 + are 



related, 



Ti = 2Sr 



it is well to examine the reasons for these relations. The wavefunctions 
for the bonding states of H 2 and Ha" are sketched in Figure 33ab. In both 
cases the kinetic energy is decreased from that in the classical wavefunc- 
tions, Figure 33cd. The decrease in the kinetic energy for electron 1 is 
obtained by examining the gradients in the vertical direction (ordinate) 
of Figure 33abcd. Here we see that H 2 + leads to a larger decrease that H 2 . 
Thus the contribution is St for H 2 but r for H^. However, for H 2 there is a 
second electron (number 2) that has a similar decrease. Thus, for H 2 the 
net is 2St as compared to r for H 2 . 



(a) H 2 * g = (XjXj. + XpXj) 




(b) HfJ <f> g = ( Xf + xp 




(c) H2 * 



c£ 




(d) Hj * ct 




Figure 2-33. 
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2.3. 4 f Comparison of Bonding in H 2 and H 2 

Although the bonding energies of H 2 and H 2 + are both determined by 
t, we see from (4l)-(44) that the value of the overlap S also plays an 
important role. From 2. 1 the form of t at large R is 

t * -(2/R)S , 

hence 

r l = 2Sr w -(4/R)S 2 . (45) 

Thus the bonding in H 2 + is proportional to S, but the bonding in H 2 is pro- 
portional to the square of S. 

At R = 1.6a , the value of S is 

S = 0.7, 

and hence 

Ha: EJ = 0.94 r 

EjJ = -2.75 r 

(46) 
H 2 + : e| = 0.67 r 

s£ = -3.33 t . 

Thus the g state of H 2 should have a bond energy about 50% larger than 
the g state of H 2 , while the u state of H 2 should be 17% less repulsive than 
the u state of H 2 . In addition, we see that the u state of H 2 should be 
about three times as repulsive as the g state is attractive. 
At R = 3ao, the overlap is S = 0. 1 and hence we obtain 
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H 2 : E£ = 0.20 t 

E* = -0.20 t 

(47) 
H 2 + : £| = 0.91 t 

£* = -1.11 T . 

Thus, at this large R, the g and u states of H 2 are five times as attractive 
or repulsive as the g and u states of H 2 . That is, at Large R, the one- 
electron bond is much stronger than the two- electron bond! This 
difference in relative bond strengths of H 2 and H 2 + for small and large R 
just results from the overlap term S that automatically arises in the 
exchange of a two-electron wavefunction. 

In the limit that S = 1, we have* 
H 2 : E| = r 

H 2 + : ef = fcr , 

leading to an H 2 bond twice that of H 2 + , the commonly expected result. 

PROBLEM: Experimentally (and theoretically), the molecules Li 2 , Na^, K^ , 
with one-electron bonds, are found to have stronger bonds than the 
corresponding molecules L^, Naa, K2 having two-electron bonds. Explain 
the origin of this effect. 



Actually, S = 1 implies R = 0, which in turn implies T = 0. 



2-A-l 
APPENDIX 2-A: ENERGY QUANTITIES FOR H 2 + 

We will consider an atomic orbital of the form 

X=^A e' ?r (1) 

centered at each of the two nuclei of H^. The coordinates are indicated in 
Figure 34. 

e 

s 

/ / 

s / 

*,' I b 

a*— "-^ 

(a) Figure 2-34. 

With ? = .1 the orbitals (1) correspond to hydrogen Is orbitals on each center. 

First we evaluate the atomic integrals for general ?, then the new 

energy quantities occurring in H 2 + . 

2-A.l Atomic Energy Quantities 

The norm of X is 2^ «, 

<X I X> = (C%) J sinBdO / d<p J r*dr e~ 2 ^ 




All 



00 



45* J r 2 dre- 2?r 



OO 



2 -1 



= i / P e" p dp 



= 1, 
where we set p = 2£r and used 



00 



/ p e~ p dp = m! . 
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Similarly the atomic potential energy is 
<x|-ilx> =-4? 3 /r 2 dre- 2?r 



= -? . 
and the kinetic energy is 

<x|-iv|x> = i<|v x f> = ie 2 <x|x> 



(2) 



(3) 



= K . 

2-A.2 Elliptic Coordinates 

In evaluating the energy quantities for diatomic molecules, it is 
convenient to use elliptic coordinates 

« = (r a + r b )/R 

71 = (r a " r b )/R (4 > 

(p = azimuthal angle about the z axis (measured 
from the xz plane) 
in place of the cylindrical polar coordinates p, <p 9 z (see Fig. 34b) or 
spherical coordinates r, , <p. (See Fig. 34b.) The geometric condition 
defining an ellipse is that the sum of the distances to the two foci is a constant, 
and hence each curve of constant £ corresponds to an ellipse. Similarly 
from the defining condition for ja hyperbola, each surface of constant r) correspoi^i 
to a hyperbola. The range of the elliptic coordinates is 

-l^Tj^ + l (5) 

< <p < 2ir 
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The volume increments in the various coordinate systems are 
dr = dx dy dz 
dr = p dp dfp dz 

2 

dr = r sine dr d0 dp 
dr = iR 3 (| 2 - V 2 ) d|dr/d<p 
The latter relationship can be derived from 



(6) 



x=|R\|(4 2 - 1)(1 - rf) cos<p 



(7) 



y = iR\|(« - 1)(1 -*]) sincp 
z=|R^ 

3 2 2 

since the Jacobian |R (I - r\ ) is just the determinant of the derivative 
matrix. 

From Fig. 34b we see that 

2 2 , R 2 

r b = P + ( z - 2") 

2 2 / Rx 2 

r a = p + (z + y ) 

2 2 2 

r = P + z , 



and hence 



and 



i/ 2 2 n 2 /Rx 2 
3(r a + r b )=r + ( y ) 



2 2 

r - r, = 2zR . 
a b 



From (4) we find 



2 2 2 2 2 

i(l + iJ ) = (r a + r b )/R 
In = (r * - r. 2 )/R 2 



2 2 2 

5 - „ = 4r a r b /R , 



(9) 
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and hence 

2 , 2 r 2 2 .. 
r = iR [ I + V - 1] • (8) 

These relations will be useful in the next section. 

In evaluating integrals over £ the following integral will be 

useful 

00 fvx m k 

r m -ax, m!e" « v (« x o) - a /„ \ m^ 

v a k=0 

For x = 1 and m = 0, 1, and 2 this becomes 

A (<*)=^e- ff (11) 

AxM^e-^l + a] , (12) 

A 2 (a)=| 3 e- a [l + a + iof] (13) 

2-A.3 Diatomic Energy Quantities 
2 -A. 3a Overlap Integral 

First we evaluate the overlap integral, 

s = <x je lx r ). 

Using (1), (6), (11) and (13), this becomes 

3 3 °° ■ +1 

S = Jdr X £ X r = (V ) (V ) (2») / « / * «* " i?V 2CR{ 

1 -1 

= e- ?R [l + (SR) + i(?R) 2 ] (14) 

2-A.3b The Exchange Potential Energy Term 

There are two terms involved in evaluating the potential energy 
of H 2 , the exchange terms 



where 



and 
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V lr s <X £ Jv J |x r > , (15) 

(16) 



(17) 





V 


1 1 
" " r a " r b 






V 






action 


we 


evaluate V* . 




First 


we 


convert v to ellipt: 


ic coordinates 




a. 

V 


_ (r a + r b } _ 
Vb 


44 




2 2 

R(l v) 



using (4) and (9). Combining (18) with (6) leads to 

vdr = - iR 2 | d£ d?7 dp , 
and integrating over (p we obtain 

V £r = ( -^ )( ^- } (27r) / |d |e " ?RS / d7? 

1 -1 

= -2?e" ?R (l+ ?R) . 

A. 3c The Penetration Term 



(18) 



(19) 



The other potential energy term V ^ has two parts 

v 11 =<x 1 l4lv + <x il-ilV' ; (20 > 

a d 

the first of which 

V^ - <X,l4 Ix,) (2!) 

a 
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is a one center integral (involving only the left nucleus) and the other 
of which 

involves two centers. This second term is the Coulomb interaction 

I |2 

between the spherically symmetric charge distribution p- = | | 
centered on the left nucleus with the charge centered at the right 
nucleus and is referred to as the penetration integral. 
As shown in equation 8 of Appendix B 

vg=A + B, 



where 



-2?r a 



B = -4» (i-) / r a dr a e 
R 



A R 



and 



3 R . _2£r 

a 



Q=4»(l,)/ r'dre 



Using 



a a 




R 

/ p"dre- ar = 4 [ 1 - (1 + «R + ^ 2 R 2 )e" aR l , (23) 

rf 



we obtain 



Q = [ 1 - (1 + 2CR + 2? 2 R 2 )e" 2?R ] . (24) 

Using (10) we obtain 

B = ^e" 2 ^ R (l + 2?R) , 



and hence 



2-A-7 



V S = " k (1 ' e " 2CR ) + ^ e " 2CR • (25) 



2-A.3d The Kinetic Energy Term 

The two-center kinetic energy integral is 

which from Appendix 1-A becomes 

= i/dr(Vx 1 - V X ). 



Since 



we obtain 



(26) 



r b = Vx + y + (z - R/2) 

/2 2 " 2 

r a = Vx +y + (z + R/2) , 

ve" a = - | e a [xl + yg + (z + |R)t ] 
a x y z 

v e " r b = .I e " rb [4 x + yg + (z-iR)t ] 
b ^ 

(where e denotes a unit vector in the x direction). 

Hence 

-?r -?r 

(Ve a ) • (Ve b ) = (cVr a r b )e^ R| [x 2 + y 2 + ( z 2 - Jl*)] (27) 

Using (6), (8), and (9) we obtain 

<M^ ■ Vx r ) = (\)(^{£)(\)e-^(t\ri 2 - 2)d|di3d^. 

R 
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Integrating over cp and substituting into (26) leads to 



3^5 



t b = i R ? / die" 51 * / dr, (? 2 + r, 2 - 2) 
1 -1 



OO 



-*nY/ (| 2 -|)e-^ d ? 

1 

= KV CR [l+a*-i(£R) 2 ] . (28) 

2 -A. 4 Summary 

Collecting together the quantities of the previous sections 
we have 

X=^7T e" ?r (1) 

t M = (x^ltlx^) =i? 2 (2) 

a 
S = <X f I V = e" ?R [ 1 + ?R + i(?R) 2 ] (14) 

t ir = <X £ |t|x r > = K 2 e" ?R [l + CR - K?R) 2 ] (28) 



2 ,-?R- - 2 



V £r = (XjjK-p -p)|X r >=-2Ce-^(l + ?R) (19) 

a d 



vg=-i[l-(l + ?R)e- 2?R ] 

2 -A. 5 Qualitative Examination of Diatomic Quantities 

The amplitude of )L. evaluated at the left nucleus is 

X r (R)=\R% e" CR 
while the amplitude of X« at the left nucleus is 

3 



(25) 



Xj (0) = \| f/w 
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Thus 

XrtR^e^x^O) (29) 



If X c were highly concentrated about the left nucleus the overlap would 



be given by 



S=e" ?R (30) 



Comparing with the correct formula (14) we see that the approximate 
form (30) has the correct exponential behavior on R but the numerical 
coefficient in (30) is correct only for R = 0. Using R = 2. 0a and 
? = 1. in (14) leads to a coefficient of 3. 33 and using R = 6. 0a leads 
to a coefficient of 19. 0, many times the value obtained with (30). 
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APPENDIX 2-B. THE LEGENDRE EXPANSION 

Consider a system such as in Fig. 34. It is often necessary to 
convert expressions involving say the distance of the electron from 
nucleus b over to a new expression involving the distance of the electron 
from nucleus a, as indicated in Fig. 35. 




a r b 
Figure 2- 35. 
The relation between these coordinates is 

r b = r a + r2 " 2r a Rcos e (1) 

A case of particular importance is to convert l/r fe over to the new 
coordinates. This leads to 



(2) 



r fe r V~l + p - 2p cos 6 

where p = R/r . If p < 1 the radical in (2) can be expanded as 
a 

, ,/ = f P * P, (cos 6) , (3) 

V 1 + p - 2p cos 0_q 

where the P* (cos 9) are the Legendre polynomials 

P (cos e) = 1 

P x (cos 0) = sin 6 (4) 

P 2 (cos 9) = \ (3 cos 2 0-1) 
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Thus (2) becomes 



(*/)' 



1 — P £ (cos0), 
r b 4=0 ( r> )* +1 



(5) 



where r, and i\ denote the lesser and greater, respectively, of r and R. 

If there is a spherically symmetric charge distribution p a (r J 

* a 

centered at a, then the total electrostatic interaction with a charge centered 
at b is 



V = /drp a (r ) -L 



b 



oo oo 



-2,1 f r" dr a p(rj 



1=0 



a*a a 






IT 



f sin Ode P. (cos 0) , (6) 



where we integrated over the cp coordinate. The Legendre polynomials 
have the property that 

o 
so that (6) becomes 



J sin 6 d0 P^ (cos 9) = 2 5 £Q , 



(7) 



V = 4tt f°°r 2 dr a p(r) i- 



= 4 * { 4" f r > r , **<*«> + L ** d ' a P»fr a )> • < 8 > 



R 



a a a a 



R 



a a a a 



The quantity 



R 



Q = 4tt / r drp(r) 



a a'a a 



(9) 



is the amount of charge within the sphere centered at a and passing 
through b. The contribution of this charge to V is just the same as if 
all this charge were concentrated at a. 
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The quantity 



6Q = 4 7rr p (r ) 6r Q (10) 



a ^a x a y a 



is the amount of charge on a sphere of radius r and thickness 6r . The 

a a 

potential within such a uniformly charged sphere is constant and equal to 
r- 5Qas implied by the second term of (8), 



r a 
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Appendix 2-C Coulomb and Exchange Integrals 
In §2.4 we indicated that 



J.. ^ K-- ^ 

i] i] 



These relationships are derived below 



(a) J i3 - > 0. 



The Coulomb integral is 

J ij = f^i *i(l)*i(l) /dr 2 ^ *j* (2) 0j(2) 
0*(1)0 . (1)^(2)0 (2) 



r !2 



Since the integrand is positive for all values of r x and r 2 , the integral 
must be positive. This integral is also denoted as 



J.. = [(f>* <b.\<b* d>.] . 



where the orbitals on the left are for electron 1 and those on the right 
are for electron 2. 

(b) K.. > 0. 



The exchange integral is 

K ij = K^)V"/ d 4 *j*<2)*i(2) 



* J. I J* 



= [or^uf^i 



To prove that K-- s* we set 

■•■J 
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*i = «i + iVi 



*i = { i + iv i 



where |. , 77- , |. , and 77- are real. This substitution leads to 
'i' '1 7 j '] 

K y = [{*i(l) - i77 i (l)}{4 j (l) + ii»j<l)} I tift) + U? i (2)}{« ] (2) - i7j.(2)}] 

We now define the charge distributions p x and p 2 as 

Pi = Ijlj + T?i^j 
P2 = «ii?j " »?i«j • 



This leads to 



Ky = [pi + ip 2 Ipi - ip 2 l 

= I Pi IpJ + i{[p 2 IPil ' [Pilp 2 ]} +[p 2 |p 2 ] 
> o =0 ^0 



and hence 



K- > . 
(c) J.. ^ K.- . 

To show that J.. ^ K-. we consider the wavefunction 

*(1,2) = 0^ - 0.0. (1) 

The electron-electron interaction energy for this state is 



<*|-L|*> = J7d Tl dT 2 -A- [<//(!, 2)| 2 ^ 0, 



i 2 " - ' * 2 r i 2 
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since the integrand is positive everywhere. Substituting the wave function 
(1) leads to 

<*iitl*> = 2 ( J irV (2) 

and hence 

J.. ^ K-. . 

Expression (2) provides the physical significance of K-. It is the change 

in electron repulsion energy upon superimposing both products of orthogonal 

orbitals A.0.. That is, the wavefunctions -— [<f)i. ± </>•&.] lead to Coulomb 

1 1 -J2 * 1 1 v 

repulsion energies of 

J ij ± Ky • 
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Appendix 2-D» Two-Electron Integrals for H 2 
2-D. 1 The Coulomb integral, J- 

The Coulomb integral, J - has the form 

J ir = (X j ^(l)x j? (l)l^ r lx r (2)x r (2)> 



(1) 



J 



Xp(2) X r (2) 



where \. and x r are * s atomic orbitals centered on the left and 
right protons. Letting 

J,(2) = J"dT lX ;(l)x»(l)^- , 



we have 



J Cr = / dT 2<M 2 ) X*(2) X r (2) • 



r x / /x r 



(2) 



(3) 



In the remainder of this section we will assume x« and X r are real. 

First we evaluate J. by expanding — — as 

* r 12 

00 ^ (k-lmhtr || II 

L. = S T . — ' < pj m l (cos X ) pi 111 1 (cos ^^e 1111 ^ 1 " ^) 



12 k=0m=-k (k+ m|)» r 



k+l k 
> 



(4) 



(the Laplace expansion) where r 12 is the distance between the points 
with spherical coordinates (r 1? 9 19 (p ± ) and (r 2 , 6 2 ,<P 2 ) , 




Note that both r x and r 2 are with respect to the left-hand Genter, 
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Since x« is spherically symmetric (about center H) the integrals over 
X and 4> x will be nonzero only when k = m = 



J C = ** f x ll dT *l*\( T nl)4 







(5) 



Breaking the interval of integration to remove the r^ yields, 



r-r 



J„(2) = 



a2 



4 ff 



a2 



/ A ( r al) r al dr al + / Xf< r al> r al r a2 dr al 
r a2 



(6) 



Letting 



X.(l) =^/-^- e" ?ral 



we find 



J l( r a2) = 



'a2 



" 2 ?r a2 
1 - e az (1 + ?r a2 ) 



(7) 



(8) 



Now we must change to elliptic coordinates in order to evaluate the 
integral over the coordinates of electron 2. 

Using elliptical coordinates £ and t\ (see Appendix 2-A) 



| = (r a + r b )/R 



1 = < r a - r b> /R 

d r = -^(£ 2 - n 2 )d^drid<p 
8 




in (8) and (3), we find 
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2 2ir i °° 
o -i 1 



(1 + 7?) 

1 °° 



i.e-? R « + n)(i + i R(l + rj) ) 



3j-»2 

"2 



= "St" / d7] /dSd-iy) 



or 



- 1 i 



1 oo 



e-SR(«-l).e- 2 ? R *(l + 4-R(i + ,)) 



0) 



j. = JW Jd./dKI-^e-^^-"), 



,3^2 

"2 



-l i 

1 OO 



+ -£#- /dij Jd«({-ij)e- 2?H * + 



(10) 



4-r,3 



-1 1 

1 oo 



+ 



£?-/d,/de««-,-)e-«H{ . 



- 1 1 



Evaluating these elementary integrals gives the following result 



j 1p = ? 



1 -2£R / 1 11 3 «„ 1 i*2„ 2 \ 

Tr - e (tit + t + t ?R + t ? R ) 



(ii) 



Problem: Evaluate ( Xi (l) x.fl) 1^- I x„(2) X„(2)> 



£\ / AjfV / I j.^ I A.j.\-/ A£< 
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2-D. 2 The Exchange Integral K 



kr 



Now we evaluate the exchange integral 

K lr = <Xj|(l) X r (2)| — | x r (l) X^ (2) >=[x^X r |x^X r ] 



(12) 



Again we can define a quantity 
1(2) = /d Tl [ Xj i(l)x r (l)-p-] 



(13) 



involving integration over the first electron. Unfortunately the 

Laplace expansion of — - will now lead to an infinite sum because 

12 

X^(l) X r (!) is n °t spherically symmetric. 

Instead we expand — - in terms of elliptical coordinates 

r i2 

using the Neuman expansion, 



12 



= A £ S (-l) m (2k + 1) 
R k=0 m=-k 



f- -» 2 

(k - | m | ) ! 



(k+ |m|)! 



x 



(14) 



x 



Pi m| [ ^] Qi m| [ ^1 *J m| toi) 4 m| (i?,) e ta <*i-*.> 



< 



5> J "k 



|m 
where P k are the associated Legendre functions and 

|m| 
Q k are the associated Legendre functions of the second kind. 

Our function X^ X r is independent of cp so the only nonzero term 
in the m summation is m = 0. 

To simplify the k summation we use the property 



/ P k (i?)P k /fa)d* = if k*k' 



-1 



k 



together with the facts that 
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P (77) = constant 
and 

p 2 (u) = l(u a -i). 

Thus since the volume element for integration is -~- (£ 2 - rj^d^ drj^d^, 
we find that by integration with respect to 77 , the only nonzero terms 
are for k = and k = 2. 

We will now skip pages of tedious algebra to the result 
(see, for example, J. C. Slater, Quantum Theory of Molecules and 
Solids , Vol. I, pp 266-273 for the details). 

It is first convenient to define 

S = e"^ R (l + £R + i£ 2 R 2 ) 

a = e + ^ R (l - £R + -H 2 R 2 ) 

1 1 t °° -t 

C= / ~ e dt - /-j-dt = 0.57722 (Euler's Constant) 



1 



GO _f 

e L 



Ei(-x) = - f - e — dt (integral logarithm) . 



With these definitions the result is 



K ir = ro?{-e" 2?R (-^- + -T-SB + 3? 2 B 2 + f? 3 B 3 ) 



+ 



or for small R 



+ -=~ [S 2 (c + in(?R)) + a*Ei(-4 £R) - 2S cr Ei (-2 £ R)]] 
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K 4r = 2^[|- - l^ 2 R 2 + (-^-+"4- in4)^ 4 R 4 ] + (higher order terms), 
Note also for small Rthat if we include terms through R 2 , 



hr =*5(|--J-S a R 2 + • • • ) 



&g r -2H7-- 2S « + • • • ; 



•• K j^ r < J i r ■ 
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APPENDIX 2-E. DETAILED ANALYSIS OF THE EXCHANGE TERMS 

In this section we will provide more detailed analysis of the 
classical and exchange terms for H* discussed in §2.2 and §2.3. 

2-E.l The Potential Energy Terms 

Just as with the total energy (see §2.2.1) the potential 
energy can be partitioned into classical and exchange terms as 
follows 



2-E-2 



where 



v„* + v +2V. 1 

V = M rr ^ v lr + _1_ 

g 2(1 + S) R 



1+S + R 



= V C * + V* , (1) 



v ir = <4|v|r>. 



v cf = v M + i W 



= i(v M + v rr ) +1 

= /drv(r)p c *(r) + i 



V s = V;£r " SVfJg 

S 1 + S 



T v 



1 + s 



(p is given in 2.3-3). Similarly 



V = V + V x 
u u 



(31 



t v * v fr - S v w (4) 
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where 

V u = T^T • (8) 

First we examine the form of V . Substituting (2. 1-4) leads to 



■,rCl ,, I 1 1 I v. 1 



The first term is just the potential energy of an isolated hydrogen atom. 
The second term is the net Coulomb interaction between a proton 
(on the right) and a hydrogen atom (on the left). As shown in App. A 
[eqn. (25)] 

^V c ^<x £ i^--^-|x 1 >= (l + i)e" 2R 



and hence 



V ci = -1 + (l + i)e" 2R (7) 



Consequently in the simple classical description (superposition of atomic 
densities), there is no bonding of H^" . [This classical description is 
equivalent to bringing up a proton to a hydrogen atom without allowing any 
changes in the wavefunction of the hydrogen atom. ] The other terms in the 
potential energy arise from interference effects. That is, they occur 
because we superimpose amplitudes rather than densities. 
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The total electron density for the g state is 

x f> x l + x r x r + 2 x i x r 

Pp. = — » 

g 2 (1 + S) 

which can be partitioned into classical and exchange parts as 

Pg = P + P g > (8) 

c 
where p is given in (2 .3-3) and 



ci- 

>P 

(1 +S) 



p x = f X l X r " S P ] ^ (9) 



Since 



and 



x 



Jdrp = 1 

Jdrp ci = l 
the integral of p A must be zero 

/drp x = . (10) 

That is, p x merely shifts density around with no net contribution to the 
total electron charge. As a result we can determine the sign of 

V x = /drv(r)p x (r) (11) 

from Fig. 2-8. Here we see that p X > near the bond midpoint 
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while p^ < near the nuclei; that is, p* leads to a shift of charge from 
the nuclear region to the bond region. Since v(r ) is much more 

negative near the nuclei than near the bond midpoint, this shift of charge 
into the bond region leads to a positive value for V (R) as shown in 
Fig. 36, 



EXER (USE : 



(a) Evaluate p* and p c * at the bond midpoint for R = 2.5 a 
[Answer: „*(») = MM , p ci ( R } = 0J86] _ 

& & IT * IT 

(b) Evaluate p* and p c ^ at the left nucleus. 

5 

[Answer: p* (0) = ^i° and p ci (0) = ^ ] . 
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2-E. 2 The Kinetic Energy Terms 

The kinetic energy of the <j> a state can be written as 

T = tjeA + trr + 2tjgr 
s 2(1 +S) 



t tJL + t Jgr 
1 +S 



where 



We will write 



where 



= t 



rr ? 



t lr = <x £ |t|x r > 



T_ rp«^J- rpX 

g g ' 



1 +s 



r t 



1 + s 
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(12) 



(13) 



T =i(t a +t rr )=t u ( i4) 

T * = t Ar - S T c * 

g 



(15) 



T t s t fr -ST cf (i e) 
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Similarly T u = T c ^ + T^ , where 

c l 
T x _ _ hr - S T _ . _^t_ (17) 

u 1 -S 1 -S 

Since T is just the atomic value of the kinetic energy (independent oLR), 
the changes in T responsible for bonding must all be contained in T*. 
Thus the plots of T and T in Fig. 12 are actually plots of T* and T x . 

The large negative value of T responsible for the bond in H* 
results from the large negative value of r t . We will now examine why 
r t is large and negative. From Appendix 1-A we know that 



- i / drfai,. v\ ; (18) 



Al- 



and 



t u = i <|vx 1 r> , 



t - - < Vv > 



Thus r t = tj, r - | S (t w + t rr ) (19) 



= A / d $ r {V XjL • Vx, - | [(Vx £ ) 2 + (V)^) 2 ] } (20) 



In order to understand the significance of the terms in (20), we will 
consider first the case where t, is modified by replacing 
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^ * VXj. in (20) by | V\ i | |v>^ |. This leads to 

an integrand of the- form 

|VX r ||Vx £ |-|[(Vx £ ) 2 + (y) (r ) 2 ] . (21) 



-r 



a 



Since X^ = e we obtain 



V ** = " X £ e r 



(where e is a unit vector in the r direction) and hence (21) becomes 

[ *r x £ ~Z (x £ + x r } ] 



However the term in brackets is just proportional to p in (9) and hence 

from (10) the resulting integral is zero. Thus it is the difference between 

the dot product term 

VX . vX 
JL r 

in (20) and the absolute value term 

in (21) that is responsible for the large negative value of r t and hence of the 
chemical bond. To emphasize this we define a function called the 
contragradience 

C(r) = |VX £ | |vx r | - VX^ • VX r (2.3-7) 

such that 

r f = -i/d 3 r C(r). (2.3-6). 

Large contragradiences lead to a large negative t. and hence strong bonds. 
As discussed in §2.3.3 and illustrated in Fig. 14, the largest values of 
C(r) occur for points in between the nuclei. 
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2-E.3 Specific Results for H* 

From App. A the explicit form of r, for H* is 

T t = -|R e K 

Thus as expected, t,— * as R— >0 and as R-+<». The minimum in t occurs 
for R = 2. Thus we would expect the maximum bonding effect to occur 
near R = 2a . Indeed this is the optimum R of the exact wavefunction of 
H^. Approximating the bond strength as 

rpX _ I 

r g ~ "TTF 
we obtain for R = 2 (where S =0.6) 

T x = -^ = -0.12h 

g (3)(1.6) 

= 3. leV 

In fact the bond energy is only O.lh = 2.5 eV, 
and the T* term does dominate the bonding. For 
quantitative considerations we should, of course, use the total E . 
From Appendix A the explicit form of r y for H+ is 

\ = e " R [|" f R+ |R 2 ] > ( 22 ) 

-2T? 
neglecting terms of order e . Combining with T t leads to 

T = e ~ R[ R"f R] (23) 

(neglecting terms of order e ^), whereas 



S = e~ R [l+R+ !r 2 1 
L 3 J 



Thus 



S = -R- (1 "R" + W ) + ° ( R^ ) 
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and hence r ~ - 



R 



(24) 



r + 



Various energies for Hj are tabulated in Table I, 



Table 2 


-1. Energy auantities for the LCAO wavefanetions of the g state of H 2 . All 
quantities are in atomic units. 


R 


s , 


AV ci 


V X 


AV total 


T X - AT total 


AE total 


00 




0.0 a 


0.0 


0.0 a 


0.0 a 


0.0 a 


10.0 


0.00201 


0. 00000 


0. 00121 


0. 00121 


-0.00151 


-0.00030 


8.0 


0.01018 


0. 00000 


0. 00536 


0. 00536 


-0.00708 


-0.00173 


6.0 


0.04T10 


0. 00001 


0.01933 


0.01934 


-0.02841 


-0.00907 


5.0 


0.09658 


0.00005 


0. 03195 


0.03200 


-0.05120 


-0.01920 


4.0 


0.18926 


0.00042 


0. 04485 


0.04527 


-0.08214 


-0.03687 


3.5 


0.25919 


0. 001 17 


0.04858 


0.04S75 


-0.09792 


-0.04817 


3.0 


0.34851 


0.00331 


0. 04837 


0.05168 


-0.11076 


-0.05908 


2.5 


0.45831 


0. 00943 


0. 04301 


0. 05244 


-0.11727 


-0. 06483 


2.0 


0.58645 


0.02747 


0.03250 


0.05997 


-0.11374 


-0.05377 


1.5 


0.72517 


0.08298 


0.01901 


0.10199 


-0. 09700 


+0. 00499 


1.0 


0.85839 


0. 27067 


0.00695 


0.27762 


-0.06599 


0.21163 


0.5 


o 96034 


1.10364 


0. 00080 


1.10443 


-0.02578 


1.07865 



a The values at R = « are V ci = -1.0, V total = -1.0, T total = 0.5, E total = 0.5, 
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Summary 

The energy € of any approximate wavefunction O is an upper bound on 
the exact energy of the ground state E , 

e > E , (3.1-3) 

leading to the variational condition: If an approximate wavefunction (and hence 
the energy) is a function of some parameter > , then the optimum wavefunction 
satisfies the (necessary) condition 

— =0. (3.1-4) 

Expanding the unknown wavefunction cp in terms of a basis 

P 

and applying the variational condition leads to a set of matrix equations, 

HC =ESC (3.1-25) 

for obtaining the optimum coefficients (i.e., wavefunction). 

More exact wavefunctions of H* are considered in §3.2, but in § 3. 3 we 
find that the description of bonding in terms of exchange energies is retained. 



In §3.5 we present an overview of three useful methods for wavefunctions: 
(a) The Hartree Fock (HF) method is a generalization of the MO wavefunction 
in which the wavefunction (ground state of a two-electron system) is taken as 

$ HF (1,2) = 0(1)0(2) (3.5-10) 
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and the orbitals cj> optimized by solving the differential equation 



H HF = (h + J , ) = € (3. 5-20) 



or the matrix equation 

S C = eSC . (3.5-18b) 

These equations are nonlinear and must be solved iteratively. 

(b) The generalized valence bond (GVB) method is a generalization of the VB 

method with the wavefunction taken of the form 

4> GVB (1,2) = a (l) b (2) + 4> b (l) a (2) (3.5-22) 

and the orbitals a and <£> b optimized. This leads to two matrix equations 



I a C a = c a SC a 



I C b = € b SC b 



(3.5-31) 



analogous to the HF equations (3.5-18) and to two differential equations 
analogous to (3.5-20). 

(c) The configuration interaction (CI) method with the wavefunction taken of 
the form 

<f> CI (l,2) = D C X (1)X (2) 

For the ground state this wavefunction can always be written in terms of 
naturaj^ojjitals {x } as 



4> CI (1,2) = E C M x^(l)x^(2) • (3.5-53) 



In § 3. 5. 5 we find that the ground state of any two electron system is 
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nodeless and symmetric 

$(1,2) = $(2,1). 

In §3.6 and §3.7 we find that the HF wavefunction accounts for all but 
about 1.1 eV of the energy for He and Hj and that a CI wavefunction with five 
NO T s accounts for all but about 0. 15 eV. The four correlating NO's for this 
wavefunction all involve one nodal plane. 

All these methods involve expansions in terms of basis sets. For He it 
is possible to obtain highly accurate HF and GVB wavefunctions with only two 
(s-like) basis functions [the double valence (DV) basis ] and for H 2 similar 
quality wavefunctions can be obtained with six basis functions (two s and one p 
on each center) , [the DVP or double valence plus polarization basis ] . 
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Chapter 3. More Exact Wavefunctions for One- and Two-Electron Systems 
In Chapter 2 we discussed approximate wavefunctions for H^ 

and for H 2 

^J? ^r ^ ^r ^1 * 

In this chapter we will consider some aspects of more exact wavefunctions 
for these molecules and for the two-electron atom, He. Our emphasis will 
be on qualitative ideas, explicit methods for solving for the wavefunctions 
will not be developed until Chi 20c. 

In §3.1 we examine the variational principle and in §3.5 we use this 
principle to examine some of the useful methods for calculating wavefunctions: 
Hartree-Fock (HF), generalized valence bond (GVB), and configuration 
interaction (CI) . 

Accurate wavefunctions for H*, He, and H 2 are discussed in sections 
§3.2, §3.6 and §3.7, respectively. In §3.3 we reexamine the nature of the 
bond in H* and in § 3. 7. 2 we reexamine the bond in H 2 . 
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§3.1 The Variational Principle 

The wavefunction for the ground state of the hydrogen atom has the form 
e~" r and the wavefunction for the ground state of the harmonic oscillator is 

However, there are very few interesting systems for which the Schrodinger 
equation can be solved exactly. Even so there are general procedures allowing 
one to obtain highly accurate (indeed arbitrarily accurate) wavefunctions. 
The powerful tool needed is the variational principle which will be discussed 
extensively in Chi 20c. Here we will outline some of the key results of this 
principle in order to clarify the qualitative discussions of Chi 20a. These key 
results are: 
a. Upper bound theorem. 

If \f/ and E are the exact ground state wavefunction and energy of a 
system 

3C^ = E o^o (1) 

and if <p is an approximate wavefunction with energy 

<*olrc|* > 
€ o = — i r > ( 2 ) 



then 



£o ^ E o . (3) 



That is, the energy evaluated using any approximate wavefunction is an upper 
bound in the exact energy of the ground state. 
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b. Variational condition 

Given an approximate wavefunction (r) depending upon some 

A 

parameter \ , 

e.g., <j> x = e 

then the optimum value of x for describing the system must satisfy the 
necessary condition 

^ = 0, (4) 



dX 



where 






Equation (4) is referred to as the variational condition . 

c. Basis set expansions 

Given a set of functions 



{Xjto, M =1,2,...P} , (5) 

the optimum wavefunction of the form 

P 

= 3 °» V (6) 

satisfies the condition 

Sh C =eSs C, (7a) 

or in matrix notation 

H C = ESC, (7b) 
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where 



H = (x 3C x ) 



S = (y Y ) 



(8) 



The set of functions (5) used for expansion of the unknown function cf) is called 
a basis set . The unknowns are the coefficients 

{C^; n =1,...P} (9) 

which are obtained by solving the matrix equations (7). 

§3.1.1 Discussion Of Upper Bound Theorem 

The upper bound theorem (3) is easy to derive. Consider that the exact 
eigenstates {^} of the Hamiltonian were known 



3C 



^i = E i^i (10) 



with i = as the ground state. 

Since the set of functions {;//.} is complete, we can expand any 
approximate wavefunction <£ as 

*. = S q^ . (ii) 

Of course, in a real problem we will not know the functions {\p} , and hence 

we will not be able to calculate the {C,}. However, the analysis in this section 

will serve to establish a relation between the approximate and exact solutions. 
From (10) and (11) 

K(t> = S C-3C^ = V C,E.^, (12) 

v i ii i ii'i 

and the energy of the approximate wavefunction $ is 
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€o ee {0 o |3c|0 (> >=SEc j *C i E i (^.| l // i > = Slql'E.. (13) 

Assuming <t> is normalized, 

1 = (<t> \<P ) = BcJ'C^jl^) = plcj 2 . (14) 

Using (14) in (13) we obtain 

e - E = S|C.| 2 (E.-E ). (15) 

i 

Since E i ^ E and | cj 2 > 0, the right-hand side of (15) is necessarily 
positive, and hence 



or 



€ ^ E 
o ^o 



(16) 



That is, the energy calculated for any wavefunction is never lower than 
the exact energy for the ground state wavefunction * Thus we say that 
e n is an upper bound on E ft . 

* 
In deriving (16) we assumed that (f> could be expanded in terms of 

the eigenfunctions of K . This requires that <£ satisfy the boundary 
conditions for 3C. For example, if the boundary 

conditions for the system described by (10) were such that all wave- 
functions are anti- symmetric, then we could not allow (j> to contain a 
symmetric part. 
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§3. 1. 2 Discussion of the Variational Principle 

Since any approximate wavefunction must yeild an energy above the energy 
of the exact (ground state) wavefunction, we have a useful criterion for improving 
approximate wave functions. Namely, if you can find some change in the wave- 
function that leads to a lower energy, then do it! And reject any changes that 
increase the energy. Ultimately if we consider all possible changes in the 
wavefunction, this procedure must yield the exact wavefunction. Normally we 
lose patience before considering all changes in the wavefunction, and instead we 
consider functions of certain restricted classes. Here our basic criterion for 
approximating the wavefunction will be to select that function of our specific 
restricted class leading to the lowest energy. If X is some variable parameter 
for the restricted set of functions being considered, then the optimum wave- 
function must satisfy 

(17) 



dE 







ax 



since otherwise a lower energy could be obtained by changing X a bit, as 
illustrated in Fig. 1. 




^X 



optt 



Figure 3-1, 
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This criterion for optimizing a wavefunction is called the variational 
principle and forms the basis of all methods we will consider for 
determining wavefunctions . It should be noted here that (17) is not 
sufficient to guarantee a minimum with respect to variations in \ 
(this requires — — - > ) and even if a minimum is found it need not 
in general be the minimum leading to the lowest energy. Fortunately 
for the types of problems we deal with, these potential difficulties can 
usually be avoided. 
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§ 3 . 1 . 2a Parameter Optimization 

Consider as an approximation to the ground state of the hydrogen 
atom the function 

<M«) = *~ ar2 (18) 

where a is a parameter. To determine the value of a minimizing 
the energy, we first calculate the energy as a function of a, 

e(oo = <«oi-iv a - 1 k>/<«> k> 



r 



3a 



-fe-) 1 • 



(19) 



The optimum value of a is given by 



dE < a ) _ ft _ 2 i 



do? 
or 



f— 1* 



a o P t = "ST = °- 283 --' 

Substituting this into (19) we obtain 

E Kpt) = "^ ="0.4244... 

(recalling that the exact energy is E = -0. 5). Thus, even though (18) is con- 
siderably different from the exact eigenfunction for the ground state of the 
hydrogen atom, by optimizing a we are able to account for 84. 9% of the 
energy. 

Exercise, (a) Derive the energy expression (19). 

(b) Consider 

-CUT 

ze " 
as an approximation to the 2p orbital of hydrogen atom. 

Find the optimum a . 

(c) How would you proceed to obtain a corresponding 
approximation to the 2s orbital of hydrogen atom ? 
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§3.1.3 Basis Set Expansions 

We will now use the variational principle to determine the best 
representation of an approximate wavefunction as an expansion (6) in terms 
of the functions of some finite basis set (5). The energy is 



E = 



(0 |gcl») 



_N 
D 



(20) 



where 



N - £ C*H C 

D = S C*S C 
li, v ii \iv v 



(21) 
(22) 



and H Ml , S (1J are given in (8). (We do not assume here that the basis functions 
are orthonormal; they must of course be linearly independent.) 

The energy (20) depends on the P parameters {C }, and thus from the 
variational principle we require that 



8E 

ac 



= 0, for ii=l, 2, ...,P . 



From (20) this leads to 



3E 


1 


3N 


d % 


D 


3C M 




1 


" 3N 




D 


_ ac 



N 3D 



ir ac, 



u 



- E 



M 



ap 

ac 



= 



and hence 



8N 

ac 



- E 



JD 

ac 



= 



/x 



(23) 
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Assuming that the basis functions {\ } and coefficients {C } are all real, 
we obtain 



and hence 



Tj (H - e.S ) C 



LV 



= 



(24) 



Sh c = eSs c 



(7) 



In matrix notation (7) becomes 



HC = ESC. 



(25) 



If the basis functions are orthonormal 



s = a 

liv \iv ' 



(26) 



the variational condition (25) becomes 



HC 



EC. 



(27) 



Thus the variational principle leads to a finite matrix equation directly analogous 
to the Schrodinger equation. Indeed if a complete set of basis functions is used, 
the solution of (17) or (27) is the exact solution of the Schrodinger equation! 
Although the wavefunction and basis functions were written as one electron 
functions, this procedure applies identically for many-electron wavefunctions . 



The more general case leads to the same equations. 

Note that if the basis functions are real then H = H and S = S 



# 
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§3.2 Accurate Wavefunctions for H^ 

The LCAO wavefunction of H^~ discussed in § 2. 2 is an approximate 
wavefunction and does not provide a quantitatively accurate description of 

H? near R . In this section we will discuss more accurate wavefunctions of 

c e 

H^. First we consider a useful intermediate level description, the MBS 
wavefunction. 



§3.2.1 Scaled LCAO Wavefunctions 

We will describe the wavefunction of H^" in terms of linear combinations 
of two orbitals, x* and x r , centered on each proton, but rather than atomic 
orbitals we will use scaled atomic -like orbitals 



X4 =V(?/7r)e^ r a (1) 

X r =7(pA)e^ r b 

The scaling parameter £ is referred to as an orbital exponent. Use of 5 = 1 
leads back to the LCAO description of §2. 2, £> 1 leads to more contracted 
orbitals, while ?< 1 leads to more diffuse orbitals. 

Using the basis set (1) the wavefunctions of H 2 + have the form 

*g = ( V \)/V2(l +S) 

, (2) 

^u = ( "\ + V /V2 ^ " S) 

just as in § 2. 2. However the energies of these wavefunctions depend upon both 

£ and R (see App.-2A for the specific dependence of the integrals on £). 

At each R we will now use the £ leading to the lowest energy. Since the 

forms of E and E are different, the optimum ? will be different for the 

g and u states, as shown in Fig.l. 
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Figure 3-1. The optimum orbital exponents for the MBS descriptions 



+ 



of the g and u states of H, . 
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As shown in Fig. 2 and Table 1, the improvement in the energy for the g state 
is quite remarkable, leading to energies close to the exact answer. For the 
u state both the LCAO and MBS energies are quite close to the exact answer. 
In discussing such wavefunctions we will use the following terminology 

(a) LCAO denotes the use of a linear combination of atomic orbitals using 
the orbital exponents of the atoms. 

(b) MBS ( minimal basis set ) indicates the smallest set of atomic-like functions 
that would describe the case of R = oo . For finite R the orbital exponents will 
generally be optimized. The results of MBS calculations will be discussed 
further after a discussion of the exact wavefunctions of H*. 



Table 3-1 
Optimum Bond Length (R e ) and Bond Strength (D e ) for the g state of H 2 + . 
All quantities in atomic units. 











Relativistic c 




Non-Relativistic 
Neglect Nuclear Kinetic Energy 


Neglect 
Nuclear 
Kinetic 
Energy 


Include 
Nuclear 
Kinetic ■, 
Energy d 




LCAO a 


MBS a 


EXACT 


EXACT i// 


EXACT i// 


R e 


2.493 


2.00 


2.00379 c 


2.00376 c 


2.00562 c 


D e 


0.06500 


0.08651 


0.102635 C 


0.10264 C 


0.101785 C 



a. Unpublished calculations, Wadt, Olafson and Goddard. 

b. D. R. Bates, K. Ledsham, and A. L. Stewart, Phil. Trans. Roy. Soc. 
A246^215 (1953). 

c. S. K. Luke, G. Hunter, R. P. McFachran, and M. Cohen, JCP 50 
1644 (1969). 

d. G. Hunter and H. O. Pritchard, JCP 46 2153 (1967). 
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Figure 3-2 The LCAO (? = 1.0), MBS (optimum ?), and exact energies for 
the g and u states of H 2 + . Note that Fig. (b) does contain three 
different lines. The vertical scale of (b) is twice that of (a). 
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§3-2.2 The Exact Wavefunction for H 9 + 

Previously we considered approximate solutions 
of the SchrBdinger equation 

3C;// R (r) = E(R)^ R (r), (3) 

where the Hamiltonian is 

x = 4v 2 - ^ - J- + * (4) 

a d 

(see Fig. 2-1 for the coordinate system). Exact solutions to (3) have also 
been obtained as will now be described. 

From §3.1 we can obtain arbitrarily accurate wavefunctions for H^* 
by expanding the orbital in terms of a sufficiently general basis 

Ix^; /i=i» 2 , . . .p}, 

p 

<Mr) = E C M x^r) , 

where the expansion coefficients are obtained by solving the P by P 
matrix equation 

HC =CE , 
with 



V = <X JLt |3C|X 1/ > 



(assuming the basis to be orthonormal). As the basis set is made more 
complete (P — » °°) the wavefunction approaches the exact wavefunction. 

Although the above procedure is practical, it is possible for K* to 
solve directly for the exact solutions. The procedure is examined in 
Appendix D. 
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§3.2.3 Comparison of Wavefunetions and Energies 

The various wavefunetions of the g and u states are compared in 
Fig. 3 f or R = 2a . For the g state we see that the shape of the LCAO 
wavefunction in the bond region is in good agreement with the exact wave- 
function; however, the magnitude of the density in the bond region is ~ 25 to 
30% low. The MBS description leads to reasonably good densities near the 
nuclei but too low a density in the bond region. Thus with MBS the shape 
of the wavefunction is not well described. 

In the u state the LCAO wavefunction is in much better agreement with 
the exact wavefunction than is the MBS wavefunction. I do not understand the 
reason for this. In addition, I do not understand why the energies of these 
wavefunetions are so similar (see Fig. 2b) despite the large differences in 
the wavefunetions. 

In Fig. 4 we compare the LCAO and MBS wavefunetions as a function 
of R, finding that the LCAO description does reasonably well for R>4a . 
Note the large difference in the behavior of the g and u states for small R. 
These differences were also manifest in the optimum exponents of Fig. 1. 
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Figure 3-3 The wavefunctions of H 2 + at R = 2. a . 
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Figure 3-4 The amplitudes of the LCAO (dashed) and MBS (solid) wave- 
functions for (a) the g state and (b) the u state of H 2 + . 
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§3.3 More on the Chemical Bond 

In § 2. 2 and 2. 3 we analyzed the bond of H+ in terms of the LCAO 
description. Now we will reexamine the bond using more accurate wavef unctions. 

With more accurate wavefunctions we still find (§3.3.1) that the exchange 
energy E x , (more specifically the exchange kinetic energy T x part of E x ) is 
responsible for the bonding or antibonding of the g and u states of H+. 

On the other hand, partitioning the energy into the total potential 
energy (V) and the total kinetic energy (T), we find (§3.3.2) that neither can 
be solely responsible for bonding. 

§3.3.1 The Classical and Exchange Energies 

Defining classical and exchange terms just as in §2.2, 

E = E e£ +E X (1) 

E cSL = (x^SClxp (2) 

but using the MBS wavefunctions, we obtain the results of Fig. 5. Thus the 

exchange eaqrgy dominates the bonding just as for the LCAO wavefunction. 

Partitioning the E x into potential and kinetic parts, V x and T 3 ^, 

E x =y x +T x (3) 

as in Fig. 5, we see that T* favors bond formation while V* opposes it, 

just as for the LCAO wavefunction. 

Thus in terms of the classical and exchange quantities the LCAO and 

MBS descriptions are quite similar. In both cases it is the large decrease in T* 

that is responsible for bond formation. Just as discussed in § 2. 3, T* 

is large and negative because the atomic orbitals are contragradient in the 

region between the nuclei. In particular the T* is similar in character for 
the LCAO and MBS descriptions. With £>1 the gradients get larger and 
favor a smaller R so that the differences in T^ for the LCAO and MBS 
descriptions are easily understood. 



3.3-2 
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,x 



Figure 3-5. (a) The total energy (E) and the components E and E for the MBS 



r + 



X 



rX 



wavefunction of the g state of Hj . (b) The T and V components of E 
All quantities are relative to R = <» . 



3. 3. la An Ambiguity 3.3-3 

There is a flaw with this procedure of decomposing the 
energy into classical and exchange parts. Adding a second basis function 
on each center, say X - and X 2r , and optimizing the coefficients (§3. 1) 
leads to 

*g= C i( X a + X ir> + C 2( X 2 £ +X 2r) 
and adding additional functions we ultimately obtain the exact wavefunction 
in the form 

OO 

H'Si c i^*i + ^- (4) 

Thus we can define optimum left and right orbitals as 

*i S ? <i X k4 (5) 

X = 2ji c \ 
r i k A kr 

and obtain an exchange energy for the exact wavefunction. The problem is 



that for the exact wavefunction there is not a unique choice for the left 

and right functions X - and X . As a result, there is some ambiguity in 

the exchange energy for the exact wavefunction. On the other hand, with 

optimized basis functions only a few functions (say, two s and one p on each 

z 

center) lead to quite accurate descriptions but with no ambiguity in the 
decomposition (5). 

§3.3.2 Potential and Kinetic Energies 

Rather than the partition (1) of the energy into classical and exchange 
terms, it has been much more common to partition the energy into total 
potential energy, V, and total kinetic energy, T, 

E = T + V . (6) 

I believe that this partition mixes up the things characteristic of bonding with 
other quantities that are nearly independent of bonding with the result that neither 
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quantity (T or V) consistently contains the bonding stuff. A good illustration 

of this is to compare the quantities for the LCAO and MBS wave functions 

of H 2 + . As shown in § 3. 3. 1 the classical and exchange energies behave 

very similarly for these two cases. However as shown in Fig. 6 the 

behavior of T and V for the MBS and LCAO wavefunctions is markedly 

different. Thus for LCAO the T(R) is always lower than T(<*>) while 

V(R) is always higher than V(<*>) . This might suggest that it is kinetic 

energy that is responsible for the bond. However for the MBS wavefunction 

T(R) is below T(<*>) only for R >2. 7a . Thus at R o = 2ao , T(R) >T(°o) and it 

total e 

would be ludicrous to assert that the/kinetic energy is the quantity dominating 

bonding! On the other hand, in the MBS wavefunction,V(R) > V(°°) for 
R = 3. 5 a . Thus although V(R) dominates the bond at R e , it opposes bond 
formation for R ^ 3.5a . Furthermore, for LCAO, V(R) opposes bonding 
for all R. 

Such difficulties have convinced me that (6) is not a useful partition of 
the energy. The key indication of this is that although the total energy changes 
monotonically from R = « to R e , the V and T for the MBS and exact wave- 
functions are not monotonic, each dominating the energy over different regions. 

Hence neither can be uniquely responsible for bonding. 

* 
Occasionally, energy curves are analyzed by partitioning the V into 

various parts 

V = V en + V 1111 + V ee (7) 

where ee, en, and nn denote electron-electron repulsion, electron-nuclear 

attraction, and nuclear-nuclear repulsion terms (ee is not present f or £L> ). 

en 
As shown in Fig. 7a, each term is monotonic, with V 

en 
decreasing with R. One might conclude from this that it is V that is 

responsible for bond formation. However as seen from Fig. 7b the 



Usually in the analysis of rotational and conformational barriers in 
polyatomic molecules. 
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Figure 3-6 The kinetic and total potential energies for (a,b) the 
g state and (c,d) the u state of H 2 + . All quantities 
are relative to the value for R = «. 
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Figure 3-7. 



The total potential energy (V) and the partition into V 1111 
(nuclear -nuclear energy) and V en (electron-nuclear energy) 
for (a) the g state and (b) the u state of H* (exact wave functions), 
All quantities are relative to the value for R = «> . 
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V en and V nn are also monotonic for the u state and again V 611 decreases 

en nn 

with R, but this state is repulsive! Thus despite similar V and V 

for g and u, we obtain radically different potential curves. My conclusion 

en 
is that V is dominated by quantities other than those responsible for 

bond formation. 
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§3.5 Overview of Theoretical Methods 

In this course the plan is to cover the important qualitative ideas in 
Chi 20a and to examine the theoretical methods in detail in Chi 20c. However, 

familiarity with the theoretical methods is important for discussing 
qualitative ideas and hence we will outline these methods. 

§3.5.1 Basis Sets 

Several methods involve solving for the optimum shape of one-electron 
orbitals <pAr) . The general procedure for carrying out such calculations 
involves selection of a basis^set 

{x u ; m = i,...p}, (l) 

suitable for describing the optimum orbitals 

i( r)= Sc iX (r). (2) 

/i=l 

Here the basis functions are fixed and hence selection of the optimum 
coefficients 

{ C ii» C 2i> #,- » C Pi"l' (3) 

serves to determine the orbital <f>Ar). This procedure is analogous to a 
Fourier expansion where harmonic functions (sines and cosines) are used 
as basis functions in (1). 

For an exact description of the optimum orbital it is generally 
necessary to use an infinite number (a complete set) of basis functions, 
however for practical reasons we must use a finite set. Indeed from 
numerous studies of molecular wavefunctions there are principles that can 
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be used to select rather small basis sets that yield quite accurate wave- 
functions. 

In evaluating the wavefunctions and energies using a basis set as in 
(1), we must evaluate integrals of the form 

<*J h IV (4a) 

<X M (1) V2)|r7-|x a (l) V 2 )> > (4b) 

where the functions may be centered at various regions of space. Thus 
an important criterion in selecting the basis is that the molecular 
integrals be practicable to evaluate. In order to obtain the best wave- 
functions with the fewest basis functions, we want to choose the basis 
functions to have shapes characteristic of the eigenstates of the 
molecular systems. 

For a Coulomb potential (i.e. the hydrogen atom), the eigenstates 
have the form 

1— ZjY 
s: e 

2s: (r - a)e zZ,r 

2p z :r cos 6e"^ Zr (5) 

2p • r sin €> cos <p e" 2 r 

-- Zr 
2p i : r sin G sin cp e 2 

3s: (r - j3r + a)e 3 ^ r 

etc. 
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(where normalization is ignored and the constants a and /3 are unimportant 
to our considerations here). In order to describe with a finite number of 
basis functions singular characteristics such as the cusps occurring near 
the various nuclei, we should include in our basis set functions having 
similar singular characteristics. Thus for a molecular system we should 
use atomic functions like (5) centered upon the various nuclei of the molecule, 

The radial parts of the functions in (5) all can be built from functions 
of the form 

r n e"^ r , (6) 

where various values of n and of the orbital exponent , £, must be allowed. 
Functions of the form (6) are preferable to the hydrogen atom orbitals (5) 
since (6) is more convenient for evaluating the molecular integrals. 
Combining functions of the form (6) with appropriate angular functions 
(Z p , the real spherical Harmonics) leads to a convenient set of one- 
particle orbitals 



n -£r 



r e 



2^(0, <P) (7) 



for use in atomic and molecular wavefunctions . These functions (7) 

# 
are referred to as Sla^e^Ju^ictions (or Slater-type orbitals, STO) 

in honor of an early exponent* of such functions. They are denoted 

as Is, 2s, 2p, etc. just as for hydrogen atom orbitals. The orbital 

exgonent^, £, in (6) is considered as an adjustable parameter and is 

generally chosen as the optimum value for the particular molecule and 

basis set of interest [rather than taken as ? = Z/n as suggested by (5a)] . 



#We will use the term function when referring to an arbitrary function 

as in a basis function and the term orbital when referring to a specific 
optimized orbital as in a HF or GVB orbital. 

*J. C. Slater, Phys. Rev., 36, 57 (1930). 
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For example, a good basis for describing the wave function 
for H 2 is to use two Is Slater functions (denoted as Is and ls f ) ? a 2s 
Slater function, and a set of the 2p Slater functions (2p_, 2p and 2p Tr ) 

Z X y 

+ 

on each center. The optimum exponents at R = 1. 4 a are ' 



?(ls) =0.965 

£(ls f ) =1.43 

?(2s)=1.16 (8) 

£(2pz) =1.87 

£(2px) = ^(2py) =1.71 

(where the molecular axis is along z). 

With this basis the CI wavefunction leads to an energy of -1. 16696 h 
(at R = 1.4 a ), 99.4% of the exact answer* -1.17447 h. Note that the 
optimum orbital exponents are significantly different from the values 
for the free atom 

? 2 S = °- 5 

^p = - 5 - 



'A. D. McLean, A. Weiss, andM. Yoshimine, Rev. Mod. Phys., 
32, 211 (1960). 
*W. Kolos and L. Wolniewicz, JCP41, 3663 (1964). 
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The second type of basis functions commonly used in molecular 
calculations are Gaussian functions where the e"^ r of (7) is replaced 



by e~ aT and n is taken as I , 



r l e' ar2 Z lm (9 } cp) (9) 

Although Gaussian functions have the wrong behavior as r -* and as 
r -» oo , they serve just as well as Slater functions in describing the 
valence orbitals and the bonds of molecules. The major advantage of 
Gaussian functions is that the molecular integrals (4) required for large 
molecules are much simpler (and less time consuming) than for Slater 
functions. 

Generally the basis sets are optimized for the atoms. If properly 
carried out, the atomic basis sets supplemented by a few additional 
functions (polarization functions) serve to provide very accurate 
descriptions of the molecular wavefunctions. 

§3.5.2 The Hartree-Fock Method 
§3.5. 2a The Basic Equations 

In §2. 2 we described the simple MO wavefunction of H 2 in which 
the two-electron wavefunction is expressed as 

$(1,2) = 0(1)0(2) (10) 

where <p is the MO 



4> = (X, + Xj/V2(l +S) 



r 



and X| and x r are hydrogen orbitals centered on the two nuclei. Now 
we will consider the case where (p is allowed to be completely general. 



3.5-6 



Thus if { x } is some basis set, we write 



M 



= Z) c x 



(11) 



with the coefficients {C } chosen so that $ in (10) leads to the lowest 
possible energy. 

The energy of (10) is 



E = <$|3C ei |$>/<*|*> 



= 2<0[h|0>+ J^ + 



R 



(12) 



where 



J <H> = (**!**) = J" d *i* <1)>(1) /d r. 



.3 tf>*(2)0(2) 



(13) 



ia 



and 



(010) = 1. 



(14) 



Applying the variational principle to (12) 



3E 



3C. 



= 



J* 



with the constraint (14), leads to 



<v< h + V e )^> = 



(15) 
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where e is referred to as the orbital energy, 

€ = <*[h|«> + J^, (16) 

and 

r 3 ^*( r 2 )^( r 2 ) 

J ^( r i) = f d r 2 (17) 

P 12 

is the electrostatic potential at point r x due to the charge density 
|$(r 2 )| 2 integrated over all r 2 . 

Substituting (11) into (15) leads to 



S (H „ - €S „)C = 0; h = 1,2, ...P. 



(18a) 



where 



V^J^ + W 



(19) 



In matrix notation (18a) is written 

HC = SCe . (18b) 

Since 



< x J J *iV = ^/"VV^'i^"^ 



H is a function of the unknowns {C } and (18) is nonlinear. Since the 
basis functions {x } are known, all integrals in (19) can be evaluated 
just once so that (18) becomes a (nonlinear) algebraic equation. 

In order for the variational condition (15) to be satisfied for all 
X of a complete set, the function 
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(h + J^ - e) <$> 
must be zero. That is, the differential equation 



(h + J J 4> = €<f> 



(20) 



or (for H 2 ) 



(-iV 2 - ± -± + j,) ( j ) =e<t> 



r a r b # 

must be satisfied in order that be a completely optimum function. 
The resulting optimum wavefunction (10) is called the Hartree-Fock (HF) 
wavejfunction^ the optimum orbital of (20) is called the Hartrjjej^Fock^ 
orbital, and this whole approach is called theJTartr^e^J^c^method, all 
in honor of the Englishman D. R # Hartree and the Russian V. Fock 
(sometimes Fok) who first developed it. 

§ 3.5,2b Solution of the Hartree-Fock Equations 

The differential equation (20) is not linear in cj> since JL depends 
upon 0. The usual approach to solving (20) is the iterative method in which 
we guess the orbital, O , evaluate J, and solve the linear equation 

(h + J^) *! = e ^ 
for a new orbital <p L . Then <j> 1 is used to evaluate a new J , and 

(h + J^ ) 2 = € 2 

is solved for a new orbital cp 2 . This process is continued until it converges, 
that is, until * « * . 
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For atoms the HF equations (20) can be reduced to one -dimension 
and solved numerically. However, for molecules the only practical 
procedure is to use a finite basis set and to solve the resulting matrix 
equations (18). These are also solved iteratively. One guesses the 
coefficients {C°} and evaluates the H of (19). With H fixed, the 
matrix equations (18) are linear and easy to solve for a new set of 
coefficients {c' } . This process is continued until it converges. 

By solving the matrix HF equation (18) for larger and larger basis 
sets one can in the limit approach the results of solving the numerical 
equations (20). Indeed by proper choice of the basis functions it is 
possible to obtain very accurate solutions for very small P, e.g., P = 2 
for He and P = 6 for H 2 . 

§3.5. 2c Historical Note 

Before real quantum mechanics (i.e. , the work of Schrodinger, 
Heisenberg, and their contemporaries in 1925, 1926), physicists (and 
chemists) were attempting to understand the structure of atoms (and 
molecules) on the basis of a many-electron Bohr atom. The idea was 
that each electron moved along a different Bohr orbit, experiencing 
electrostatic interactions due to all the other electrons but satisfying 
various (postulated) rules in order to obtain agreement with the periodic 
properties of the elements. 
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Hartree in 1923 suggested* approximating this problem by assuming 
that the average interaction with the other electrons leads to a net potential 
that is a function only of the distance from the nucleus. He then tried to 
determine the form of V(r) [same V(r) for all orbitals] by fitting to the 
experimental energies of the orbits of various electrons (e.g. , from 
X-ray data). 

After quantum mechanics Hartree realized that he could convert 
this idea into quantum mechanics and actually solve for the potential and 
orbital. He thus started directly with (20) (and its generalization for more 
electrons) and began solving for atomic wavefunctions. These equations 
are called the Hartree equations [Hartree, a properly modest English 

gentleman, continued to call them the Self-Consistent Field (SCF) equations], 

# 
Slater in 1930 pointed out that Hartree 's equations could be derived 

using the variational principle; thus putting Hartree f s ad hoc approach 

on a more fundamental basis. Slater also pointed out that for 

many-electron atoms there are additional terms (we call them exchange 

terms) that should be in the wavefunction (from the Pauli principle, see 

Chapter 4), but Slater showed that these terms were of the same size as 

the intrinsic (correlation) errors in the Hartree approach and did not 

pursue them further. 



*Proc. Camb. Phil. Soc. 21, 625 (1923). 



tproc. Camb. Phil. Soc. 24, 89 (1927), outlined the mathematical techniques 

is 



and ibid. , p. Ill, reported Hartree wavefunctions for He (e =0. 917- h) 



and approximate calculations for Rb. 
# Phys. Rev. 35, 210 (1930). 
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V. Fock* in 1930 included the Pauli principle and derived the 
corresponding variational equations obtaining Hartree T s equations but with 
additional exchange terms. These equations are now known as the 
HartrjjejJFock equations. (Hartree referred to them as SCF with exchange . ) 

The approach of using a finite basis set for obtaining HF wave- 
functions (rather than solving numerically, a la Hartree) is sometimes 

called HF Roothaan in honor of the early leader in the development and 

f 
application of this procedure. We will make no such distinctions, although 

Roothaan T s paper is amazingly complete and Roothaan is rightfully 
credited with the development of the basis set expansion (BSE) approach. 
An early application of BSE was by C. A. Coulson [Proc. Camb. Phil. 
Soc. 34, 204 (1938)] who concluded that BSE was not practical and that 
SCF orbitals would not prove to be very useful for molecular structures. 
Basically Coulson showed that the same effort required to obtain an 
accurate HF wavefunction would, if applied to other forms of the wave- 
function (with electron correlation), yield far better energies. The point 
missed by Coulson is that for larger systems these other methods quickly 
become much more cumbersome and expensive than HF. Roothaan's work 
came at just the right time. Application of the BSE approach for larger 
molecules depends upon electronic computers, the development of which 
was just starting in 1951. 



*Z. Physik 61, 126 (1930). 



^C. C„ J. Roothaan, Rev. Mod. Phys. 23, 69 (1951); G. C. Hall, 
Proc. Roy, Soc. A205, 541 (1951). 
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§3. 5. 2d Interpretation, Correlation, and the Koopmans 1 Theorem 

The variational condition (20) has the form of a Schrodinger 
equation for a particle moving in the potential 

_ 1 i 

V "" F~~ FT + J d 
a d 

(assuming for the moment, H 2 ). 

This is just the classical potential that would be obtained if the second 

electron were replaced by its classical potential, J . . Thus the Hartree- 

Fock orbital is the eigenstate of the motion of an electron in the 

average potential due to the other electron. Indeed Hartree originally 

derived his equations from just such classical considerations. In the 

above derivation however, we applied the variational principle and 

found that the best possible orbital satisfies such an equation. 

Of course, in the real molecule the electron motions will be 

such as to keep — as small as possible and — and — as big as 

Fl2 r a r b 

possible, while also minimizing the kinetic energy. Thus at instants for 

which one of the electrons happens to be close, to the left nucleus, we expect 

that the other electron will tend tobenear the right nucleus. Such instantaneous 
correlations in the motion of the electrons are ignored in the Hartree- 

Fock wavefunction. Both electrons move in the same orbital independently 

of the instantaneous position of the other electron. Hence the error in 

the Hartree- Fock wavefunction is called the correlation error. 

The energy e in (20) is called the orbUalejoergv^. From (16) 

it has the value 

e =(0|.hk) + J # 
= E, - E, 
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where 

E 2 = 2<0|h|0) + J^ 

is the energy of the two electron system with both electrons in <j> while 

E x = <0|h|0> 

is the energy with only one electron in $ . 

Thus e is just the negative of the ionization potential . There are two 
approximations here. One is using the Hartree-Fock energy for the 
two-electron molecule, leading to too high an energy for the two-electron 
system. The other error is in describing the ion with the optimum orbital 
$ found for the two-electron molecule, leading to too large an energy for 
the ion. These errors often tend to cancel yielding ionization potentials 
within ~ 10% of the exact value. This approximation* of using the orbital 
energy to approximate the ionization potential is often called the Koopmans 
theorem . Although not strictly the theorem that Koopmans proved 

[T. Koopmans, Physica, 1, 104 (1933)] •, we will also refer to this 

t 

approximation as Koopmans Theorem . * 



Note that there is an s at the end of this name and that the Dutch oo 
sounds about like our long o. 

^Shortly after this work Koopmans switched to economics. He is on the 
Economics faculty at Yale and in 1975 won a Nobel prize for his work in 

optimization theory in economics. 

t 

+ The first application of this approximation was by D. R. Hartree, 

Proc. Camb. Phil. Soc. 24, 111 (1927). 
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§3.5.3 The Generalized Valence Bond Method 

In §2.3 we described the simple VB wavefunction of fi> in which 
the two-electron wavefunction is expressed as: 

$(1,2) =(x^X r + X r X £ ) , (21) 

where \* and x r are atomic orbitals. 

We will now consider wavefunctions of the VB form (21) 

*(l,2)=* a « b+ ^ a (22) 

but where the orbitals <£ a and $ b are allowed to be completely general. 
To obtain the best such orbitals we will apply the variational principle, 
requiring the orbitals to lead to the lowest possible energy. The optimum 
orbitals are called the generalized valence bond orbitals and the resulting 
wavefunction is called the gej^aUzjedj^aJence^^ (GVB) wavefunction. 

§3. 5. 3a The Basic Equations 

The energy of the GVB wavefunction (22) is 

E=g + I (23) 

where 

N =[(a|h|a) + <b|h|b> + J ab ]+[2(a|h|b)S ab + K ab ] (24) 

D= 1+S ab 
S ab = < a ! b >- 

Here and in most of the following we will use orbital subscripts (e.g., 
a and b) to denote orbitals (e.g., <£ a and <£ b , respectively). Just as in 
§3. 5. 2 we will consider that a and <p, are expanded in a basis { x } 
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^a = I>V C 



jua 



a 






and require that 



(26) 



#: ; =0 < 27a > 

juia 

3^- =0 . (27b) 

jiib 



condition (27a) leads to 

<X M l(H a - ea )k a >=0 , (28) 

where 

e a =E - <b[h|b> (29) 

and H a is an operator taking care of all other terms. To solve (28) we 

expand (p with (26) leading to 
a 

Z<X /x l(H a -e a )| V C Mi = , (30) 

V 



* The form of H a is 



H a = (h + J b + K^) + P b h + h P b - EP b 

where P b s |b)(b| is a projection operator. However, the explicit form 
is of no importance here. 



3.5-16 



which in matrix notation is 



H a C = £ SC (31a) 

- -a a a ' v ' 



where the elements of the H a and S matrices are 



H a , = <X fl |H a |x > 



ill ii (32) 

<x M |h|b><b|x l/ > -E<x /1 |b><b|x i ; 



and 



For the HF wavefunction the corresponding matrix is 

H = <x (h + JJx > (33) 

The presence of the second (exchange) term in the wavefunction (22) 

leads to the other terms in (32). Of these additional terms only the K fe 

term would be present if the orbitals c£ a and b were orthogonal. 

Equation (31a) is the condition for orbital (p to be optimum, 

a 

there is a similar equation 

5 b C b = e b SC b (31b) 

to solve for the optimum orbital <&. Thus in the GVB method we must 
solve self- consistently for two orbitals whereas in the HF method we 
have just one orbital and hence one equation to solve. Otherwise both 
involve similar computational procedures. 

As shown in Fig. 8, the GVB orbital of H 2 corresponds closely to the 
VB orbital. 
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§3, 5.3b Interpretation 
If (28) is sati 
set, then <p satisfies the differential equation, 



If (28) is satisfied for all basis functions x of a complete 



The H a in (34) can be written as 

H a = h + v b (35) 

where V b contains all the terms depending upon orbital # b . We can con- 
sider V, as the average potential* seen by the electron in <p due to the 
D a 

electron in 0,. The average potential is not just the Coulomb potential J fe 
due to the electron on <£, (as would be expected classically), but also con- 
tains other terms arising from the quantum mechanical form of the wave- 
function. However, these additional terms are not chosen arbitrarily and 
indeed were determined through the variational principle as the terms 

required in order that (6 be the optimum orbital to place in the two-electron 

a 

wavefunction. Thus we can consider the potential V, as the quantum 
mechanical generalization of the classical Coulomb term for the 
interaction between electrons in overlapping orbitals <£ a and $,. 

The operator, H a , in (35) is equivalent to the Hamiltonian 
for an electron moving in the potential due to the nuclei (contained in h) 

plus a potential V fa due to the electron in orbital ck . Since orbital 

a 
<f> is an eigenfunction of H , we can interpret (p as the eigenstate of 

an electron moving in the average potential (V,) due to the other electron 
of the system. Similarly, of course, we can interpret <h as the eigen- 
state of an electron moving in the average potential (V ) due to the other 

a 

electron. Thus with this interpretation we can describe the two-electron 

*Note well that V^ is not a local potential (that is, a function of r). Rather V b 
contains integral operators and upon operating on a , puts a inside an integral. 
Even so we can consider V b as the effective potential due to cj> h asseenby </> 
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system in terms of two one-electron systems each of which contains 
the average potential due to the other electron. Such a description 
of a multi-electron system in terms of electrons moving independently 
of each other will be termed an independent particle interpretation (IPI). 
We will find, especially for larger molecules, that such IPPs will be 
very useful in understanding the wave functions. 

It is important to note here that the IPI comes about from the 
one-particle Schroedinger equation, such as (34), arising from applica- 
tion of the variational principle to a special type of wavef unction (22). The real 
electrons of a molecule are quite indistinguishable, and we do not imply that one 

of the electrons moves in one orbital (say 6) while the other electron 

a 

moves in the other orbital (<£,). What we say is that the orbitals satisfy 
a one-electron Schroedinger equation for which the field term is the 
average potential of an electron in the other orbital. This is not the 
equation describing the motion of one of the real electrons. However, 
considering the eigenstates of two fictitious, distinguishable electrons, 
we do obtain the optimum orbitals for the many-electron wavefunction (22). 
It is really the orbitals which are distinguishable here, not the electrons. 

§3. 5. 3c GVB Natural Orbitals 



In order to obtain another view of the GVB wavefunction we will define 
the GVB natural orbitals* <j> and <p as the sum and difference of the GVB 
orbitals <p t and 4> T 

g g (36) 

« u = «Y-^)/d u , 

where 
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D = V2(l + S) 
g J (37) 

D u = V2(l - S) 
and S is the overlap of the GVB orbitals 

Rearranging (36) leads to 

20 =D + D (b 
^r g^g u^u 

^£ g g umi 
and hence the total wavef unction becomes 

Thus we may view the GVB wavefunction in terms of fa and 6 where there 

is always one electron in fa and one in 6 or one may view this wavefunction 

in terms of GVB natural orbitals 6 and <f> where part of the time both 

g u 

electrons are in <p and part of the time both are in fa . Here fa resembles 
the bonding orbital and <j> the antibonding orbital. The equivalence of these 
two descriptions may be clear in Figure 9, where Figure 9d and Figure 9g 
are equivalent. The first natural orbital fa has a good kinetic energy but a 
bad two-electron energy. Mixing in a small amount of <p causes an increase 
in the kinetic energy, but this is more than compensated by the decrease in 
the electron- repulsion energy, leading to the optimum wavefunction Fig. 9d 
or g. 



(a) 
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Figure 3-9. 
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Requiring that the total wavefunction be normalized leads to 

= C g Vg " C u^u (38) 



where from (37), 



2 



C u D u (1 - S) 

U U (39) 



v2 



C g D< g (l + S ) 

That is, if the overlap between the two orbitals is nearly zero (Ha for 

R -* °° ), then the two configurations come in with nearly equal coefficients. 

On the other hand, for H 2 at R = 1. 4 a S = 0. 8 and hence 

Eh. * o.n 
c g 9 



*In this section we use g and u for the orbitals as appropriate for Ha. However, 
the discussion does not depend upon inversion symmetry and all results apply 
also to a GVB pair for an unsymmetric system . 
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§3.5.4 Electron Correlation 

In a real atom the electrons are expected to move somewhat in concert 
so that they avoid getting too close to each other while remaining close to the 
nucleus. That is, their motions are somewhat correlated. On the other 
hand, in the HF wavsfuncticn, 0(1) (f)(2), each electron is placed in the same 
orbital and hence the probability of either electron being at a particular 
position is independent of where the other electron is. That is, the electrons 
in the HF orbital are uncorrected in their motions. For this reason the 
difference between the HF energy and the exact energy is called the 
correlation error. For the ground states of two-electron atoms (H~, He, 
Li + , etc.) the correlation error is about 1.1 eV. In addition, for H 2 (at R ) 
the correlation error is 1.1 eV. Although this correlation energy is small 
compared to the total energy of these systems (e.g. , 1.5% for He), it is 
comparable to many quantites of interest. 

In the GVB wavefunction for H 2 

<^(!) * r (2) + * r (l) *j(2) , (40) 

one electron is in 0. while the other electron is in <£ r regardless of which 
electron is in which. Hence, there is static correlation in the sense that 
the orbitals for each electron are in a slightly different region of space, 
and hence on the average the electrons stay farther apart. However, the 
presence of electron 1 at a particular location of the </>- orbital does not 
affect the probability of electron 2 being at any particular position in orbital 
, and hence we may consider that the GVB wavefunction does not provide 
for instantaneous correlations among the motions of the electrons. Since the 
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GVB wavefunction is the most general wavefunction involving just two spatial 
orbitals, we may consider that all correlation error beyond GVB involves 
instantaneous correlation of the electrons. When important to distinguish 
these effects we will refer to the latter as dynamic electron correlations and 
the difference between HF and GVB as static electron correlation. 

Now consider the description of correlation in the natural orbital (NO) 
representation of the GVB wavefunction 

*(i?2) = C g <V 1} <V 2) " c u KM <V 2 > • («) 

Assume that electron 1 is at some position R on the right side of the molecule 
and consider the likelihood of electron 2 being at equivalent positions R or L 
on the right and left sides of the molecule. In the HF wavefunction 

$ HF (R,R) = <p(R) tf>(R) 
<D HF (R,L) = 0(R)«(L), 



and since 



<MR) = <ML), 



we have equal probabilities 



$ HF (R,R)| 2 = |$ HF (R,L)| 2 (42) 



of the electrons being on the same or opposite sides 
In the GVB wavefunction (40) we find 



* GVB (R, R) = ^ (R) r (R) + cp r (R) ^ (R) 
3> GVB (R, L) = <fy (R) r (L) + r (R) fy (L) 
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and hence 

<S> GVB (R,R) <$ GVB (R,L) , (43) 

that is, we obtain the static correlation referred to above. Using the NO 
form of the wavefunction we obtain 

$ NO (R,R) = C g 4> g (R) 4> g (R) - C U 4> U (R)0 U (R) 

<I> NO (R, L) = C g g (R) g (L) - C u U (R) u (L) 

= C g g (R)0 g (R) + C U U (R)0 U (R) 

(using the symmetries of $ and ) and hence 

g u 

<$ NO (R,R) < $ NO (R, L). (44) 



Just as in (43). Comparing (41) with the HF wavefunction 

# HF (1,2) = g (l)* g (2) (45) 

we see that in order to obtain effective electron correlation, the second NO 



must have a shape similar to that of the first (dominant) NO but with an extra 
nodal plane bisecting the first NO. This allows maximal difference between 
* (R,R) and * (R, L) and hence maximal electron correlation. 
We will later find such arguments in terms of nodal planes to be useful in 
describing other electron correlation effects. 

The above discussion should be made clear in Fig. 9. 
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§3 .5. 4a Ionization Potentials 

In general we expect the correlation error to increase with the number 
of electrons (since there are more and more complicated interrelationships 
ignored). Thus the ionization potentials calculated from HF and GVB should 
be too small. On the other hand we can get an approximate IP from 
Koopmans' theorem. The Koopmans IP is the energy difference between 
the self-consistent energy of the N-electron system, E N , and an energy of 
the N-l electron system, E N _j, obtained using orbitals from the N-electron 
wavefunction. Thus the description of the ionic (N-l electron) state is non- 
optimum leading to too high a value for E N _ j and hence too large a prediction 
of IP. However, the IP calculated using self consistent wavefunctions of the 
N and N-l electron systems should be too small. Hence there is a cancelling 
of errors such that the Koopmans theorem value of IP is usually rather good 
(within ~10%). These effects are indicated in Fig. 10. 
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§3.5.5 Configuration Interaction Wavefunctions 
Starting with a one-electron basis 

{ XjLt ; M = 1,2-..P} 

we can construct a two-electron basis 

{x^iJx^a); h,i/ = i,2,...p} 

by combining all products of one-electron functions. In terms of this basis 
the most general wavefunction is 



$(1,2) = S C x u (Dx 7 .(2) . (43) 



The terms in (46) are called configurations and the resulting wavefunction is 
called an configuration interaction (CI) wavefunction. 

Applying the variational principle §3.1, the optimum coefficients for 
(46) are solutions of equations 

J (H W - E) V = (47) 

analogous to (3. 1-7) except that each simple index fj. or v is replaced by a 
combined index a 77 or \iv. For a complete basis (P = 00) the resulting CI 
wavefunction is the exact wavefunction of the system. 
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§3. 5. 5a Permutational Symmetry 

Because the electrons are identical the Hamiltonian must be 
invariant (unchanged) upon permutation (interchange) of the electrons 

3C(2,1) =3C(1,2). (48) 

[Recall that 3C (1,2) = h(l) + h(2) + JL.] 

As a result of this permutational symmetry the exact eigenstates of 3C can 
always be taken as either symmetric 

* s (2,l) =* s (l,2) (49) 



or antisymmetric 



# a (2,l) = -**(!, 2) (50) 



under permutation. 

The proof (see App. A) is quite analogous to that in §2. 1 where we found that 

for a system with inversion symmetry, all eigenfunctions are either g or u. 

Later when we discuss the Pauli principle and spin we will find that 
symmetric spatial wavefunctions <& s are allowed only for singlet (S = 0) 
spin states and antisymmetric spatial wavefunctions ^ a are allowed only 
for triplet (S = 1) spin states. 

§3. 5. 5b The Nodal Theorem 

Next we will show that the lowest state of 3C(1, 2) [assuming 5C is 
symmetric (48)] is always a symmetric wavefunction, \J/ S (1,2), (49). 

As shown in §1.2 the ground state of a system is nodeless, that is, 
the wavefunction of the ground state has the same sign everywhere. For 
a one -electron system this means that 
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cannot be the ground state* whereas 




can. 



The nodal theorem applies also for many electron systems, such as 



3C(1,2) = -§V*- iv^ + v( ri ) + v(r 2 ) + -±- 

1 12 



(51) 



(the proof is exactly as in §1. 2). We will now use the nodal theorem to show 
that the ground state of any two-electron system must be a symmetric wave- 
function. 

Letting & = r^ in an antisymmetric wavefunction 



a 



^ ( r 2>£i) = ~P (£!,&) 



leads to 



* *(*!>&) = -**(£» &) 



and hence 



^(r^r,,) =0, if r x = r. 



(52) 



* except for R = oo, where this state is degenerate with the nodeless state. 



o. 5-28 



[For example, typical one- dimensional cases are illustrated in Fig. 11. 1 
Thus every antisymmetric function has a nodal plane whereas symmetric 
functions need not have nodal planes. Since the nodal theorem implies that 
the ground state is nodeless, then the ground state must be permutationably 
symmetric, (49). Later we will find that a symmetric spatial wavefunction 
must correspond to a singlet spin state and hence the ground state of a two- 
electron system must be a singlet state. 

In the case of a sufficiently singular potential it is possible for the 
lowest wavefunction having a node to be as low as the lowest nodeless wave- 
f unction. However, in three dimensions our Hamiltonian (51) is not this 
singular and hence the inequalities apply. 




a - 



z, a 

a) Antisymmetric Wavefunction 

(triplet state) 




b) Symmetric Wavefunction 
(singlet state) 



Figure 3-11. Illustration of nodal patterns of two-electron systems 
(in one dimension). 
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§3.5. 5c Natural Orbit als 

In App. B we show that the CI wavefunction (46) for the ground state 
of any two-electron system can be written as 



*°(1,2) = £ C X (1) x (2) (53) 

/i=l ^ ^ ^ 

(that is, double occupied orbitals only) where the orbitals {x }, called 

*■ ji 

natui^or^itals, are orthonormal 

( Vv ) = V • (54) 

Since (53) has only P terms rather than P* as in (46), it is obviously easier 
to interpret. 

The density of electrons in a two-electron system is defined as 



p(l) = / d 3 r 2 *(1,2)**(1,2). 
Thus using (53) leads to 

p(l) = S C*C X (l)*x (1) (X (2) I X (2)) 

= S Pfi | x „(i)| 2 

where p = |c I . 

Since there is a total of two electrons in the system 

fd*r x p(l) = 2, 
the coefficients must sum to two, 

Tj p = 2. 

Consequently, in terms of natural orbitals, the total density of the CI wave- 
function is just the sum of the densities of the natural orbitals weighted by a 

population p that sums to two. 
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§3.6 Wavefunctions for He 

In this section we will illustrate the HF, GVB and CI methods by 
describing the wavefunctions for He atom. 

§3.6.1 HF Wavefunctions for He Atom 

First we consider various approximations to the HF wavefunction 

<£> HF (1,2) = 0(1)0(2), 

where the HF orbital (p is expanded in a basis set. 

§3. 6. la MBS 

The simplest description of He atom is to place two electrons in the 
Is orbital of He + 

X = e"^ r (1) 

where £ = 2.0. The total energy in hartrees is just 

E = 2e x + J = -2.75, 

where 

€ x = -2.0 

is the energy of He + and 

tl — ' ~Q Qj — JL a Civ , 

is the Coulomb interaction of the two electrons (see App. C). 

This description can be improved by optimizing £ (leading to the 
MBS description). As shown in App. C the energy has the form 
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E = ? 2 Ti + ?V L , (2) 



where 

T x = (2) (|) = 1 
V l = -2Z + | 

(i.e. , Tj_ and V x are the kinetic and potential energies for the case of 
C =1). Requiring that 3E/3£ =0, to obtain the optimum £," leads to 

= 1.6875 
for He. Since the optimum £ for the one-electron atom is 

C = z, 

we can interpret the C^pp as an effective charge that has decreased from Z 
due to the presence of the second electron. It is as if the second 
electron partially shields the nucleus; hence the quantity 

5 

a = 



16 

is sometimes referred to as the shielding correction. 
Substituting (3) into (2) leads to 

2 __ 

E OPT = " ri\ = + — ^OPT = " ^OPT (4) 

= -2.84766 h . 

Xhis energy is the same as if there had been two non- interacting electrons 
each experiencing the Coulomb field due to a nucleus of charge 

? OPT = z " Jig • 
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The exact energy for He atom is 

-2.9037 
Thus the above simple wavef unction accounts for 98. 5% of the exact 
energy! Since the correct energy of He + is 

-2.0 , 
the use of the calculated E of (4) leads to a predicted IP of 0. 84766 or 
94% of the exact value. Use of the Koopmans theorem leads to 

-IP = <4> | h | 0> + J^ = \ C^pr - Z C OPT + I £ OPT 

= -k ^QPT + U ^OPT = " °' 89648 » 

99. 2% of the experimental value (a better value is obtained because we 
describe the ion badly). 

§3, 8. lb Bigger Basis Sets 

The results of using various-sized basis sets for HF calculations 
on He are listed in Table 2. In the cases P = 1, 2, and 3 extensive 
optimization of the parameters was carried out, leading to quite short 
expansions. Thus with P = 2 we are within 0. 000007 h = 0. 0002 eV = 0. 005 kcal/mol 
of the HF limit (P = »). With P = 3 the energy is correct to 8 decimal 
places (comparing to the HF limit). Bear in mind here that the exact 
(nonrelativistic) energy for He is -2.903 so that even the exact HF 
wavefunction is off by 0. 042 h = 1. 1 eV = 25 kcal/mol. 

The HF orbitals in these various approximations are plotted in Figure 15, 
Note that even though P = 3 and P = 12 lead to the same energy (to 6 decimal 
places) there are still noticeable changes in the orbitals. 

The conclusion here is that two suitably chosen basis functions are 

adequate for describing He. Such a basis is referred to as double^zeta (DZ) 
or double valence (DV). 



Table 3-2. Parameters for HF wavefunctions of the ground state of He. E is the total energy, 
e is the orbital energy. The orbital exponents are shown in parentheses while the 
expansion coefficients are not. All quantities are in Har tree. atomic units. 



P 
E 
e 



2.847656 
0.89648 



2 a 
-2.861673 
-0.91792 



ns(0 C. ls(l. 68750) 1.0 



ls(1.453) 0.84289 
ls(2.906) 0.18159 



3 a 
2.861680 
•0.91795 



ls(1.450) 1.36211 
2s(1.723) -0.28189 
2s(2.641) -0.10724 



12 b 
-2.861680 
-0.91796 



ls(3.0 
ls(1.4 
2s(3.0 
2s(1.4 
3s(3.0 
3s(1.4 
4s(3.0 
4s(1.4 
5s(3.0 
5s(1.4 
6s(3.0 
6s(1.4 



0.45742 
0.00000 
0.24427 
0.12985 
0.13657 
0.11340 
0.09451 

-0.08686 
0.00819 
0.02546 
0.02767 

-0.00267 



P. S. Bagus, T. L. Gilbert, H. D. Cohen, and C. C. J. Roothaan, unpublished, 1966. 
C. C. J. Roothaan, L. M. Sachs, and A. W. Weiss, Rev. Mod. Phys. 32, 186 (1960). 
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Q 

3 



CL 
< 



.6 



.4 



.0- 



0.8 



0.6- 



0.4- 



0.2- 



0.0 



-4.0 



p= 12 



He 
HF 




-2.0 0.0 2.0 

Z-C00RDINATE (BOHR) 

Figure 3-12. Comparison of Hartree-Fock 

wavefunctions for Helium using 
various (optimum) basis sets, 
p indicates the number of func- 
tions in the basis set. 
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§3.6.2 The GVB Wavefunctions for He 

For He atom the GVB wavefunction is the optimum wavefunction of the 
form 

S GVB (1,2) = * a (l)* b (2) + * b (l)* a (2) (5) 

The GVB orbitals of He are shown in Fig. 13, where they are compared 
to the Is orbital of He + and to the HF orbital of He. We see that <£ GVB 

, „ staU „ «„ «. ,. — o, .... »« ta , v u M JL «-.. 

Thus we envision He as having (i) one electron in the Is orbital of He + (this is 

orbital 4> ) experiencing an effective nuclear charge of Z ~ 2 and (ii) the 
a 

second electron in an orbital (#, ) experiencing an effective charge of Z ~ 1 

(nuclear charge of 2 but shielded by the <J> electron). 

a 

This type of correlation is referred to as in- out correlation since when 
one electron is closer to the nucleus, the other tends to be farther away. 

This GVB picture is somewhat different from the HF model where both 
electrons are in the same orbital and one cannot relate the description so 
simply to that of He + . 



Describing both <j> and <p, as simple exponentials and optimizing the 
exponents leads to effective charges of 

£ b = 1.189 

r = 2.183 
a 

as expected from the simple picture. [J. N. Silverman, O. Platas, and 
F. A. Matson, J. Chem. Phys. 32, 1402(1960).] 
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A more extreme case is H~, the GVB orbitals for which are shown 

in Fig. 14. Here the first electron (in <p) is very similar to a hydrogen Is 

a 

orbital and the second electron is barely bound, leading to a very diffuse <£, 
orbital. As shown in Table 3, the HF wavefunction for H~ yields an energy 
of -0. 487, higher than the energy of the hydrogen atom (implying that H" is 
not stable with respect to H plus an electron). The second electron cannot 
leave since the HF orbital is doubly occupied; thus either both electrons stay 
or both leave. The GVB wavefunction yields an energy of -0. 513, correctly 
accounting for the stability of H~ (the exact energy is -0.527). 

§3. 6. 3 CI Wavefunctions for He 

The results of several CI calculations on He are shown in Table 4. 

Analyzing in terms of natural orbitals leads to the results in Table 5. 

Here we see that the 2s, 2p , 2p and 2p natural orbitals provide the 

x y z 

dominant electron correlation effects. These are the only natural orbitals 
containing just one nodal plane. Plots of the 2s and 2p NO's are given in 
Fig. 15, where we see that the higher NO's are concentrated in the same 
region as the Is orbital but with the additional nodal plane (circular for 
<£~ and planar for each ^n^ 
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Table 3-4- Energies for CI wavefunctions of the ground state of He atom. 



Number of Basis Functions' 
s p d f g 



Energy 



5 

5 4 

5 4 3 

5 4 3 

5 4 3 

Pekeris 

HF 

GVB 



-2.87896* 

-2.90039 a 

-2.90258 a 

-2.90307 a 

-2.90320 a 

-2.90372 b 

-2.86168 

-2.87800 



b 



A. W. Weiss, Phys. Rev. 122, 1826 (1961). 
C. L. Pekeris, Phys. Rev. 115, 1216(1959). 



Table 3-5. Analysis of He CI wavef unction in terms of Natural Orbitals. 



Natural 
Orbital 



Energy Lowering 
Millihartrees 



% of Total b 

Correlation Energy 



2s 


16.30 


2p 


19.51 


3s 


0.88 


3p 


1.63 


3d 


1.80 


4s 


0.09 


4p 


0.26 


4d 


0.36 


4f 


0.35 


Totals 


41.18 







46. 


46 


2. 


09 


3. 


87 


4. 


28 


0. 


22 


0. 


62 


0. 


86 


0. 


83 



98.05 



N. SabelliandJ. Hinze, J. Chem, Phys, , 50, 684(1969), 



b 



Total correlation energy = 0. 0420 hartree. 
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NATURAL ORBITALS FOR He 




0.4 




-6.0 



-3.0 0.0 3.0 

Z-COORDINATE (BOHR) 



6.0 



Figure 3-15. The natural orbitals for He 
(Sabelli and Hinze, loe. cit . , 
Table 4). 
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Using the five dominant NO'S: Is, 2s, 2p x , 2p 2p z , leads to the 



wavefunction 



CI 

* = C ls *ls *ls " C 2s *2s *2s " C 2pl V V + V V + V V 1 (6) 

xx y y z z 



where 



C lg = 0.99599 

C 9o = +0.06160 
zs 

C 2p = +0 - 06188 • 



This wavefunction has an energy of -2. 8975, just 0. 17 eV above the exact 
(nonrelativistic) energy of -2. 9037. This wavefunction is sufficiently 
accurate that for the purposes of this course we will consider (6) as the 
exact wavefunction of He. 

§ 3. 6. 3a Interpretation of the CI Wavefunction 

To interpret the wavefunction (6) we will consider one by one the 
effects of adding any one of the four correlating terms to the dominant 
(first) term. 

The wavefunction 

C ls*ls*ls - C 2s*2s*2s < 7 > 

can be rewritten in a GVB form 
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where 

^a =l/ Vls + VV2s 

(9) 

v^0 b = vc7;0 ls -vc^0 2s 

are GVB type orbitals. Thus from §3.6. 2 we see that the 2g natural orbital 
in (7) builds in r^diajjsoj^ejation, increasing the probability of the second 
electron being at larger r when the first electron is at smaller r (and vice versa) 
Similarly the wavefunction 



C ls Vis - C 2p%/2p Y 
can also be rewritten as (8) where 

V2> a = vct;* 1s + ^ip> 2p 

x (10) 

In (10) we see that when one electron is in the +x direction, the other tends 
to be in the -x direction. Similar results occur for the 02p and 02 p terms. 
The three correlations resulting from the terms involving p orbitals are 
grouped together and referred to as angular correlation. 

The three terms of (6) involving p orbitals can be written as 

^2px*2px + ^2 P y *2py + *2pz*2pzl = 

= R(l)R(2)[sin0 1 cos cp 1 sinO, cos cp 2 + sii^sinc^ sin0 2 sin<^ 2 

+ cos d 1 cos 9 2 ] 

= R 2p (1)R 2p (2) cos ^2 9 (11) 
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where R(i) is the radial part of the orbital 0. and 12 is the angle between 
electron 1 and electron 2. Combining with the first term of (8) we obtain 

C ls R l S ( : ) R ls( 2 ) - C 2p R 2p (l)R 2p W COBfl »- < 12 ) 

With this form we see that the magnitude of the wavefunction is increased 
(with respect to <£ ls 0i s ) *o r | e i 2 | > 90 ° and decreased for |fl 12 | < 90°. 
Thus (12) effects an angular correlation of the electrons. 

Each of the four dominant correlating orbitals has one nodal plane 
not contained in* <f>^ and the correlation effect is across this nodal plane 
(increased probability of electrons being on opposite sides). Starting 
with the $ ls orbital, there are just four possible orbitals orthogonal to 
ls but containing a single nodal plane, namely the above four. All 
additional correlating terms will involve two or more nodal planes 
(leading to higher energies) and all are relatively unimportant, leading 
to a total energy contribution of 6. 2 mh = 0. 17 eV = 3. 9 kcal. For the purposes 
of most of our considerations of molecules an energy error of 0. 1 eV 
is acceptable, and we will completely ignore these smaller terms. 
Thus we will consider (6) as the CI wavefunction of He. 



*Of course 0^ has no nodal planes; however, we have worded this so 
as to be appropriate also for correlations of more complicated orbitals 
than 4>i S ' 
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§3.7 Wavef unctions for H 2 

In this section we will discuss the HF, GVB, and CI wavefunctions 
for H 2 . 

§3.7.1 HF Wavefunctions for H 2 

In Figure 16 we show how quickly the HF wavefunctions for H 2 converge 
as a function of basis set size (P). The major effects in the orbital shape are 
in the bond region. 

In Figure 17 we compare the MO wavefunction (P = 2, £ = 1) with the 
MBS wavefunction (P = 2, £ opT ). Here there are significant changes near 
the nuclear and bond regions. 

Comparing the energies in Table 6 we see that P = 6 leads to an energy 
within 0. 00152 h = 0. 04 eV = 1 kcal of the HF limit. We consider this as a 
good level of accuracy. The P = 6 basis has two (optimized) s functions **>■ 
on each H and an (optimized) p function on each H. Such a basis is referred to 
as double valence (for the two sets of s functions) plus polarization (for the 
p functions) and will be denoted as DV-P. 

With even the best of these HF wavefunctions, the energy is 0. 04081 h = 
1.1 eV above the exact (nonrelativistic) energy of H 2 , about the same as the 
correlation error of He (and other two-electron ions). 

The HF potential curve using the P = 6 basis of Table 6 (optimized at 
each R but restricted so that £ lg = ? 2s ) is shown in Fig. 18. Just as with 
the MO wavefunction, the HF wavefunction at large R leads to very serious 
errors. Thus at R = 6 a with the P = 6 wavefunction, the energy is 



*More commonly, double valence is referred to as double zeta. 
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Figure 3-1 6, Hartree-Fock wavefunctions for 
H a with R = 1 . 4 ao . p indicates 
the number of functions in the 
basis set. 
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E = -0.82199 (already far above the dissociation limit, E = -1.0), and the 
orbital energy is e = -0. 32170 (way off from the correct value at large R, 
e = -0. 50). For R = *> the HF wavefunction leads to an energy of -0. 71542 
which is 7.744 eV above the dissociation limit.* 

§3.7.2 The GVB Wavefunction for H, 

The GVB wavefunctions and energies for several (optimized) basis sets 
are given in Table 7. A quite adequate description (0. 2 kcal from the limit) 
is obtained using a single (optimized) s function and a single p function on 
each center. Even the MBS is only 4. 1 mh = 0. 11 eV above the limit. 

The GVB orbitals as a function of R are shown in Fig. 19. At large R 
the orbitals are atomic-like, but for smaller R the GVB orbital gradually 
becomes more contracted about each nucleus. These readjustments in the 
orbitals are such that the contragradience in the bond region is about the 
same as for the VB wavefunction. From 1 to 6 a the GVB orbitals lead to 
a much greater overlap than the VB orbital as shown in Fig. 19. For example, 
at R = 1 . 4 a , S GVB = 0. 804 as compared to S VB = 0. 753. 

§3. 7. 2a Energy Analysis 

The GVB orbitals for H 2 are compared with the VB orbitals in Fig. 21, 
where we see that the orbitals readjust in such a way as to maintain the large 
contragradience in the bond region while concentrating the orbitals more about 
each nucleus. The GVB energy curves are compared with other energy curves 
in Fig. 22. 



*W. A. Goddard, J. Chem. Phys., 48, 5337 (1968). 
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Table 3-7. Energy and wavefunctions for GVB calculations on K, at 1. 4 a . 
Only the 0^ orbital is given, the <6 orbital is the mirror image. 
The same basis occurs on both centers with the orbitals on 
the left first (the basis functions on the right have no exponent 
listed). A p basis function with + coefficient is positive 
toward the second center. All quantities in hartree atomic units. 



E 


-1.147777 


€ 


-0.6877 


s 


0.79700 



1.151345 

0.68348 

0.80093 



6 U 
-1.151526 
-0.68472 
0.80420 



ns(?)C ls(1.2005) 0.91287 

Is 



ls(1.1909) 0.88890 ls(1.3129) 0.77499 



0.12303 


2p J2.0928) -0.00672 
z 


2s(l. 


,1566) 


0.11116 




Is 0.13631 


2p z (1.9549.) 


-0.00310 




2p„ 0.03006 


Is 
2s 

2 P. 




0.12161 
0.04199 
0.03769 



Using ls( 1.262) and 2s(l. 191) basis functions on each center leads to 
E =-1.147804. 

b Using ls(1.3092), 2s(1.1273), 2p (1.700), 3d 2 (2. 37) basis functions on 

z z 



each center leads to E = -1. 151887. 
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Figure 3-20. Comparison of overlap, S = 
<0£ | cf> r ) for the VB and GVB 



wavefunctions. 
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Some of the energy parameters of the HF, VB, and GVB 
wavefunctions are compared in Table 8, while the ? t as a function 
of R is given in Fig. 23. For 5 = 1. all three wavefunctions yield 
an R e far too large (14% to 19%). Optimizing £ leads to errors of 
only 1% to 2% in R and improves the calculated bond energies by ~20%. 
It is characteristic that GVB leads to too large an R while HF leads 
to too small an R. 

Using the form 

*a*b + *b*a (1) 

for the GVB wavefunction, we can define classical and exchange terms 
much as for the VB wavefunction 



E GVB s< Vr |3C KV 
iL GVB =£j GVB" ^GVB > 

etc. This leads to the results in Fig. 24, where we see that the exchange 
term still dominates the bonding. In particular for R>R e the 

E x is very nearly the same for VB and GVB. Thus the main improvement 

Co- 
here is in the classical term, E . Similarly, in Fig. 25 we see that 

it is the exchange part of the kinetic energy that dominates the bonding 

energy. Again for R >R^ we see only minor changes in T* between 

e 

VB and GVB. 
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Figure 3-21. Comparison of GVB (two basis functions, ? ,= 1. 2005) 
and VB (£= 1) orbitals for H 2 (R = 1.4a ). 
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Figure 3-22. Comparison of the energy curves for MO, HF, VB, and GVB 

wavefunctions of HLj. Only two basis functions were used for HF 
and GVB. The results for both ? = 1 and ? . are shown. (The 
C op t as a function of R are given in Fig. 2 6.) The energy is 
relative to the energy of two hydrogen atoms. 
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Table 3-8 . Comparison of results on H 2 for approximate wavefunctions 

using two basis functions. All quantities are in atomic units; 
the energies are relative to two hydrogen atoms at R = °° . 



HF 



VB 



GVB 



Exact 



5 = 1. 





R e 


1.603 






e 


- 0. 0990808 


^OPT 




R e 


1.385 
-0.128231 






bPT 


1.1931 


R = 1. 


4a 


E 


-0.128189 






^OPT 


1.1895 



1.643 


1.568 




0.115971 


-0.118651 




1.414 


1.431 


1.401 


0. 139083 


-0.147938 


-0.174470 


1.1661 


1.1937 




0.139049 


-0.147777 


-0.174470 


1.1695 


1.2005 
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Figure 3-23. ? . for the :hF, VB, and GVB wavef unctions . 
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Figure 3-24, Comparison of the classical and exchange energies 



for the VB and GVB wavefunctions of H 2 . 
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Figure 3-25. The kinetic and potential parts of E x for 



the VB and GVB wavefunctions of H 2 . 
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§3.7. 3 CI Wavefunctions for H 2 



In §3. 6. 3 we found that in He there are four important correlations each 
involving a correlating natural orbital having one nodal plane 



2 s radial 

2 Px ) 

2p \ angular 



(2) 



For H 2 the HF orbital is nodeless and again we can find four correlating natural 
orbitals each with one nodal plane. These are illustrated in Fig. 26, where 
the names la la u , etc. will be explained below. 

As R -» 0, the H 2 orbitals in Fig. 26 change smoothly into (we say that 
they correlate with) the He orbitals in (2): 



H 2 


He 


g 


Is 


2a -> 
g 


2s 


^u- 


err 


In - 
ux 


>2p x 


uy 


' 2p y 



(3) 



and hence the correlation effects are closely related: 
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Figure 3-26. Correlating orbitals for H 2 . Long dashes indicate 
nodal planes, solid lines are positive amplitudes 
and dotted lines negative amplitudes. 
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H 2 He 

left-right (la u ) — angular-z (p ) 
starboard-portside (lrr Y ) ~~ angular-x (pj (4) 

UX X 

up-down (l7r uy ) — angular-y (P ) 
in-out (2a g ) - radial (2s) 

The five dominant natural orbitals for H 2 are shown in Fig. 27, which 
should be compared to Fig. 15 for He. With these five dominant NO'S the 
wavefunction 



^V^V'VV'V^/Vk (5) 

- Cl *u [01 V>ux + 'VV 

leads to an energy of -1. 1699 h. Comparing to the exact energy of -1. 1744 h, 
we see that wavefunction (5) accounts for all but 4. 5 mh = 0. 12 eV = 2. 8 kcal 
of the exact energy. This is quite adequate for our purposes and we will 
ignore all other terms. 

A more complete analysis [E. R. Davidson and L. L. Jones, J. Chem. 
Phys. , 37, 2966 (1962)] of CI calculations on H 2 (for R = 1.4 a ) in terms of 
NO'S is given in Table 9. 

For the molecule at R e the dominant correlation is left-right. This 
becomes even more so for larger R. Thus at R = °° the exact wavefunction is 



*(l,2) = Cl A <p -c la «j t 

g g g u 1 u u 
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NATURAL ORBITALS FOR H 2 (R = |.4a ) 




Figure 3*27. 



The natural orbitals of H 2 for 
R = 1. 4 a . [From G. Das and 
A. C. Wahl, J. Chem. Phys. , 
44, 87 (1966)] . 
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Table 3-9. Analysis of H 2 CI wavefunction in terms of NO. a 



Natural Energy Lowering % of Total , 

Orbital Millihartrees Correlation Energy 

la 
u 

Iff 
u 

2a 
g 

Iff 

g 

g 

2a 
u 

16 
g 

4a g 

Totals 39.43 96.60 



a E. R. Davidson and L. L. Jones, J. Chem. Phys . 37, 2966(1962) 
Total correlation energy = 0.04082 hartrees. 
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where 



*la =( x lsje +x lsr)/^ 2 
*la u = (- x ls4- x lsr)/^ 2 



That is, only left-right correlation is present at R = °° . For R < 0. 8 a 
in-out correlation becomes more important than left -right correlation. 

3.7.3a Notation 

For diatomic molecules orbitals are classified in terms 
of their dependence upon cp [the angle of rotation about the 

molecular axis (z)]. Thus 

a ==> independent of cp 
ir x ==> cos(p 
7r y ==> sin cp 

where cp is referenced with respect to the xz plane. 
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§3.8 Open Shell Wavefunctions 

In §2. 2 we found that the second and third states of H 2 have the form 



$ = <j> <b - <ft <b . 
g u u. g 



(1) 



where the orbitals and <p u are orthogonal. Such wavefunctions (with 
orthogonal orbitals) are referred to as open shell wavefunctions. We will 
occasionally deal with such wavefunctions and will analyze some aspects of 
the wavefunctions here. 

The energies of the wavefunctions are 

E = -g- (2) 



where 



= <* g lv<*ul*u> ± <^ , gl*u><*ulV 

= 1 (3) 



*n = <4> g 4> n \x\<l> g <t> n ±4> u 4> g ) . (4) 

The first term of (4) is 

(^g^J^kg^u) =<g|h|g)<u|u) + <g|g><u|h|u> + <gu|— |gu) 
= (g|h|g) + (u|h|u) + J^ , 
where 

<* g l*u> = (5) 



and the second term is 
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±<*g# u |3e|tf> u tf> > = ±{<g|h|u)(u|g> +<g|u)(u|h|g) +<gu|— |ug)} 



= ±K gu- 

Thus the energy (2) becomes 



E = <g|h|g) + <u|h|u) + J gu ±K gu (6) 

= <g|h|g> +<u|h|u> + (gg|uu) ± (gu|gu) . 



Thus 



E - E = 2IC - _ . (7) 

gu 

Since K > 0, the 3 $ state is always below the # state. 

The above analysis shows that the wavefunctions 

<b cb ± (b d> 
^g^u "u r g 

lead to an electron repulsion energy 

E ± K^„ . 
gu 

Thus the significance of the exchange integral K is that it is the change 
in the energy npon superimposing the exchanged wavefunction <#> <pf 
on 4>W<pW . See §2.2. 

g U 
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Appendix 3-A. Permutational Symmetry, Two Electrons 

Summary: Since the Hamiltonian 3C(1, 2) for a two-electron system 
is invariant under permutation of electrons 

3C(2,1) = 3C(1,2) (1) 

the exact eigenstates of 3C are each either symmetric or antisymmetric 
under permutation. 

Proof: Consider that # is an exact eigenfunction of 

3e(l,2)¥ (l,2) = E ¥ (l,2). (2) 

Renumbering the electrons this becomes 

3C(2,1)¥ (2,1) = E «r (2,l). (3) 

But using (1) in (3) leads to 

K(l,2)* (2,l) = E * (2,l). (4) 

Thus from (2) and (4) both tf (l,2) and ¥ (2,1) are eigenfunctions of 
3C(1, 2), both with the same energy. There are two possibilities here: 
(i) there are two (or more) different (linearly independent) states with 

energy E or 
(ii) there is only one state with energy E . 

In case (ii) it must be that * (2, 1) is proportional to * (1, 2) 

* (2,1) = A* (l,2). (5) 

But interchanging 1 and 2 in (5) leads to 

* (1,2) = x* (2,l) (6) 

and substituting (6) into (5) leads to 
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* o (2,l) = x 2 ^ (2,l). 

Thus 

X = ±1 . (7) 

That is, for a nondegenerate state the wavefunction must be either symmetric 
(X = +D 

* s (2,l) =* s (l,2), (8) 

or antisymmetric (x = -1) 

* a (2,l) = # a (l,2) (9) 

under permutation of the electrons, respectively. 

Assuming now case (i) we define new functions 

* s (l,2) = * (i, 2 ) + *(2,l) 
* a (l,2)s* (i, 2 ).-*(2,l) 



Applying 3C we obtain 

3( * S = E * S 

and hence the exact eigenf unctions of 5C are again either symmetric or 
antisymmetric . QED 
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Appendix 3-B. Natural Orbitals 

A general CI wavefunction for the ground state of a two-electron system 

P 

* s (l,2) = T, c x (l)x y (2) (1) 

can always be rewritten in terms of doubly-occupied orbitals 

* s (l,2)= S^yD^), (2 ) 

where the natu^l^oj^itals { x } are linear combinations of the original basis 
functions {x }. 



Proof: 



Since \p s is symmetric, the coefficient matrix is symmetric. 



C = C 

If we choose new basis functions {% ; M = 1? . • . P} that are linear 
combinations of the old basis function {x .,; M = 1, . . . P} 



x = y v x 



V 

then the wavefunction (1) becomes 



* b = V c V x (1)V X (2) 
jLu/,077 

aT7 



where 
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C = V y C V 
ot) ^ lie \lv vr) 

The wavefunction * s is unchanged by this transformation of the basis, 
but in the new basis the CI expansion coefficients are different. 
In matrix notation the new coefficients are given by 

C =VCV (3) 

Since C is a real symmetric matrix, there is always some transformation 

V for which the transformed matrix C is diagonal. Thus there is always 

a particular choice of basis functions such that 

P 
* S = I 5^X^(1)^(2) (2) 

M=l 

With this basis there are only P terms in the CI expansion rather than 
P as in (1). Thus (2) is a much simpler wavefunction. To find the V 
leading to the natural orbitals we must first solve the CI equations to 
find C. Hence the natural orbitals do not help us solve for the CI wave- 
functions. However, having obtained a CI wavefunction, we will immediately 
transform to the natural orbitals in order to discuss and interpret the 
wavefunction. 
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Appendix 3-C. Evaluation of Energy Quantities For Atoms 



Here we consider the evaluation of the various energy quantities 
for a two-electron system with both electrons in the same Is orbital, 



* ls =Ne 



-Cr 



(1) 



where the orbital exponent C is variable. 
3.C.1 One Electron Quantities 

First the normalization coefficient, N, is obtained from 



1 = <0 1<8 l* 1q > =N 2 Jdxdy dz e" 2?r 



OO 



2 r 2 -2Cr 

= 4ir N J r dr e * 



4ir 



OO 



(2?) 



2 f * -o 

3- N J p dpe H 



8tt 
8C 3 



N 



so that 



N=, 



(2) 



The nuclear attraction terms are 



V^-<ls|-f |ls> = 



-ZN 2 / r dr e" 2 ^ r J J sin d d<p 







ZN 4?r 
2 — 



(3) 
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The kinetic energy term is obtained most simply by noting that 

zt , is = -^isK 

(where e is a unit vector in the r direction), and hence 

T ls~ < ls l f l ls > = < ls l-i v2 Us> =i<v0 lg - vtf> ls > 

= i^ 2 <0 ls |0 ls >=ic 2 (4) 

To check these quantities consider He + where £ = 2. In this 
case (3) and (4) lead to 

E ls = T ls + V !s = *?* - ? Z 
where Z = 2. Optimizing £ leads to 

and hence 



2 
= _i»Z 

Is 



E«_ = -|Z 



both of which are correct. 
3.C.2 Two Electron Quantities 

For He we also need the two-electron interaction term 

J ls,ls = <*ls*lsir7 2 l ls*ls> = <*ls' J lsl*ls> ' (5) 
where 

J ls (l) = /dx 2 dy 2 dz 2 jjL |0 ls (2)| 2 (6) 

is the Coulomb field evaluated at r x due to an electron (called 2) in 
the Is orbital. 
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The complication in evaluating such integrals as J- - is that 

■l s , is 

the integrand of (6) depends on r 12 . The usual solution is to use the 
Laplace expansion 



i r < 

FT, = 2 "TFl P jP (coS 



12 * ' 

4 =0 r> 



where 

*<= r x if r, < r 2 

= r 2 if r 2 < r 2 

and oppositely for r^. With (7), (6) becomes 

°° °° T <T 2 ,2 ^ 2TT 

J ls (1) = S {/ "Tn r 2 dr 2 |* ls (2)| J sinfldGP^cosfl)/ d*} . (8) 
I = ° r > o o 

The integration over 5 is zero unless ft = 0, so that (8) becomes 

00 2 - ri 2 , l2 °° 



J ls (l)=4»/ £^|^ is (2)| 2 = 4^{f /r 2 dr 2 |0 ls (2)| + /r 2 dr 2 |0 ls (2)|} 
o > x o r i 



(9) 



Before proceeding to evaluate J ls (l)» one should notice 
that (9) is a well-known result in electrostatics. The quantity 

Q 1 = 4v J r 2 2 dr 2 |0 ls (2)| 2 



is just the part of the charge distribution inside the point r x . According 
to (9) the total contribution of this spherically symmetric charge distri- 
bution inside r x is the same value 



as if all that charge (Q x ) were localized at the nucleus. Letting 

Q 2 (r) * 47rr 2 |4> ls (r)| 2 , 

the quantity 

Q 2 (r)dr 

is the charge on the spherical shell of radius r and thickness dr. 
According to (9) the contribution of this charge to the potential is 

pQ 2 (r)dr 

and is the same for all r x inside r. 

Grunging on, .we find 

r x 2 Pi 

Q 1 =47rN 2 f r 2 2 dr 2 e" 2Cr2 = ^r f 9 dpe" p 
o Wo 

where p x = 2£r 1# Integrating by parts this becomes 

Qx-Il-e^Ml + Pi+iPi 2 )] • 
Similarly the second term of (9) is 

I 2 = 4*N Jr 2 dr 2 e" 2 ^ 2 = C^e^ 1 + e" Pl ] . 
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Thus 

J ls (1) = rT Q * + * 2 = rT (1 " e " 2 ^ 1} " ?e " 2 ^ ri ' (10) 

For large r x this becomes 

J, = -±- (11) 

Is r 2 



as expected, and for small r 2 we obtain 

J ls (D=!:-K S r 1 2 + 0(r 1 8 ) 



Thus J- has the form in Fig. 28. 
Is 




Figure 3-28. The Coulomb potential J lg (r) 
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(12) 



Using (10) in (5) we obtain 





, 3 f 1 



2C 



= 4£ { -ff 2 - 3 

(2C) («) («) 

- 8 S 



(13) 



Problem 1. Carry through the above analysis for a wavefunction 
of the form 

<t>(l)6(2) 
where 



and 



0(1) =Ne" ?r * 



0(2) =N T e" 771 * 2 



with £ and f? different. 

3.C.3 Qualitative Analysis of J^ - 

Defining the average size, F, for the 0- orbital as 
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^ls'r'ls' — 



we see from (3) that 



- 1 



An approximate value of J lg ls can be obtained by assuming each 
electron is at its average radius, r = ~-, and averaging over the distances 
between these electrons, assuming each to be on the sphere of radius r. 
If the instantaneous location of electron 1 is taken to define the z axis, 
then the average position of electron 2 will be approximately in the 
xy plane. This leads to 



r 12 = ^r 



and hence to 



1 * ?= 0.707 ? 



ls, ls ^ — y?jj 
The exact value is 



J lsls = 0.625? 
so that the above estimate is only about 10% high. 
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APPENDIX 3-D . The Exact Wavefunction of H 2 



r+ 



In order to solve for the exact wavefunction of Hf we use elliptic 
coordinates 

« = < r a + r b) /R 

-n = ( r a - r b )/R 
<p = azimuthal angle 

as defined in Appendix 2-A. With elliptic coordinates, the Hamiltonian for 
H* becomes separable (expressible as a sum of terms each depending on a 
different variable), and hence the exact wavefunction of HjJ" can be factored 
into terms each depending upon different variables, 

*U,t?,<P) = A(|)Mfo)*(?) . 

For the ground state of H 2 + at R = 2a , the resulting (unnormalized) 
wavefunction is [D.R. Bates, K.Ledsham, andA.L. Stewart, 
Phil.Tran&Roy.Soc.A246,215 (1953)] 

A (|) = (1 + |)°- 34679 (1 + 0. 01686 + 0. 00046 V 1 * 48501 ^, 
M(r)) = [1.1450P o (t?) + 0. 29844 P 2 fa) + 0.011461 P 4 (??) 

+ 0. 000184 P 6 fa) + 0. 000002 P 8 (rj)] , 
$((p) = 1, ( 5 ) 

where 6 = (? - 1)/(1 + & and P^fa) are Legendre polynomials 
[P = 1 > P 2 = i(3r ? 2 -1),...] 

For comparison the MBS wavefunction in elliptic coordinates is 



t -5r a -Cr 



^MBS = ^ e -a + e 



») 



y> R «. y> R •* 



e 



"^ * (6) 
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corresponding to 

A(0 = e" 1 - 23 * 
M(rj) = 00511(1.2377) 

at R = 2a (the optimum £ at R = 2a is 5 = 1.23). At R = 2, the MBS 
wavefunction leads to £ = 1.23, somewhat more diffuse than the 1.485 
for the exact wavefunction. The optimum energy and bond lengths for 
the exact wavefunction* are listed in Table 1 . 

Exercise : expand the M(if) expression for the MBS wavefunction in terms 
of P^ (77) and compare with the Bates et al . wavefunction. 

*When we say "exact" here we are referring to the exact solutions of 
(3) and (4). However, (3) and (4) do not lead to an exact description of H*. 
The two main assumptions here are (1) the neglect of the nuclear kinetic 
energy terms (referred to as Born-Oppenheimer breakdown) and (2) neglect 
of relativistic effects. Inclusion of nuclear kinetic energy leads to 
corrections of order — -, where M is the proton mass (in Hartree atomic 
units, e.g., 1/2M = 0. 0003 h = 0.007 eV). The actual correction to 
the energy at R e from the nuclear kinetic energy terms (see Table 1) 
is +. 00085 h = 0.023 eV = 0.53 kcal, and from the relativistic effects 
is . 000005 h = 0.00013 eV = 0.003 kcal. In order to compare with 
experiment such terms must be included (actually, for H* the experimental 
results are not yet precise enough to require these corrections). However, 
in this course we will generally ignore such effects and will refer only to 
results of nonrelativistic, fixed nuclei calculations. 



